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PREFACE 

The present volume covers those parts of the theory of the stability of rotating 
gravitating liquids that seem to be of primary importance in determining the 
evolution of such systems. Apart from the intrinsic mathematical and dynamical 
interest of the subject, the problem is also of considerable interest from a cosmo- 
gonical standpoint, as its solution is the sole source of theoretical information on 
the question of how an isolated unstable rotating mass will develop. The 
important conclusion is reached, contrary to Jeans's views and to those still 
largely current amongst astronomers, that the dynamical evidence is entirely 
adverse to the so-called fission process of formation of binary systems. The 
work accordingly removes theoretical foundation from this process as playing 
any role in the evolution of binary systems. In this way the study indirectly 
assumes astrophysical value in that in disposing of the fission hypothesis it lifts 
what seems to the writer to have been one of the major obstacles to progress 
with the greater problem of stellar evolution. 

My own interest in the subject of rotating gravitating liquids, as best I can 
remember, began twenty years ago with Professor H. F. Baker's lectures at 
Cambridge, when he used to conclude his course on celestial mechanics with a 
rapid survey of those parts of the present subject that could be dealt with by 
elementary methods. The astronomical importance of the problem became clear 
to me from studies of the origin of the solar system, and in particular it was 
from the problem of the origin of satellites that T was first led to suspect the 
validity of the fission process this, incidentally, going to show how even con- 
jectural studies can at times contribute suggestively to more precise matters. 
Investigation of the details of Jeans's researches soon disclosed numerous places 
where he had fallen into error, and J accordingly resolved to see if after rectifying 
these the subject could be brought into some more coherent form, not in conflict 
with the strong indications already available from other related fields. The 
requisite study was, however, far too protracted to be undertaken in any urgent 
or immediate way, and spread over several years it has had to take second and 
lower place to many other items that have occupied me, both of a scientific and 
non-scientific character, during much of that time. The present book is the out- 
come of this work, and while it deals for the most part only with a single line 
of attack aimed at settling the main cosmogonical question that of the fission 
process it is hoped that the discussion may serve to sum up and clarify progress 
to date and at the same time provide an introduction to the subject that may 
give at least initial help to others wishing to proceed to more recondite matters 
in which there is so much scope for further research. 

Numerous references are supplied in the Appendix, and although these do not 
claim to be comprehensive there should be no difficulty by means of them in 
tracing almost any item of the literature of the subject. 

ST. JOHN'S COLLEGE R. A. LVTTLETON 

CAMBRIDGE 

April 1951 



Chapter I 
INTRODUCTION 

The subject of the forms of relative equilibrium of a rotating mass of homogeneous 
gravitating liquid had its inception with the discussion by Newton (1687) of the 
figure of the Earth. In this it was simply assumed that a possible figure of the 
free surface would be that of an oblate spheroid with its least axis coincident 
with the axis of rotation, and it was not until Clairaut's work many years later 
that the validity of this postulate was examined. In the first instance Clairaut 
gave a proof resting on an approximate expression for the potential of a spheroid, 
but meanwhile Maclaurin (1740) produced an accurate demonstration of the 
possibility of the spheroidal form, and this led Clairaut also to publish an exact 
solution in place of his former one. Tt was rigorously shown by these authors 
that a spheroid is a possible equilibrium form whatever its eccentricity of meridian 
section provided it possesses an appropriate quantity of angular momentum. 

That an ellipsoid with three unequal axes, the least coinciding with the axis 
of rotation, is also a possible form of relative equilibrium, provided the angular 
momentum is greater than a certain amount, remained undiscovered until Jacobi 
(1834) pointed it out in a letter to the French Academy. Jacobi himself does 
not appear to have published the result, and it seems first to have been referred 
to publicly by Poisson shortly after Jacobi's communication to the Academy, 
There is perhaps something of an element of surprise about Jacobi's result in 
view of the symmetry that might be expected to be associated with any form 
produced by a rotational field, and the fact also that the Jacobian figures exist 
only if the angular momentum exceeds a certain amount no doubt contributed 
to the series being overlooked for so long. 

The first member of this Jacobiaii series is also a Maclaurin spheroid, but 
thereafter, for greater angular momentum, the equatorial axes are always different, 
and the elongation of the figure continually increases with the angular momentum. 
The limiting final form on this series, as infinite angular momentum is approached, 
has infinite longest axis, while the axis of intermediate length tends to equality 
with the third and least axis, both of them approaching zero. 

Up to this stage researches had not gone beyond the question of the existence 
of possible equilibrium forms, nor did they do so until the problem was taken 
up by Poincare (1885) in the paper that has since become celebrated. His 
investigation developed a method for studying the difficult question of the 
stability of the spheroidal and ellipsoidal forms, and necessarily involved also 
the consideration of the existence and properties of other equilibrium forms. 
There seems little doubt that a thorough discussion of these questions can be 
made only by using ellipsoidal harmonic analysis, which by the time Poincare 
commenced his work had already been extensively developed by Lame and 
others, and lay ready to hand. 



2 INTRODUCTION 

As is well known, where rotating systems are concerned, the stability of steady 
configurations of relative equilibrium is a more complicated question than for 
statical systems, and in fact two different kinds of stability may be involved, 
usually designated by the terms 'secular' and 'ordinary'. The foimer contem- 
plates the presence within the system of friction (vanishing with the relative 
velocities), whereas the latter is independent of dissipative action. Poincare* 
was able to show that if the mass of liquid is regarded as passing slowly along 
the Maclaurin series in the direction of increasing angular momentum, it becomes 
secularly unstable at a certain degree of flattening. This particular configuration, 
called a 'form of bifurcation', coincides with the beginning of the Jacobi series. 
The instability enters for a deformation of the surface involving a certain second 
order harmonic function. Moreover the Jacobi series can be regarded in its 
initial stages as the result of making a second-order harmonic deformation of 
the critical member of the Maclaurin series. This second-order deformation may 
be pictured as the addition of a stationary wave on the surface of the spheroid. 
This wave has two places of greatest elevation and two of least elevation on 
the equator of the spheroid, while the displacement vanishes at the poles, and 
it has the effect of transforming the spheroid into an ellipsoid with unequal axes. 
Thus the Jacobi series may be regarded as branching off the Maclaurin series 
and coming into existence as a result of a surface deformation involving pre- 
cisely the harmonic terms through which the Maclaurin series has become unstable. 

For greater values of the angular momentum along the Maclaurin series, 
whose members are all possible equilibrium forms, Poincare' showed that further 
separate points of onset of secular instability occur involving the higher order 
harmonic deformations. The proof of these results for the spheroids depends 
for the most part on properties of the ordinary zonal and tesseral harmonic 
functions to which the ellipsoidal harmonics reduce, in terms of suitable co- 
ordinates, when the ellipsoid has two axes equal. The corresponding investiga- 
tion for the series of Jacobi forms involves the general ellipsoidal harmonics, 
but in a closely analogous way Poincare was able to show that the Jacobi figures 
first become secularly unstable for a certain third-order harmonic deformation, 
and that for greater elongations and angular momentum, further separate con- 
figurations of instability occur, corresponding successively to harmonics of 
fourth-order, fifth-order, and so on. 

The occurrence of these instabilities accordingly suggests the possibility that 
at each point of bifurcation there branches off a new series of configurations of 
equilibrium given initially in each case by the addition of a surface deforma- 
tion expressed by the same harmonic through which the instability first appears. 
As the Jacobi series is described in the direction of increasing angular momentum, 
instability first enters for a certain third-order harmonic, and the initial forms 
of the branch series will thus be obtainable by adding a surface deformation 
of this kind. The surface wave involved in this case has three places of greatest 
elevation and three of least elevation in the equatorial plane, and vanishes at 
all points in the plane orthogonal to this defined by the two shorter axes. The 
amplitudes of the deformations supposed added are assumed infinitesimal 
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throughout, but if a small finite amplitude is adopted for the purposes of 
illustrating the nature of the new equilibrium forms, the resulting surface 
(Fig. 16, p. 110) was considered by Poincare" to bear marked resemblance to that of 
a pear. Accordingly the initial figures of the series became known as the 'pear- 
shaped' figures, or sometimes the 'piriform' figures. The difficult problem that 
accordingly next arose was to ascertain whether or not the pear-shaped series 
is initially stable or unstable. This in essence is Tchebychef's problem, and it 
gave rise to extensive investigations by Poincare*, Darwin, Liapounoff, and 
Jeans, to mention the chief authors. Minor contributions on various technical 
points connected with it were made by Schwarzschild, H. F. Baker, and others. 

The properties of the fluid mass having been assigned, as uniform, gravitating, 
and if necessary viscous, the general problem of possible equilibrium forms and 
their stability can be stated as a purely theoretical one, but Poincare' to some 
extent, and Darwin almost entirely, were interested in the problem from its 
possible cosmogonical applications. The general form of the pear-shaped figure 
undoubtedly gave rise to the notion that if the mass were stable and evolved 
by equilibrium forms along this series, with the furrow continually deepening 
as the figure elongated, the final result would be two detached masses rotating 
in circular orbital motion about each other. Thus it seemed plain to Darwin 
that the dynamical theory, if it could be established in accordance with these 
ideas, would give a strong theoretical basis for supposing this to be the method 
of genesis of double systems in the celestial universe. Indeed, Darwin eventually 
announced that he had proved the pear-shaped series to be initially stable and 
hence that this would in fact be the course of development. 

Now in order to establish that the pear-shaped figures are stable it is sufficient 
to show that they are secularly stable. On the other hand, if the series were 
shown to be secularly unstable, it would in general require further investigation 
to decide in what manner the system is next likely to develop. At all events, 
Darwin approached the problem by a method that aimed only at settling the 
question of secular stability, as also did Liapounoff, though this author appears 
to have been attracted solely by the theoretical problem, and in no wise inter- 
ested in the astronomical implications of the results, which pass unmentioned 
in his numerous papers. Subsequently, Jeans also concerned himself with 
establishing the secular stability or otherwise, apparently under the impression 
that ordinary stability was not relevant to the problem. 

Had these writers been able finally to conclude that the pear-shaped figure 
was secularly stable, their treatments would have represented a complete 
solution of the immediate problem (though there would have still remained the 
question of how far the pear-shaped series continued to be stable with increasing 
angular momentum, just as this was necessary for the Maclaurin and Jacobi 
series). But in point of fact Liapounoff, and later Jeans, concluded that the 
figure was quite certainly secularly unstable, and moreover Jeans claimed that 
Darwin's original investigation, when properly interpreted with certain technical 
errors set right, itself in reality led to this same result. As far as the objective 
of Liapounoff 's inquiry was concerned he had attained it, but Jeans's expressed 
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intention was to obtain theoretical evidence of the course of dynamical stellar 
evolution, albeit under what have since proved perhaps unduly idealized condi- 
tions. The desired information, however, is not conclusively provided by the 
knowledge that the system is secularly unstable, for this means that the pear- 
shaped form itself never comes into existence, and if the continuation of the 
Jacobi series remained ordinarily stable, as might possibly happen, the rate of 
departure from the critical Jacobi figure might not take place at all rapidly if 
frictional effects were small. For example, the lunar orbit is secularly unstable, 
but ordinarily stable, and its evolution under frictional forces proceeds extremely 
slowly and would cease altogether if friction were absent. 

Thus to complete the information that may be derived from the consideration 
of small deformations of the system it is necessary to determine whether the 
Jacobi series remains ordinarily stable or not beyond the critical Jacobi figure. 
It is automatically so before this stage is reached as a consequence of its secular 
stability. This question requires a totally different treatment from that of 
secular stability, since it is necessary for its solution to study the actual periods 
of possible small oscillations of the system, and not merely the manner in which 
some single entity, such as the moment of momentum, changes along the initial 
stages of the pear-shaped series. The determination of the ordinary stability 
of the Jacobi series has been undertaken and solved by Cartan, who has 
succeeded in proving that for displacements involving the third-order harmonic 
through which the ellipsoids first become secularly unstable, they simultaneously 
become ordinarily unstable. 

With this information available the nature of the development beyond the 
critical Jacobi figure can be studied with more certainty. For if the free surface 
receives a displacement involving third-order harmonics, and any general physical 
disturbance may be assumed to contain such terms in its expression, the ampli- 
tudes will begin to increase exponentially at a rate independent of the amount 
of friction present. The system can no longer oscillate about an equilibrium 
form, since none exists of a stable character, and instead a dynamical motion 
must ensue until the system succeeds in finding its way to a new state of steady 
motion. The equations of small motion of the system permit this development 
to be followed only so long as the velocities and displacements involved in it 
remain small, but with increasing amplitude the approximations leading to 
linear equations of motion become less accurate. The system must, however, 
reach eventually some other steady condition involving no further dissipation 
of energy, and the interesting question arises as to what this final configura- 
tion will be. Unfortunately it is not possible to investigate the question in 
detail by anything approaching rigorous means, but it may well be, as has 
always been maintained, that the result will be a division of the original mass 
into two detached portions. However, if this view is correct, there is necessarily 
an important difference from Darwin's ideas on the course of development, for 
it can be shown not only that the pieces must be of considerably different sizes, 
but what is more important still, that they ,must separate to infinity. The final 
steady state would then consist of two separate unequal stable masses receding 
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with constant relative velocity, the original excess angular momentum causing 
the instability reappearing now as orbital angular momentum. 

From the point of view of cosmogony the main question of interest is to obtain 
as rigorous a demonstration of this course of development as possible. It would 
also be of interest to give as complete an account as possible of the whole evolu- 
tion of the problem, but the literature of the subject is so extensive and much 
of it of an exploratory character, that it would scarcely be practicable in a single 
volume to give more than an outline of it. Yet it has seemed worthwhile to 
give the full mathematical discussion of such parts as are essential to establishing 
the extreme plausibility of the course of evolution described above. To do this 
we begin by discussing the subject of stability with particular reference to 
rotating systems. This is followed by a discussion of the spherical, spheroidal, 
and ellipsoidal forms, together with certain of their properties that can be 
established by simple means, as illustrations of the dynamical theory. Next is 
developed the ellipsoidal harmonic analysis required for further progress, when 
the necessary properties of Lamp's functions are derived, and then using this 
mathematical technique an account is given of Poincare*'s investigation of the 
secular stability of the Maclaurin and Jacob! series. This is followed by an 
account of Cartan's discussion of the ordinary stability of the Jacobi forms. 
Finally, the subsequent development of the system is considered arid its possible 
cosmogonical implications discussed. 



Chapter II 
STABILITY 

STABILITY OF STATICAL SYSTEMS 

Equilibrium configurations 

Let us consider a mechanical system whose position can be specified by n 
generalized coordinates q-^q^y > ( ln> an< ^ ^ s motion at any instant by the general- 
ized velocities q^q^ ...,</ n , where the dots denote differentiation with respect to 
time t. Also let us suppose that the forces acting on the system are derivable 
from a potential energy function V depending on the coordinates only, so that 

V = V(q t ), 

The kinetic energy T will in general depend on both the coordinates and the 
velocities, and will be a homogeneous quadratic function of the latter. Thus 

T = r(g,,4,). 
Motion of the system will take place in accordance with the Lagrangian equations 



which possess the energy integral 

T+ V = constant. (2) 

Hence the possible equilibrium configurations, within the limits of the number of 
coordinates adopted, are determined by the n equations 

^ = (-l,2,...,n), (3) 

which are simply the conditions that V is stationary. 

These equations may have certain admissible solutions, each of which will 
correspond to a possible equilibrium state, and we may denote a particular one by 

q i = a i (i = l,2,...,n). (4) 

Moreover, it can readily be seen that if V is an absolute minimum in the con- 
figuration then the equilibrium is stable. For, in accordance with the energy 
integral (2), if a slight disturbance from equilibrium occurs, it follows, since T 
cannot become negative, that V cannot increase above its equilibrium value by 
more than a very small amount. This means in turn that none of the coordinates 
can deviate by more than a small amount from its equilibrium value, and there- 
fore that during the motion the system must always remain in the immediate 
neighbourhood of the equilibrium configuration, which is what is meant by 
stability. 
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Linear series of configurations 

If the description of the system contains a parameter p, say, not itself dependent 
on any of the generalized coordinates, which for any reason is slowly changing, 
or is so regarded, the potential V will in general depend on /z, and the values of 
g t - giving the solutions of (3) will also depend on p. Thus the particular solution 

(4) may be written / v / , x 

V y 0< = a<(fO ( = 1,2,. ..,%), 

where the a/s are now functions of p. 

When IJL undergoes a small change dp, the solution will become 

q i = a,(/A + d/i) = a i + -.d f jL 9 (5) 



and hence the original configuration will give place to an adjacent configuration 
for the new system with slightly different JJL. Thus starting from a given 
equilibrium position, a continuous series of other possible configurations is 
obtained as /LI slowly varies. Such a set of configurations Poincare has termed 
a 'linear series'. The essence of the idea is the presence within the system of a 
changing parameter which, however, varies sufficiently gradually for the change 
not to affect equilibrium at any stage. 

Stability 

Let us consider next what may happen to the stability of the system as it 
moves slowly along a linear series. We examine this question first by diagram- 
matic means. If we take (q^ /z) as rectangular coordinates with the /i-axis 
vertical, then for different values of V on the right-hand side, the relation 

V(g i ;n)=V (6) 

will correspond in two dimensions to a family of curves, or in three dimensions 
to surfaces. With more than two coordinates, q t , the relation will represent 
hyper-surfaces according to the number of dimensions involved. As there is 
only one value of V for a given point (q^^) these surfaces are essentially non- 
intersecting. 

The condition that the tangent plane to a given F-surface shall be per- 
pendicular to the /z-axis is simply that dV = for dp = 0, which is equivalent 
to the relations (3). The equilibrium configurations therefore correspond to the 
points where the F-surfaces are horizontal. In Fig. 1, let V l and V 2 represent 
two surfaces of the family and suppose V 2 > V v Then the point P represents an 
equilibrium configuration and the corresponding value of /z is p l9 as shown. The 
heavy line joining successive equilibrium points P lt P 2 , . . . represents the linear series. 

Slight displacements of the system, without change of /x, will be represented, 
for /LI = /Ltj say, by points such as P/ in the tangent plane at P r In the circum- 
stances postulated, such changes will involve an increase of potential energy if 
the F-surfaces are concave downwards as shown, since the general direction of 
V increasing is as indicated. Thus the condition for stability can be stated in 
the form that the F-surfaces must be concave downwards. 
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Clearly, if the general direction of increasing V were in the present case the 
other way, stability would require the surfaces to be concave upwards. The 
rule is therefore that for stability the concavities must be towards the direction 
of V decreasing. 



/u-axis 



Direction of 
gradual change 
of H 




V increasing 



Exchange of stabilities 

As the representative point of the configuration passes along the linear series 
no new feature arises so long as the F-surfaces remain concave in the same 
direction. However, it may happen that they gradually develop in one of the 
ways shown in the diagrams of Fig. 2, which illustrate some of the possible 
forms that the F-surfaces may take. The analytical discussion given later covers 
all possible cases. At present we consider these four examples. 

(i) As the linear series is described in the direction indicated by the arrow, 
the concavities eventually become the other way and additional downward 
concavities appear. Accordingly, at the critical point C another linear series 
BOB cuts across the original series A 1 CA 2 . Such a point is termed a 'point of 
bifurcation', and the corresponding equilibrium position a 'form of bifurcation'. 
If the series A^C was at first stable, the change in direction of the concavities 
means that beyond C the continuation CA 2 of this series must be unstable. 
That is, at C the original series loses its stability. On the other hand, in the 
case depicted, the new series BCB has its concavities the same way as the series 
A^C, and hence the configurations of this series will be stable. There has occurred 
a transfer of stability to the members of the new series. For each value of JJL 
corresponding to the portion A C t there exists only one equilibrium form, and 
this is stable, but for each value beyond this there are three possible forms of 
which two are stable and one unstable. 

(ii) Here, assuming the initial series A^C to be stable, the continuation CA 2 
of this series must again be unstable, but the series BCB must now also be 
unstable. In this case there is therefore a disappearance of stability at C. 

(iii) Here two configurations of equilibrium are possible to begin with for each 
value of /x, and if that on the series A^C is stable, that on CA 2 is unstable. Not 
only is stability lost at (7, but for values of ft beyond that of C there are no 
equilibrium forms possible at all. 
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(iv) Here only one configuration is possible for values of p corresponding to 
the series A^C, and again there are no equilibrium forms possible beyond that 
corresponding to C. 

Cases (i) and (ii) are of special importance in the actual problems with which 
we shall be concerned. In the application of these considerations to physical 
systems it is convenient always to choose for //, some parameter that increases, 




V increasing 




V increasing 




V increasing 




V increasing 



0") 



or decreases, monotonically as the system gradually evolves. We shall see that 
one or other of such quantities as the angular velocity or angular momentum 
can usually be postulated to satisfy this requirement. Assuming such a choice 
to have been made, it is seen that in (ii) the members of the unstable series 
BOB cannot actually come into existence, even though the value of /A would 
permit it theoretically, if the system is originally in a stable configuration on 
A^C. It could, of course, if necessary, be regarded as placed in one of these 
unstable positions. 
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It may be remarked that the term 'exchange of stabilities', which is often 
used in connexion with cases such as (i) and (ii), is to some extent a misleading 
one, because stability is not necessarily transferred to the new series. There is 
always a loss of stability of the original series, but the new series, if one exists, 
may or may not be stable. 

Following Jeans, these possible cases can be represented schematically by the 
following idealized diagrams: 




0") ( )V > 



Fig. 3. The heavy lines denote stable series and the dotted lines unstable series. The 
arrows indieate the direction of evolution. 

Rules for graphical determination of stability 

These considerations therefore lead to the following rules for determining 
diagrammatically how the stability of a system changes at a point of bifurcation: 

If in evolving along a stable series represented by a vertical line a point of 
bifurcation is reached, stability will be lost by the first series and transferred 
to the new series if the curve representing the latter turns upwards. This is 
case (i) above. 

If the new series turns downwards, as in (ii), configurations on it cannot come 
into existence, and the continuation of the original series is unstable. 

If the original series itself turns downwards, as in (iii), or ceases altogether, 
as in (iv), there are no further equilibrium forms either stable or unstable for 
values of the parameter immediately beyond that corresponding to the point 
of bifurcation. 

It is highly important to understand that the number of coordinates used to 
describe the system must remain the same throughout for these results to be 
valid. Continuous systems, such as liquid masses, require for their complete 
description an infinite number of coordinates, but by restricting the mass to 
special forms the number of coordinates required can be made finite. Thus, if 
only ellipsoidal forms are admitted, then but two coordinates are needed, for if 
a, b t c denote the semi-axes these must always be related by abc = constant. 
If only spheroidal forms are admitted, then a = b also, and one coordinate 
suffices. If a diagram of the above kind is drawn for the Maclaurin series no 
point of bifurcation appears because no other equilibrium series is possible 
within these restrictions. On the other hand, if a diagram is drawn for the 
ellipsoidal series, it is found to have a point of bifurcation at the point where 
the Maclaurin series crosses it. 

It is for the present reason that the possibility of the pear-shaped figure only 
became demonstrated when the full number of freedoms could be coped with 
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by the method of ellipsoidal harmonics, though something of the kind had already 
been vaguely conjectured by Kelvin and Darwin. 

Analytical discussion 

The condition for a point of bifurcation on a linear series can readily be deter- 
mined analytically. 

Supposing q i = a t is an equilibrium configuration when the value of the para- 

8V 
meter is /*, so that = for this configuration, let us consider the value of the 

<%i 
potential energy F^J/LI) in an adjacent configuration given by (aj + St/ 

As far as second-order terms in the small quantities 8q it fyt, we have 



(7) 



where summation is made over repeated suffixes. 

For a fixed value of 8/1, the quantities 8q t may be regarded as the coordinates 
of the system, and the conditions that (a^ + S^j/x-f S^u,) represents a new 
equilibrium configuration are the n equations 



8V 



= (i=l,2,...,n). 



Writing V rs for 8 2 V/8q r 8q K and V r/t for 8 2 V/dq r 8p,, these conditions, written at 
length, have the form: 



=0 



The solution of these may be written in dcterminantal form as follows: 



where 



A = 



v nl 



I V V 

Ml M 2 



These equations determine the requisite step 8^ (in a space of n + 1 dimensions) 
from a known equilibrium configuration a t to a new equilibrium configuration 
for a given increment 8/z, of the parameter /*. They will fail to determine such a 
step only if 



= 0, 



(10) 
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and this is the analytical condition for a point of bifurcation or a point where a 
series ceases with increasing p. 

Returning to the diagrams of Fig. 2 for a moment, we see that in case (i) if 
the original configuration is just below (7, there are three possible adjacent 
equilibrium configurations for a small increment S/i, so that for such a system 
A must be zero at (7, as otherwise a unique adjacent configuration would neces- 
sarily exist. Similarly in case (ii) A must vanish at C. In cases (iii) and (iv) the 
linear series are such that Sfjt, must itself be zero if the configuration is to remain 
on the series, and hence in these cases too A must be zero at C. 

Whenever this happens the next stage is indeterminate to the present order 
of approximation and could therefore be found only by proceeding to a higher 
order, when evidently the equations corresponding to (8) would no longer be 
linear, and their solution would give more than one set of values for the 8q it 
These solutions would presumably be real for positive S^u, in a case such as (i), 
and for a negative SJJL in a case such as (ii), whereas for (iii) and (iv) their solution 
would be imaginary and therefore not correspond to an admissible step to a 
real adjacent configuration. 

Stability 

Let us examine next how the vanishing of A affects the stability of the original 
series of configurations (/ t = a L . 

The parameter \L is now to be kept constant, so that the variable part, &V 
say, of the potential energy for adjacent configurations is 

8V = Wrs^r^s- 

By means of a real non-singular linear transformation of the coordinates, now 
8q t , this expression can be transformed to a sum of square terms only, and 
denoting by 1 ,0 2 ,...,0 n a new set of coordinates chosen in this way, we shall 
have 



wherein 6 1 ,6 2 , ...,& are certain constants depending on the original position of 
the system and on /*. Let us denote the determinant of this transformation by 
m, so that ra will be essentially non-zero. By a known theorem, the discriminant 
of the new quadratic form will be simply m times that of the original form. 
Thus we must have 
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which reduces to 



6 1 6 2 6 3 ...6 re = 



(11) 
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For the original configuration to be stable, the potential energy 8V must be 
positive definite, and, by a known theorem, if this is so it means that b lt 6 2 , ..., b n 
must all be positive. Whatever their values, positive or negative, the quantities 
6 1 ,6 8 > ...,6 n are called the 'coefficients of stability' of the system. They are not 
unique, since a quadratic form can be transformed to square terms only in an 
infinite number of ways; but however the reduction is made, the numbers of 
positive and negative coefficients in each of the resulting expressions are 
invariable. 

If the original system is stable, 6 l5 6 2 , ...,6 n will all be positive initially, and a 
change from stability to instability will occur, as \L changes, only when one of 
these coefficients vanishes and changes sign. The values of p, for which this 
occurs are accordingly given by the same equation as that determining the 
points of bifurcation, namely, A 

Tt may therefore be concluded that at a point of bifurcation there occurs a 
loss or gain of stability. If the system is describing a part of a linear series that 
is stable, it must lose its stability for the configurations on the continuation of 
this series. If a branch series passes through the point for which A = 0, it will be 
stable if its representative curve turns upwards so that new equilibrium con- 
figurations on it become possible as the parameter /z continues to increase. If it 
turns downwards these configurations will be unstable and cannot come into 
existence. 

We give next some simple examples of statical systems illustrating the fore- 
going theory. 

EXAMPLES OF STABILITY OF EQUILIBRIUM OF STATICAL SYSTEMS 

(i) A heavy uniform rod AB of mass M and length 2a is supported by two equal 
crossed strings BC.AD of lengths 2b attached to two points C and D at distance 2a 
apart at the same level, and the system is 
constrained to move in the vertical plane 
through CD. 

Denoting by 26 the inclination of the rod 
to the horizontal, consideration of the depth 
of the mid-point of AB below the level CD 
gives immediately 

U = V/2Mg = - cos (6 a - a 2 cos 2 0) 
where g is the acceleration of gravity. Equilibrium positions are therefore 

S iven by rfC7 _ sin (fc 2 - 2a 2 cos 2 0) 

~dB~-~ (& 2 -a 2 ^os 2 0)* " ' 

that is, sin = 0, or cos 2 = b 2 /2a 2 . 

The relevant solutions are 

= 0, or = cos~ l (b/aj2). 
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The first corresponds to the horizontal position of the rod and is always a possible 
equilibrium configuration. The second, when it is a real angle, gives two sym- 
metrical positions in which the rod is inclined to the horizontal. 
To examine the stability of these various positions, we have 

d* U cos0(6 2 -2a 2 cos 2 0) a 2 sm 2 0cos0(36 2 - 2a 2 cos 2 6) 
' 



dB 2 (6 2 -a 2 cos 2 0)* (6 2 - a 2 cos 2 0) 3/2 

If = 0, this gives U ee = (6 2 - 2a 2 )/(6 2 - a 2 )* and the position is therefore stable 
if 6 > a ^2, and unstable if b < a J2. Thus the second inclined positions only exist 
when the horizontal position is unstable. If cos 2 = 6 2 /2a 2 , we find 

~2a~ z , 

and since this is essentially positive, these inclined positions, when they exist, 
are always stable. 

If, therefore, the system is imagined to begin with a/b small, and gradually 
increasing, the series of equilibrium positions will provide an instance of case (i) 
of the foregoing theory. The horizontal position is the only equilibrium form 
possible, and is stable, until a ^2 becomes greater than 6, when two additional 
stable equilibrium forms become possible, while the horizontal position becomes 
unstable. 

(ii) A heavy bead of mass m can slide on a smooth vertical circular wire, radius a, 
and is attached to a second particle of mass M by an inextensible string that passes 
over a smooth peg at the highest point of the circle, the mass M hanging freely. 

Denoting by the angle between the two portions of the string, as shown, 
then we find 

U = V/2ga = M cos - m cos 2 0. 
Equilibrium positions arc therefore given by 

- = M sin 4- 2m cos sin = 0, 
du 

that is, = 0, or = cos~ l M/2m, 

which gives real only if M < 2m. 
To consider stability, we have 




= -Jf cos0+ 2m(2cos 2 0-l). 
u* 

Fig. 5. 

If = 0, U ee M -f 2m and is positive if M < 2m. 

If cos = M/2m, U ee = (M 2 - 4w 2 )/2w and is negative if M < 2m. 

If therefore M/m is taken as parameter and regarded as gradually increasing 
from zero, there are at first three possible equilibrium positions of which = 
is stable and the two symmetrically placed positions on either side of the vertical 
are unstable. For M/m > 2, only the former position is possible and it is then 
unstable. The system is therefore an example of case (ii) of the theory. 
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(iii) A heavy bead can slide under gravity on a smooth wire in the form of a conic- 
section with major -axis horizontal and minor -axis vertical. 

For the parameter let us take the eccentricity e of the curve, and suppose 
that the latus rectum 26 2 /a remains fixed. Then for e< 1, equilibrium positions 
are at the ends of the minor axis. The lowest position is stable and the highest 
unstable. For e> 1, however, the curve becomes a hyperbola and no equilibrium 
positions exist. We have therefore an instance of case (iii). 

(iv) A heavy particle can slide under gravity on a smooth 
circular wire in a vertical plane and is attached by means 
of an elastic string to a point A of the circumference at 
angular distance OL from the highest point. 

Supposing OL is initially zero and the string sufficiently 
short for the particle to be in contact with the curve, 
then if a is gradually increased, so that the particle is 
drawn up the curve, there will be only one possible equi- 
librium position for each value of a, and this position 
will eventually reach the highest point. If any further 
increase in a occurs, the system will cease to be a 
statical one, so that the linear series will end abruptly. 




STABILITY OF ROTATING SYSTEMS 

Where rotating systems are concerned, the configurations of importance are 
those of relative equilibrium in which the whole system rotates steadily about a 
fixed axis through the centre of mass as if it were a rigid body. In such a state 
there are no relative motions of the parts, so that as far as internal influences 
are concerned no dissipation of energy can occur and the system remains steady. 
If relative motions do occur for any reason there will in general no longer be a 
unique angular velocity of the system, though the direction of the angular 
momentum vector may provide a fixed direction at the centre of mass. If so, we 
can adopt a set of rectangular rotating axes at this point of which the third, say, 
is fixed and the other two rotate round it. The positions of the particles can 
then be referred to this rotating frame. 

Let us suppose the axes Ox,Oy,Oz to be rectangular and rotating about Oz, 

defined in this way, with angular velocity o> = 7\ not necessarily constant. 

If a particle of mass m has coordinates (#, y, z) in the rotating frame, its space 
velocity at any instant will be (x a>y,y + wx y z). Hence the kinetic energy of a 
system of such particles is given by 

T = $Zm{(x - a>y) 2 + (y + a>x) 2 + z 2 } 

= %Zm(x* + y 2 + z 2 ) + <*>Zm(xy - yx) + \u>*I*m(x 2 + y 2 ). 

The angular momentum round the z-axis relative to a fixed frame of reference 
is given by % = % m {(y + a)X ) x -( x - w y) y } 

yx) + wl,m(x 2 + y 2 ). 
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If we write 

T R = iX 



then T H is the kinetic energy relative to the rotating frame, and H R is the angular 
momentum relative to this frame. Also / is the moment of inertia of the system 
round the z-axis, and is independent of the frame in which the rectangular 
coordinates x, y happen to be measured. Expressing T and H in terms of these 

quantities we have 

(12) 



(13) 
and, by eliminating H Rt we have also 

T = T R + a>H-w 2 I. (14) 



Let it be supposed next that the system is capable of being described by 
n+ 1 coordinates, consisting of n generalized coordinates q i that fix the configura- 
tion relative to the rotating frame, and the azimuthal coordinate 0, which will 
be assumed not to occur explicitly but only as if/. Assuming, as will be the case 
for natural systems, that the internal forces are derivable from a potential 
function F, and that the only external force is a couple G round the z-axis, the 
Lagrangian equations of motion for the system are 

~ -' ..... > 



-*- 

BT dT 

But from (12) ---.- = arid ,. = H, so that (16) becomes 
dip doj 

%-* 

which expresses the equality of the rate of change of angular momentum round 
Oz with the external couple. 

For a freely rotating system G will be zero and // constant. Tf, on the other 
hand, the system can always be described by means of n generalized coordinates 
relative to axes constrained to rotate uniformly with angular velocity co, the 
time not occurring explicitly in the description, an external couple will be required 
of amount given by (16). For our purposes there are two main cases to con- 
sider. First, we shall take aj = constant, and second, the case // = constant. 
We shall also show in what circumstances the two cases may be regarded as 
equivalent. 
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Systems rotating with constant angular velocity 

To transform the Lagrangian equations (15), we treat the <//s and <//s as 
independent in the usual way, so that we have 

. _ &r . dx _ fa dx d dx 

X ~~ (Ji ' " '' 



mu 

Then 

^ *. 
so that 

, 
and 



R R ^ - - 

Hence , (- ^ - -- -^ = 22m # O l/ -w _- 

frji \ 8< h y 



Tf R V ( &/ 8X\ V ( ty fa\ 

- , R - = S?n \x - J . -y -I = Lra (x ~--y Q - , 
%<h \ ?><h <V \ %ii <W 

fy . to 
^-7/ a 

q t J Bq 



summed also over j. (18) 



CiC ?/ 

The factors _. ^-r are real geometrical quantities associated with the system 

%y. ^) 

and can be regarded as of known analytical form when the generalized coordinates 
have been selected. We may denote them by p ijy so that we have 

ft, d9) 

^ 3l ' 

and & = 0, for each value of i, no summation being involved. 

Substituting the value of T given by (12) in the left-hand side of equations 
(15), we have 



The term on the left multiplying w has just been calculated, and inserting this 
in terms of the & ; 's, the equations become 



These are the equations of motion of the system relative to uniformly rotating 
axes. They differ in two important respects from the usual equations referred 
to unaccelerated axes. First, by the presence of the term io> 2 /, which appears 
added to the potential energy. The force terms resulting from this are sometimes 
called the ordinary centrifugal forces. Second, by the terms in the velocities, 
w Pij ( li> which depend also on the rotation w; these are usually called gyroscopic 
terms, but sometimes also the compounded centrifugal forces. Their presence is 
of special importance, for they can, if sufficiently strong, render a system stable 
in a certain sense (to be explained later) when the usual energy conditions would 
indicate instability. 
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The terms F t have been added in (21) to represent any external forces that 
might be acting. It will be convenient later in discussing the possible effects of 
friction to identify these terms with the frictional forces. 

Conditions for relative equilibrium 

If a configuration is one of equilibrium relative to the uniformly rotating 
frame, then we shall have T R = and all the velocities q t = 0. Hence in the 
absence of any external forces the conditions are simply 

^-(F-io,2/) = (=l,2,...,n). (22) 

% 

Thus the only difference in the conditions for equilibrium, as compared with 
non -rota ting statical systems, is the appearance of the term Jo> 2 / added to the 
potential energy. That is, to find the equilibrium position we simply use 
V Jeo 2 / instead of V and express the condition that it shall be stationary. We 
may term this function the total mechanical potential. 

Conditions for stability 

If we multiply equations (21) by q t and sum for all i we obtain 



Since j8y = )8^, the terms in a* vanish identically. To reduce this equation 
further, we have, by Euler's theorem, since T R is a homogeneous quadratic 
function of the <y/s, 



and differentiation of this with respect to the time gives 

dT R _8T R 
~ ~ 



Also by direct total differentiation of T R we have 



.,., 

Whence 



dt - 
dT d 



so that (23) becomes - , (T R + V - Jco 2 /) = F i q i . (24) 

dt 

If there are no external forces, other than the couple maintaining the constant 
rotation, which is not of course included in the jPJ's, this equation possesses 
the integral ^ + y _ ^ = oonstant (26) 
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Secular stability 

For equilibrium we have seen that V Jco 2 / must be stationary, and this value 
may always be taken as zero. If, further, the configuration is such that F-^o> 2 / 
is an absolute minimum, the position is said to be secularly stable. 

Let us consider now small motion in the neighbourhood of such a position and 
suppose that the coordinates are chosen, as is always possible, to vanish for the 
equilibrium form. They can also be transformed to make both T R and V \<*PI 
sums of square terms only. This is always possible by means of a real linear 
non-singular transformation since T R is necessarily positive definite. We may 
therefore write, using q i now to denote the new coordinates, 

2T R = a^gf, where a t > 0, necessarily, "I 

and 2( V - -oj 2 /) == btf*, where b t > in this case.) 

If the system is given a slight displacement, the integral (25) will become 

', (27) 



where c is a small positive constant. Since T R cannot become negative, this 
means that each coordinate can never exceed a certain amount. For instance, 
</! can certainly never exceed (2c/6 x )* and may not reach this value. Thus, in 
the absence of friction, the system must oscillate in the immediate neighbourhood 
of the equilibrium position, the amplitudes of the various terms expressing the 
oscillation all being small. The system will therefore be stable in the ordinary 
sense of the word used for statical systems. Such (rotating) systems are in fact 
said to be ordinarily stable. 

If friction proportional to the relative velocities is present, as it must always 
be to some extent in natural systems, these oscillations will die away. For the 
frictional forces may be regarded as external forces acting in such a way that 
the terms F^q^V^q^ ...,F n q n are necessarily all negative, since friction always 
operates to decrease the kinetic energy. Equation (24) then implies that 

(T R + F ^co 2 /) = a negative quantity. 
at 

Hence, however slight the friction, T R +V \w 2 ! must continually decrease so 
long as any of the velocities q t are different from zero. Thus both T R and 
V |co 2 / must become zero, which means that all the coordinates fall to zero, 
and the system returns to the equilibrium position. It is not sufficient for T R 
alone to fall to zero, since until V Jo> 2 / also becomes zero the system cannot be 
permanently in equilibrium and the velocities would not remain zero. Such a 
system, in which all the coefficients of stability are positive, is said to be secularly 
stable. We therefore see that if an equilibrium configuration of a rotating system 
is secularly stable it is also necessarily ordinarily stable. Such systems are some- 
times said to be 'thoroughly stable'. 

The difference between the two kinds of stability will become more clear from 
the subsequent discussion of unstable systems, for which the distinction is of 
great importance. 
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Ordinary stability of rotating systems: small motion 

Let us consider next the equations of small motion in the neighbourhood of an 
equilibrium configuration of the rotating system. Using the quadratic forms (26) 
for T R and V |<o 2 /, the equations of motion (21) become, for small values of 
the </^'s 

+ V/2 + "(021 '/I + 023</3+ ) = -. /0 , 

(28) 



For oj = these equations reduce to the usual form for principal or normal co- 
ordinates. The free oscillations, F l = 0, of any one of the coordinates is then 
independent of that of any of the others. But for rotating systems this no longer 
holds even though both T n and V iw 2 / are reduced to the sums of square terms 
only. For if any one coordinate, say q lt varies, the presence of the gyroscopic 
terms in general brings about consequent changes in all the remaining coordinates. 
Investigation of the free oscillations, or small motion, can however be made 
in the usual way; that is, put F t = and assume a solution of (28) of the form 

<h = <ho* xt (i=l,2,...,n), (29) 

where the </ io 's are small constants, to be determined, and the exponential term 
is the same for every coordinate. This will be a solution provided the following 
system of linear equations is satisfied: 



0, 

0, 



(30) 



The expression of the consistency of these linear homogeneous equations in the 
(7io' s gives rise to a deterrninantal equation for the possible values of A, namely 



(31) 



>j8 21 A a. A 2 4- 6 
P21 " 



The small motion of the system accordingly depends on the roots of this equation. 
If A is one of the roots, and in general the roots will be non-multiple, the n equa- 
tions (30) are then consistent, so that any nl of them can be solved for the 
ratios of the constants <? 10 :</2o : : (7/o- ^ ne a ^ least of these constants must 
remain arbitrary, but the remainder will be expressible as multiples of it, the 
multiples being different for the different values of A. The complete solution of 
(28), with F i = 0, will then be given by the linear sum of the 2n solutions obtained 
in this way, one for each of the 2n roots of the equation (31) in A. 
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Stability 

If all the roots in A are not pure imaginary, terms increasing exponentially 
with the time will occur for a general small displacement of the system. To show 
this, we notice first that (31) is in effect an equation in A 2 . For since J; - = j8 j , 
if A is written for A, the rows and columns are interchanged. Had w not been 
present, only the terms in the leading diagonal would have occurred, and the 
roots would be simply 

A 2 = -Vi> -AK > -&>, (32) 

so that stability would require b l9 b 2 , ...,& to be all positive. But now the deter- 
minantal equation no longer has simple roots like this, and in general roots in 
A 2 other than purely real ones will occur in pairs of the form 



and these will give individual values of A of the form 

A = (p iq). 

The time factors in the coordinates, when put in real form, with these values of 

A will be of the form A .. , M . n . t . \ 
A e* 1 * cos (qt + a) + B e-' jt cos (qt + jS), 

and a general oscillation, as opposed to one specially started, will contain all 
these terms. Accordingly, if p 4= 0, one or other of the amplitudes A e lft or 
Bcr* tl will increase indefinitely. The motion will not remain small and the system 
will be ordinarily unstable. 

The, condition for the system to be ordinarily stable is therefore, that, regarded as 
an equation in A 2 , the determinantal equation (31) must have, all its roots real and 
negative. 

Stated in this form, the condition is the same as for non-rotating systems, as 
expressed by (32), but it is particularly to be noticed that for rotating systems 
this no longer requires all the coefficients of stability b lt b 2 , ...,/> w to be positive*. 
For if Af, A 2 ,, ..., A 2 are the (supposed) real and negative roots of the equa- 
tion in A 2 , then by considering the term of highest degree and the constant term 

in (SI) w. got 6lV .A = 1 <V..,,A?Ai...A*. (33) 

Hence, since a 1} a 2 , ..., w are all essentially positive, the sign of the product 
Af A|... A 2 , which must be positive for stability, will be the same as that of the 
product &! b 2 . . . b n . Thus, as far as this requirement is concerned, ordinary stability 
may hold even though some of the fc/s are negative, but the number of negative 
coefficients of stability must nevertheless be even. We thus have the following 
important result (for systems involving a finite number of degrees of freedom): 

// a finite system becomes secularly unstable as a result of the changing of sign 
of one (or an odd number) of the coefficients of stability, it must simultaneously 
become, ordinarily unstable. 

The foregoing discussion does not prove that a system can be ordinarily stable 
when V ^aj 2 I is a maximum with regard to some of the coordinates, but only 
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that negative coefficients of stability, if even in number, are not inconsistent 
with ordinary stability. To establish the result the roots of (31) would have to 
be more fully investigated, and this would require detailed knowledge or assump- 
tion as to the values of the /? t -/s and w. But it is not necessary for our purpose 
to consider the question with full generality; instead we show below, for the 
case of two degrees of freedom, that a system can be ordinarily stable when 
V o> 2 / is a maximum. 

Systems of an infinite number of degrees of freedom 

The present discussion concerns systems of a finite number of degrees of free- 
dom, whereas for a liquid the number of freedoms is infinite. It has been shown 
by Hilbert that the theory can, however, be rigorously extended to infinite 
systems, and the occurrence of imaginary or complex solutions, so that terms of 
the form e p/ appear in the coordinates, can again be regarded as indicating 
instability as for finite systems. The equation for A will in general have an 
infinite number of roots and only in special cases will it be representable in 
algebraic form, but it so happens that in the problem with which we shall be 
concerned the period equation when suitably derived breaks up into an infinite 
series of equations each of algebraic form. 

It is important to remember, however, that results derived above, such as that 
concerning the number of coefficients of stability that change sign and its relation 
to stability, will not necessarily carry over unchanged for infinite systems. To 
avoid any possibility of error arising in such a way, the safer plan is to examine 
each particular problem on its own rather than to infer its properties from results 
established only for finite systems. 

Condition for loss or gain of ordinary stability 

Suppose we have a configuration of a rotating system that is secularly stable, 
and therefore also ordinarily stable, let us consider what may happen as the 
system evolves along a linear series. In general all the roots Af, A|, ..., A^ 
of (31) will be different, so that if we regard them as arranged in order of nrngni- 
tude, that is, Af < Af< ... < AJ, the first root that may change sign will be Af. 
But this is equivalent to the product 6 1 6 2 ...6 n changing sign. Hence ordinary 
stability is certainly lost when one, or an odd number, of the coefficients of 
stability changes sign, whereas it may continue to hold if an even number of 
coefficients of stability change sign simultaneously. 

Secular stability cannot hold after ordinary stability is lost, but it may hold 
right up to this point and usually will, since in general only one coefficient of 
stability will change sign at a time. To sum up, the relation between secular 
stability and ordinary stability may be stated in the following four propositions: 

(i) // an equilibrium configuration is secularly stable it is also ordinarily stable. 

(ii) // it is secularly unstable, it may be ordinarily stable or it may be ordinarily 
unstable. 

(iii) // it is ordinarily stable, it may be secularly stable or unstable. 

(iv) // it is ordinarily unstable, it is necessarily secularly unstable. 
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These are not, of course, four independent propositions; for example, (iv) 
obviously follows from (i). 

It has seemed worthwhile to go into the relation between the two kinds of 
stability in some detail because it is over this very question that former writers 
appear to have fallen into serious error. For instance, Jeans says:* 

'Further, wo see that as tho physical conditions of a system gradually change, secular 
instability necessarily sets in before ordinary stability. Thus for problems of cosmogony 
it is secular instability alone which is of interest. A system never attains to a con- 
figuration in which ordinary instability comes into operation, since secular instability 
must always have previously intervened.' 

The present discussion makes it perfectly clear that these statements are 
incorrect. Jeans's subsequent discussion of the consequences of his solution of 
the problem rest on the same misapprehension as to what is involved. 

The nature of secular instability 

If we have a configuration that is ordinarily stable but secularly unstable, let 
us consider how it will develop. In the complete absence of friction it will, if 
slightly disturbed, simply oscillate about the equilibrium position. But if friction 
is present, we shall have by (24) 

d (T R + V - i w 2 /) = - ve quantity. 
ctt 

Tn equilibrium we have T R 0, and also V %aj 2 I is a maximum (with respect 
to some of the coordinates) which may conveniently be taken as zero. We can 
therefore imagine the system displaced to an adjacent configuration in which 
V lw z l is negative, so that initially 



where c is a small positive constant. In the subsequent motion either T R9 
V \<*PI, or both these quantities will be different from zero, and hence in any 
event the system will move. But when motion occurs the presence of friction 
causes the left-hand side to decrease further, and after a small interval t we 
shall have ^ + F _^ 2/ = _ c _^ whcre ^ >0 

Thus even if T R should vanish subsequently, V '- ,W 2 / would then be still further 
from the equilibrium value, further motion would therefore necessarily develop, 
and the left-hand side would continue to decrease. Since T R cannot be negative, 
this means that V -\<iPl becomes more and more negative, and the system 
therefore departs more and more from the original equilibrium configuration. 

It is important to notice that the rate of departure clearly depends on the 
amount of friction present and vanishes with it. Hence the departure from 
equilibrium cannot be assumed to be necessarily a rapid one, though however 
small the friction the amplitude of the motion must increase indefinitely given 
time enough. 

* Astronomy and Cosmogony, p. 199. 
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The nature of ordinary instability 

Ordinary instability means that the expressions for the initial motion of the 
system involve exponential terms of the form e pt where p is real and depends 
only on the geometry of the system, as it were, and will in general be a quantity 
of standard magnitude. The rate of departure from the equilibrium position will 
therefore in general become rapid, and the system will continue to move until, 
if friction is present, a new configuration of relative equilibrium is reached. 
But the equations of small motion indicating the instability can reveal only the 
initial stages of this motion and will not in general give any indication of the 
new steady state to which the system will eventually move. 

Systems of two degrees of freedom 

The following investigation, due to Lamb, shows how the foregoing theory 
applies for systems of two degrees of freedom, and establishes that a system 
can be ordinarily stable while secularly unstable. 

Neglecting friction to begin with and absorbing w into the j8's and the positive 
quantities a v a 2 into the coordinates, the equations of small motion may be put 
in the form n . , _ 



Assuming in the usual way </ x = g 10 e A ', q 2 = #2o eA ' an( ^ eliminating the ratio 
<7io : feo. we find 



& 1 6i = 0> (34) 

The two values of A 2 will be real if 



or W + 2b : + 26 2 ) + (6 X - 6 2 ) 2 > 0. (35) 

There are three cases to consider: 

(i) b and b 2 both positive. Condition (35) is certainly fulfilled since both terms 
are positive for any value of /?. Also both roots of (34) are negative since 



The system is therefore ordinarily stable whatever the value of /?. 

(ii) 6 X and b 2 of opposite sign. The first form of (35) shows, in the present 
case, that the requirement is necessarily satisfied. But the roots of (34) must be 
of opposite signs since AfA| is now negative. The positive root then leads to 
terms of the form e pt , and hence it follows that the system is ordinarily 
unstable. 

(iii) 6j and b 2 both negative. In this case (35) will obviously hold if )8 2 is large 
enough. If we suppose 2& x + 26 2 -f ]8 2 > 0, then 

A 2 A 2 >0 and A 2 + Af<0. 

Both roots are therefore essentially negative and the system ordinarily stable. 
The modified potential energy function is however now a maximum. 
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These results are in accord with the earlier investigation, and the last case 
demonstrates the possibility of ordinary stability holding in a secularly unstable 
system. 

Effect of friction 

Let us consider how the motion will be modified in case (iii) if frictional forces 
proportional to the velocities are present. It will be shown that in a certain 
sense the system is now unstable. 

The equations of motion may be written 



where /n,/i 2 >/22 are small quantities, the coefficients of the so-called 'dissipation 
function' given by ^(/ii'/i + 2/i2'A < /2+/22</l)- This function is always positive, so 
that work is always being done against the system, reducing its energy. 
Accordingly 



/u>0> ^^ 

The equation for the periods now becomes 

A 4 + (/u+/2 2 )A 3 + (& 1 + & 2 + j32+^^ = 0, (36) 

which differs from (34) by the presence of the various terms in the coefficients 
of the dissipation function. 

The original roots being + iX l9 iX 2 we may suppose the modified roots to be 

Pi i(*i + **i), P* i(X 2 + r 2 ), 

where Pi,p 2 ,r l ,r 2 must all be small, since they vanish with the/^'s. The sum of 
the four roots gives at once 



and the sum of the reciprocals of the roots gives, to the first order in the small 
quantities, k> - . 

^P\ , f/'2 __ _JU _/22 

A? A! \ 6 2 >U> 

since b t and b 2 are both negative. Hence p 1 and p 2 are of opj)osite signs, and 
therefore if the system is disturbed one of the oscillations will die away, and 
the other will increase exponentially at a rate depending on the/^'s. 
Eliminating p lt say, we find 



Hence if X l < A 2 , p 2 is negative and 2h is positive. Accordingly that oscillation 
increases (p 1 positive) corresponding to the smaller frequency (here Aj), that is, 
the one corresponding to the longer of the two original periods. The amplitude 
of the oscillation of shorter period decays. 

Stability when the angular momentum remains constant 

So far we have been considering dynamical systems in which w is maintained 
constant and oscillations can be regarded as taking place relative to these 



26 STABILITY 

rotating axes. If, however, we have a freely rotating system in relative equilibrium, 
so that cu is constant for this state, it is conceivable that if given a small displace- 
ment, the subsequent average rate of rotation would be slightly different from 
co, but that in a frame rotating at this different rate, the motion is nevertheless 
expressible by real harmonic terms. The position is that whereas an ignorable 
coordinate, such as an orbital angle, can increase indefinitely without necessarily 
producing great change in the general form of the path, any palpable coordinate 
must change only slightly under disturbance for a system to be stable. 

For instance, a particle describing a circular orbit under the inverse square 
law will be at relative rest in a certain rotating frame. If it is given a slight 
tangential impulse, its total energy will change and therefore also its mean 
motion round the centre of force. In the original frame it will therefore depart 
further and further from its equilibrium position and an investigation of its 
stability based on F Jco 2 / would suggest instability. 

On the other hand, the condition 



with co assumed constant is always a sufficient condition for stability, but not a 
necessary one. Accordingly, if stability is indicated by means of this function, 
then stability is assured, but if instability is indicated it may not mean a 
physically real instability, and further examination of the problem may be 
necessary. A more general criterion will next be obtained. 

Freely rotating systems 

Let us refer the system to a set of rectangular axes Ox, Oy, Oz rotating about 
Oz with angular velocity co, as yet unspecified and not assumed constant. Then 

we have seen that m 

T = 



If the z-axis is taken as the invariable angular momentum axis, then H is constant 
for all oscillations or motions of the system. Suppose now that we choose co 
always so that the relative angular momentum H R is zero. This defines a certain 
co, which will not necessarily or usually be constant with time, but which reduces 
to the usual rigid-body angular velocity when the system is in a steady state. 
We then have 



(37) 

// = co/, (38) 

and elimination of the variable co gives 

T = T R + H*I>2I. (39) 

For a freely rotating system, if V is the potential of the external forces, the 
energy equation T + V constant, now assumes the form 

T R + V + # 2 /2/ = constant, (40) 

which is similar in form to (25) except that the constant H now appears in place 
of co, and the term H 2 /2I replaces the former term Jco 2 /. 
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Exactly similar arguments to those already given for systems with constant 
rotation now show that systems of constant angular momentum will be secularly 
stable if F + 7/ 2 /27 is an absolute minimum, and therefore also ordinarily stable. 
Otherwise, the system will be secularly unstable and may or may not be ordinarily 
stable. 

The importance of this criterion was first pointed out by Schwarzschild. 

Cases in which the two criteria are equivalent 

It can readily be shown that if the small displacements of the system from 
the equilibrium configuration are such that the moment of inertia / is unchanged 
to the first order of small quantities (that is, changed only in second order) the 
two criteria are equivalent. 

For since H is invariable, we have S(o>7) = 0, and hence 

oa) = - 87, 
and is also of second order of smallness. We then have that 



correct to the second order of small quantities. Hence the condition for stability 
may be obtained by consideration of 8F Jo> 2 87, which is exactly as if the 
function V |o> 2 7 were used with to treated as constant. Thus the conditions of 
secular stability are equivalent in these circumstances. 

Where ordinary stability is concerned, since to is unaltered except for second- 
order changes, the terms in the first-order oscillation equations, corresponding to 
(28), that involve 8a> can be neglected, since they are of second order. Thus 
constant to may again be assumed in considering small oscillations, so that the 
conditions for ordinary stability may also be derived on this basis. 

EXAMPLES OF STABILITY OF ROTATING SYSTEMS 

(i) A particle of unit mass moves under a central attractive force which at distance 
r is of amount \L\r n , where n is regarded as a slowly changing parameter. 

Here V p\(n I)/-"" 1 , and 7 = r 2 . Hence if h denotes the orbital angular 
momentum, stability will depend on the function 

u-v + n _*^ l+ . 

Configurations of relative equilibrium are therefore given by 

dU _ ?_h 2 _Q 
~dr ~ r n r 3 ~ ' 

This gives h 2 = /n/r 71 " 3 , which is simply the condition for a circular orbit of 
radius r. 
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To examine its secular stability we have 

. n ^ , W _ * 



Hence such an orbit will be secularly stable if n < 3. 

The conditions implied by this calculation would correspond physically to 
the presence of a resisting medium rotating round the centre of force in such 
a way that the particle had zero velocity relative to it only when moving in a 
circular orbit. If n< 3, a slightly non-circular orbit would have its radial motion 
gradually damped out, and the system would be secularly stable. 

It may easily be shown that this is also the condition for ordinary stability. 
For if in the radial equation of motion 

r-r6 2 = 



we put r = a + , where is small and a is the radius of the undisturbed orbit, 
and eliminate the angular rate 6 by means of r 2 6 = h + 8h, say, this gives 



and the condition that the coefficient of is positive reduces to n < 3. 

It is to be noticed that because of the presence of the term 8h/a* on the right- 
hand side, the oscillations are not about r = a as mean value but about a slightly 
different value (if S7&4=0), and the mean motion in the disturbed path will there- 
fore differ to the first order from the original value. Thus if measured in a frame 
rotating with the original angular velocity, the motion would appear to be 
unstable, though in fact it would not represent a true instability in the sense 
that a palpable coordinate increasing indefinitely would do. 

(ii) A particle of unit mass moves under gravity on the inner surface of a spherical 
bowl of radius a that is made to rotate round its 
vertical axis with constant angular velocity CD. 

If 6 denotes the inclination to the vertical of 
the radius drawn to the particle, the function 
appropriate to the system will be 




The equilibrium positions are therefore given by 

- = ga sin w 2 a 2 sin 6 cos = 0. 
du 

Hence = 0, or cos = g/aw 2 , a position that 
exists only if aw 2 > g. 

The position = is evidently ordinarily stable always, since if friction is 
ignored the rotation of the bowl has no influence. In the position given by 
cos 6 = g/aa) 2 , the particle revolves as a conical pendulum and no friction is 
required to maintain it in this position. Also such a motion is well known to be 
ordinarily stable. 
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To consider the secular stability of these positions, we have 

- = ga cos 6 - w 2 a 2 (cos 2 6 sin 2 6). 
(iu* 

For 9 = 0, this gives U e9 = ga co 2 a 2 >0, and hence the lowest position is stable 
so long as acu 2 <f/. But when the upper position becomes possible, the lowest 
position is secularly unstable. Also, for cos = g/acu 2 , U 00 = o> 2 a 2 sin 2 0> 0, and 
hence the upper position is always secularly stable. At to 2 = g/a, there occurs a 
transfer of stability to a new series of configurations, and the system therefore 
provides an example of case (i) of the theory of exchange of stabilities (p. 8). 

Initial motion of the particle near the lowest position 

A discussion of the small motion in the neighbourhood of = 0, when this 
position is secularly unstable, affords an interesting example of the results 
reached on pp. 24-25. 

The equations of small motion relative to horizontal rectangular axes Ox, Oy, 
rotating with the bowl may be written 

x 2ajy co 2 x = 2kx n 2 x, 
y + 2a)X a> 2 y = 2ky n 2 y, 

where & is a small positive constant representing the friction, and n 2 = {//a. 
Writing = x + iy, these equations can be combined into 

+ (2iw + 2k) + (n 2 - o> 2 ) f = 0. 
Assuming f oce A ', the period equation is 

A 2 + (2tw + 2k) A + (n 2 - o> 2 ) = 0, 
or 

Hence to the first order in k we have 

A = i< 



Hence x + iy = 

Relative to fixed axes at the lowest point the coordinates are simply the real 
and imaginary parts of (x + iy)^**. Clearly, if a><n, that is if aw 2 <g, both real 
exponential terms decay, and the lowest position is stable, as already shown. 
But if w> n, that is aa> 2 > g, the first exponential term gradually increases, so 
that the distance from the origin must increase. The term of increasing amplitude 
corresponds to the oscillation of longer period, 27r/(n CD), in the motion that 
would obtain if k were zero. The amplitude of the term of period 27r/(n + a>) 
always decays. 



30 



STABILITY 



(iii) A particle moves on a wire in the form of a curve possessing an asymptote at 
x = a, and is attracted towards the z-axis with a force of amount A# 3 . 
Here 



U = V- 

Equilibrium positions are therefore given by 
dU 



dx 



= Xx 3 a> 2 x = 0, 



that is, x == 0, or x = <J(a) 2 /X). Hence there is a 
real position of equiJibrium, other than x = 0, 
only if a) z < Aa 2 . 

To examine the stability of this position, we 
have 




Fig. 8. 

Hence the position at x = is unstable if a) > 0, and the position at x = ^/(o> 2 /A) is 
stable. Thus if at slowly increases there is a linear series of equilibrium con- 
figurations that terminates at x = a, where at 2 = Aa 2 . 

Where ordinary stability is concerned it is readily seen that the position 
x = is ordinarily unstable, and the other position ordinarily stable, and this 
is what would be expected for a system of one degree of freedom. 

The fundamental assumption of applied mathematics 

It may perhaps be worthwhile, in conclusion of this chapter, to mention 
explicitly the general assumption that underlies almost the whole of applied 
mathematics. It is the assumption that the accurate solution of approximate 
(usually linear) differential equations of motion is approximately equal to what- 
ever would be the solution, if this were obtainable, of the accurate (non-linear) 
equations governing the system. There is certainly no rigorous general mathe- 
matical justification of this step, but it has long since represented the standard 
procedure of applied mathematics, and indeed, quite apart from its seeming 
naturalness, to certain minds at least, is forced on us by the insuperable diffi- 
culties usually lying in the way of any other method of solution. There is, 
however, the a posteriori justification arising from the now numerous problems 
in which a subsequent comparison with observation or experiment has been 
made showing considerable measure of agreement. But nevertheless, from the 
point of view of logical procedure, the step is an assumption without as yet any 
rigid mathematical demonstration of its validity. 



Chapter m 
THE SPHERICAL FORM 

If a gravitating mass of liquid of uniform density has no angular momentum it 
is evident intuitively that its equilibrium form will be spherical. 

However, to show that the sphere is a possible form, let p denote the pressure, 
p the density, and V the gravitational potential at any point of the liquid. The 
condition for hydrostatic equilibrium is 



dp = pdV. 
Since p is constant this has the integral 

p = p V + constant. 



(1) 



(2) 




Fig. 9. 



For a uniform sphere V is constant over the surface, and hence the surface 
condition, p = 0, will be satisfied for an appropriate value of the constant. 
Moreover, when (2) holds, (1) must also hold at all points of the liquid, and 
hence the sphere is a possible equilibrium form. 

That it is the only form has been proved by Lichtenstein, who shows by simple 
considerations that any non-rotating equilibrium form must be symmetrical with 
regard to planes through its mass centre. The sub- 
stance of the proof is as follows. Suppose the mass 
divided up into elementary cylindrical columns of 
infinitesimal cross -sections and with their axes all in 
the direction of the normal to some plane through 
O the mass centre. This plane may be taken as 
Oxy and the direction that of Oz. The mid-points 
PuPg, ... of the columns will define a certain surface. 
If this surface is plane, the property certainly holds 

that the configuration is symmetrical about it. Assuming therefore that it is 
not plane, there must be some point M on it at which the value of z is a maximum, 
say. (The argument would apply equally well with slight changes if it were a 
minimum.) If, as in Fig. 9, AB is the column through M , then in relation to 
every other column the point A is at greater gravitational potential than the 
point B. Hence for the whole mass V(A) > V(B), and therefore by (2) p(A) >p(B), 
whereas for an equilibrium configuration the pressures at points on the surface 
must be equal. Hence there cannot be a point such as M at which z is a maximum, 
or minimum. The surface must therefore be plane and the form symmetrical 
about it. For a non-rotating system every plane through O must therefore be 
a plane of symmetry, and the equilibrium form must be a sphere. Also, since 
the density is uniform and fixed, this sphere must be unique. 

It may be mentioned here that the proof remains valid if the mass is rotating 
about any axis parallel to Oz, for the contribution of the centrifugal force to 
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the potential is the same at A and B. Hence for rotating systems the plane 
through the mass -cent re perpendicular to the axis of rotation must always be 
a plane of symmetry of an equilibrium form. 

Stability 

It is also intuitively evident that the spherical form will be thoroughly stable 
for small displacements, but this can be readily established by means of spherical 
harmonic analysis, and the proof affords a simple illustration of the general 
method that is employed to deal with rotating systems. 

Suppose a denotes the spherical radius and that after a small displacement 
at the free surface, without change of total volume, the radius at the point of 
spherical polar coordinates (r, 6, <f>) is given by 

r = a(l+6(0,<)] (3) 

where s is everywhere small and is expressible as a sum of surface harmonics, 

tllUS oo 2,M1 

= S AySMO.ft, (4) 

M l w=l 

there being 2n+l independent functions S, and where A are constants small 
compared with unity. It will be convenient to write this more shortly as 



To the first order of small quantities, the gravitational potential at external 
points due to a uniform mass whose boundary is the surface (3) is, by a well- 
known result, 



where g is the value of surface gravity on the original sphere r = a. On the 
surface (3) the value of the potential apart from a constant can be written to 
the first order 



To find the difference V V Q of the potential energies in the displaced and 
spherical forms, we adopt the usual method of supposing the free surface to be 
r = a(l + ks), where & is a constant, and that a layer of depth adk.s is added 
to this, the final form being gradually built up by increasing k from to 1. This 
gives for the work required, correct to the second order, 



where the double integral is to be evaluated over the surface of the sphere 
r = a. By the orthogonality properties of the surface harmonics this reduces to 

da. (8) 
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The coordinates of the system, infinite in number, can be regarded as the -4 n 's, 
there being 2n + 1 such coordinates involved in a general displacement of order n. 
The coefficients of stability are accordingly 



and these are essentially positive for all n > 1. The potential energy of the sphere 
is therefore an absolute minimum and the configuration therefore secularly stable. 
If n = 1, the three coefficients of stability vanish, as would be expected, for 
the harmonics are then x, y and 2, and correspond to small translations without 
change of shape of the sphere as a whole parallel to the coordinate axes. For 
such displacements the spherical form is evidently neutral. 

Small oscillations 

Tf we assume that oscillations are brought about by external gravitational 
disturbances the resulting motion will be irrotational. Tf <t> is its velocity potential, 
the equation of continuity will be Laplace's equation, 

V 2 <D = 0, 
since p is uniform. The appropriate general solution will be of the form 

*> 2/j 1 r n 

<i> = s s B;~s:(9,<i>), 

terms with negative powers of r being excluded, since those would lead to 
infinities in the velocity. As before, it will be convenient to write this solution 



Since the oscillations are considered small, the coefficients B n will all be small. 
At the free surface, we shall have, to the first order, 



and therefore aA n = B H . (10) 

fi 

We have already seen that the potential W , say, at the free surface is, apart 
from a constant, given by (7), so that substituting in the pressure equation, 
namely ^ 

-= W + constant, 
p t't 

and expressing the condition that p is constant at the surface, we get, equating 
coefficients of the variable functions 8 to zero, 
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Hence for the part A%S% = 8%, say, of the displacement we have 

2n(n-l)g 
~ 



It is evident therefore that for displacements of order n, the 2n + 1 independent 
surface harmonics S% are normal coordinates for the system. Putting now 
s n oc e lV ', if the displacement contains only $, the corresponding frequency 
o- n , say, is the same for all 8% and is given by 



'a "a""' 2n+l ' 

where G is the constant of gravitation. 

The frequency is therefore always real and the system accordingly ordinarily 
stable, as follows in any case from the secular stability already established. 
It will be noticed that the periods depend only on the density and not on the 
linear dimensions of the system. 

The kinetic energy T of the motion at any stage can be calculated from the 
formula 



taken over the surface, where -- denotes normal differentiation. Since O is small. 

ov 

this may be written with sufficient accuracy 



(14) 

These coefficients are, of course, essentially positive, as we should expect. 

The equation of energy for the motion is T + V = constant. But consideration 
of secular stability requires only the calculation of V, for T must automatically 
always be positive definite. If V is also essentially positive, none of the co- 
ordinates can increase beyond a small amount settled by the energy of the 
initial disturbance. 
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THE SPHEROIDAL AND ELLIPSOIDAL FORMS 

To demonstrate the existence of the Maclaurin spheroids and Jacobi ellipsoids 
as possible equilibrium forms we require first an expression for the gravitational 
potential of such figures at internal points. If we consider an ellipsoid referred 
to its principal axes as coordinate axes, its equation may be taken as 



and it will be assumed throughout that a^b^c. By a well-known result, the 
gravitational potential at internal points (x, y, z) of a uniform mass bounded by 
(1) is the quadratic expression in x,y,z 

W= - 



where A 2 = (a 2 4- A)(6 2 + A)(c 2 + A), and p denotes the density and G the constant 
of gravitation. This may be written more shortly 

W 7rGp(otx 2 + fty 2 + yz 2 8), (3) 

, f 00 dA n , f 00 d\ 

where a = a6c /^--vnr, ^ = 6c /T - 2 



A (4) 

and 8 = a6c -r-. 

Jo A 

The potential energy of the whole mass, which will be denoted by F, is given 
by the volume integral 

F=-i \pWdr 

taken throughout the ellipsoid (1). To evaluate this, we have in the first place, 
by the moment of inertia formulae for a uniform ellipsoid, 



\x 2 dr = \Gp.\TTobc, \y 2 dr = lb*.%nabc, \z 2 dr = Ic 2 .$7rabc. 
Hence V = f 7r 2 #/> 2 (i a a 2 + J-j86 2 + iy c 2 - 8) 



To reduce this, we have 



2 A: = 
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and therefore 

a 2 6 2 c 2 / 1 1 1 \ 

c 2 + A / 



Hence V = - iVr 2 #P 2 a 2 fe 2 c 2 ( 1 + - A - ---,- ) -r- 

Jo \ A aA/ A 



But A/ A vanishes for both A = and A = oo, so that we arrive at the well-known 
result 



-. 

(5) 
Jo 



A 

for the total potential energy of the ellipsoid. 

For a spheroid, a = b, if we write c 2 = a 2 (I e 2 ), the foregoing integral can be 
evaluated in elementary functions, to give 



-e 2 ) 1 / 6 -- . .(6) 



From the method of calculation, it is seen that the zero of energy in both 
these results corresponds to infinite dispersion of the matter. 

To express these energies in dimensionless form it is convenient to adopt as 
unit of energy the quantity 



and if the unit of length is chosen so that r = (a6c) 1/3 = 1, then (5) becomes 

V= - iV -A- for ellipsoids, (7) 

Jo A 

sin ^ f 
and (6) becomes 7 = - f ( 1 - e 2 ) 1 / 6 ' for spheroids. (8) 

6 

Condition for a configuration of relative equilibrium 

The equations of relative hydrostatic equilibrium of a liquid mass referred to 
rectangular axes Ox, Oy, Oz rotating round Oz with constant angular velocity CD are 



(9) 



fa ' 
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where p is the pressure, p the density (supposed uniform), and W(x,y,z) is the 
gravitational potential of the mass in the equilibrium configuration. This system 
of equations has the integral 

? = W + %a) 2 (x 2 + y 2 ) + constant, (10) 

and if the configuration is one of relative equilibrium it must be possible to 
choose the constant to make p vanish on the surface of the distribution. When 
this is done, equation (10) determines p everywhere, and the equilibrium condi- 
tions (9) are each automatically satisfied at all points of the liquid. Hence the 
necessary and sufficient condition that a given configuration, with gravitational 
potential W, shall be an equilibrium form with angular velocity o> is that the 



(x* + y*) (11) 

shall be constant over its surface. 

If we denote this expression by <f>(x, y, z), say, and the equation of the surface 
by /(#, y, z) = 0, the foregoing condition is simply that d<f> for all displace- 
ments such that / = and df simultaneously. That is 



for all dx,dy,dz on the surface/ = such that 



Now for a displacement on f(x, y, z) = 0, any two of dx, dy, dz may be regarded 
as arbitrary, and hence the condition is equivalent to 



dx dij By 
when f(x, y, z) = 0. 

General ellipsoidal forms 

To demonstrate the existence of spheroidal and ellipsoidal forms it is now 
only necessary to show that the conditions can be satisfied for surfaces of these 
types. Thus, for an ellipsoidal form, by means of (10), the pressure is given by 



P " 
where G is a constant. The surfaces of equal pressure may accordingly be written 



X y 

and for the surface + r5 + -s = 1 to coincide with one of these, the conditions 
a 2 6 2 c 2 
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(12) become 



and these are equivalent to 



and 

The first of these is a relation between the axes. The second determines (for 
given p) the appropriate value of w in terms of the axes, or for given a> is another 
relation between the axes. 

Substituting for j8 and y from (4) in the last term gives, after a little reduction, 



which form shows that we must always have b > c. On the other hand it is evi- 
dent from the first form for w involved in (14), that we may have either a^b or 
6>a, but it will be convenient throughout to assume always a ^6, there being 
no physical difference between the two series. Thus in any case the rotation 
must always be about the least axis, and hence there are no equilibrium con- 
figurations having the form of prolate spheroids. 
From the definitions (4) of a, j8, y we have 

abc dX 

A)^ + A)-A' 

and hence the first of conditions (14) may be written 

j?-* 

and this is a general relation between the axes of all ellipsoidal forms. It is 
clear that there may be two possible ways of satisfying it. Either, we can take 
a = b t so that the first factor vanishes; or, supposing it possible, as in fact it is, 
with a =|=6, we may make the integral factor vanish. The former solution corre- 
sponds to the Maclaurin spheroids, and the latter to the Jacobi ellipsoids. We 
proceed to consider in turn these two modes of satisfying the relation. 

The Maclaurin spheroids 

If a = 6, the integrals for a, j3, y, and 8 are expressible in terms of elementary 
functions, and the formula (15) for the angular velocity reduces to 



where e is the eccentricity of the meridional sections. It will be noticed that 
w 2 /27rGp is independent of the size of the spheroid and is a function solely of 
the degree of flattening. 



THE MACLAUBIN SPHEROIDS 

To calculate the angular momentum H of the system, we have 

H = 

and M = 

H* 6 
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= (a6c) 1 / 3 , then 



.(1-e 2 )- 2 



(18) 



The kinetic energy T at any stage is Jfa 2 <o 2 , and with GM 2 r~ l as unit of 
energy this may be written 

, 3 CJ 2 CD 2 

rp _ _ . MQ^ 

lO'r 2 '27rGp' { * 

The following table gives a series of values of the various quantities for the 
Maclaurin spheroids. 

TABLE I. MACLAUIUN'S SPHEROIDS 
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0-1890 


0-6249 


0-4898 


0-1664 


0-99 


9210 


0-2710 


0-1551 


0-7121 


-0-4509 


0-1717 


0-9912 


1-9622 


0-2597 


0-144.8 


0-7277 


0-4436 


0-1719 


1-0 


oo 


0-0 


()() 


00 


00 


()() 



This table is based on values given by Lamb and taken originally from Thomson and 
Tait, but the values of I' and T have been computed by the author and also the whole 
of the values appertaining to the three spheroids of eccentricity 0-9299, 9529, and 0-9912. 
The first of these is the form of maximum angular velocity, the second represents the 
form at which ordinary stability ceases, and the last is the form possessing maximum 
kinetic energy. 

Reference to this table makes clear the following results, which can also be 
rigorously established by analytical means from the foregoing theory. 

(i) Every oblate spheroid is a possible equilibrium form if the angular 
momentum of the mass is suitably assigned. 

(ii) The maximum possible value of the angular velocity is given by 
w 2 = 2irGp x 0-2247 and occurs for e = 0-9299, o/c = 2-7198. 
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(iii) For each value of cu 2 less than this there are two possible spheroidal forms, 
but possessing different angular momentum. For values of o> 2 greater than this 
no spheroidal forms exist. 

(iv) The angular momentum increases monotonically and indefinitely as e 
increases from to 1, so that for given angular momentum there is only one 
possible spheroidal form. 

(v) The potential energy, as it evidently must, increases monotonically with 
e from -0-6 for the spherical form to zero for the infinite disc, which is the 
limiting form of the mass for e = 1. 

(vi) The kinetic energy begins at zero for e = and rises to a maximum value 
of 0-1719 at e = 0-9912 and then decreases to zero for the infinite disc. 

The Jacob! ellipsoids 

To demonstrate the existence of ellipsoidal equilibrium forms, with three 
unequal axes, we consider the second way of satisfying condition (16), namely 



which, since A 2 = (a 2 -f A) (6 2 + A) (c 2 + A), may be written 



P 



- = 0. (20) 



If in this equation we write Ba for a, 6b for fe, 6c for c and 2 A for A, where 6 is 
any constant, the equation is unaltered. Accordingly it must represent a relation 
between the ratios a: b: c only. That is, given a: b say, it determines a: c; or given 
a and b it determines c, and, as can be shown, a unique c, and then any multiples 
0a, Ob, 0c of these values also satisfy it. It is therefore always permissible to 
choose the unit of length so that the product of the axes of the ellipsoid abc = 1. 

That a solution of (20) always exists for any ratio of a: b can readily be shown. 
For if on the left-hand side we put c 2 = 0, the integral remaining is essentially 
positive, whereas if we put c 2 = a 2 6 2 /( 2 + 6 2 ) the integral is essentially negative. 
Since this second value of c 2 is necessarily less than both a 2 and 6 2 , it follows 
by continuity that a root must exist less than the smaller of these, and hence 
that there is always some value of c, less than both a and b, for which the equa- 
tion is satisfied. This means that for given a: b, there is not only always a 
corresponding c, but its value leads to a real value of cu, by (15). 

The angular velocity, from (14), may be put in the form 



- - 
2nG P ~ Jo (o + A) (6 + A) ' A ' 

so that, as for the spheroids, the non-dimensional quantity a> 2 /2w6r/0 depends 
only on the ratios of the axes and not on their absolute lengths. 
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The angular momentum H of the system can be expressed in dimensionless 
form as follows: 



abcr 

The potential energy V at any stage is given by (5), and the kinetic energy 
is now TQM(a 2 + b 2 )w 2 , which in the same units as for the spheroids and with 

abc 1 becomes 

/<i" 

(23) 



For the ellipsoidal forms, the integrals involved can no longer be evaluated 
in terms of elementary functions. They are, however, expressible in terms of 
Legendre's elliptic integrals of the various kinds, as Darwin has shown, but 
numerical values are much more awkward to obtain than for spheroids. Their 
computation appears to have been tackled only by Darwin, if we except some 
much earlier values given by Plana, several of which Darwin considered to be 
wholly erroneous. Table IT, calculated originally by Darwin, gives a set of 
values for possible ellipsoidal configurations beginning with the spheroidal 
member of the Jacobi series. 

TABLE II. JACOBI ELLIPSOIDS 



a 


b 


C 


o> 2 /27rG> 


II 


r 


r 


197 


1-197 


OG98 


0-187 


304 


0-585 


0805 


21 6 
270 
383 


1-179 
1-123 
1-015 


0-698 
690 
0-692 


0-187 
0-186 
0-181 


0-304 
0-306 
0313 


(0585) 
(-0584) 
-0-581 


(0-0806) 
(0 0809) 
0-0817 


601 


0-924 


0677 


0-166 


0-341 


-0-561 


0850 


88.58 
-899 


0-8150 
0-811 


6507 
0-649 


1420 
0-141 


3896 
392 


(-0-552) 
-0 551 


(0 0894) 
0901 


2 340 


0-702 


0-607 


0-107 


0-481 


-0-519 


0-096 1 


3 129 


0-588 


543 


0-066 


0-639 


-0467 


0-1006 


5 041 


0-452 


0-439 


026 


1 009 


-0355 


0993 


Lf> 


00 


00 


0-0 


LO 


0-00 


000 



In this table the values enclosed in brackets have been interpolated from 
Darwin's values and not directly computed. The values in the sixth row given 
to an extra place of decimals correspond to the ellipsoid of bifurcation at which 
the pear-shaped series crosses the Jacobi series. 

The table shows that for given a and b, there is only one value of c satisfying 
(20), a result first obtained by Meyer. It is also seen that the angular velocity 
is greatest for the spheroidal member and thereafter diminishes steadily as the 
series is described in the direction of increasing angular momentum. The figure 
of bifurcation common to the Maclaurin and Jacobi series can readily be found 
by putting a = b in (20) and writing c 2 = a 2 (l e 2 ). The integral can be evaluated 
in simple terms and the condition reduces to 



(24) 



which has for solution e = 0-8127. 



42 THE SPHEROIDAL AND ELLIPSOIDAL FORMS 

With increasing angular momentum, b and e tend to equality, but with b 
always greater than c until the limiting form with infinite angular momentum 
which consists of an infinite cylinder with a = oo,6 = c = 0, and abc = 1. 

As for the spheroids, the potential energy increases continually as the ellipsoid 
elongates with increasing angular momentum. The kinetic energy, however, 
rises to a maximum of about 0*1010 for a certain ellipsoid slightly more elongated 
than the member for which a : b : c = 3-129 : 0-588 : 0-543. 

Analytically there are, of course, two Jacobi series branching off the Maclaurin 
series, but they are geometrically and physically identical, and involve only an 
interchange of a and 6. It is therefore only necessary to consider the branch 
for which a ^ b. 

STABILITY OF THE SPHEROIDAL AND ELLIPSOIDAL 
FORMS FOR RESTRICTED TYPES OF DEFORMATION 

Having established the existence of the Maclaurin and Jacobi series, the 
next question that arises is that of their stability. But we are not in a position 
to consider general displacements of such systems for as yet we have no means 
of evaluating the changes in gravitational potential energy and moment of 
inertia for a general deformation. If, however, we restrict the displacements 
simply to those for which the free surface remains a spheroid or an ellipsoid, 
with the least axis always coinciding with the axis of rotation, the potential 
energy and moment of inertia are given by expressions already available. If 
under such restricted conditions a system is shown to be stable, nothing more 
than this will have been achieved, but if on the other hand a system is shown 
to be unstable, then in spite of the restricted freedoms involved, the system 
can be regarded as unstable from a physical standpoint, and there is no strict 
necessity to consider any other kinds of displacement since a general dis- 
placement may be assumed to contain contributions to all possible types of 
deformation. 

The discussion of stability in these restricted conditions that we are about to 
consider also provides interesting illustrations of t the application of the criterion** 
of stability developed earlier and will help to prepare for the more general 
discussions given in the subsequent chapters. 

Stability of spheroidal forms for increasing angular velocity 

Suppose we consider the Maclaurin series and take the angular velocity of 
the coordinate axes as the gradually increasing parameter. Also let us suppose 
that the deformations contemplated are such that the form always remains 
spheroidal with the axis of symmetry coincident with the axis of rotation. 
The system then has only one degree of freedom and the single coordinate 
required for its specification may be taken as e the eccentricity of the sections 
through the axis of rotation. This example is, of course, artificial from a physical 
standpoint, since angular velocity is not the physically significant factor in a 
free system. 
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The stability or otherwise of the system then depends on the total mechanical 
potential 



where V is given by (8), and the moment of inertia by 

/ = f Jfa 2 = f Mr*(l -e 2 )- 1 / 3 . 



We thus have 



( i _ 



+ 2t( 1 - e 2 )- 1 



where k = a) 2 /27rGp and remains constant during the deformations. In the 
present example a> will not be the angular velocity of the spheroid at all times 
during the deformations, but only of the equi- 
librium form in the neighbourhood of which 
the motion is taking place. It is of course the 
angular momentum in a fixed frame that 
remains constant. It will be convenient to 
use e to represent the eccentricity during the 
deformations and also the eccentricity of the 
particular member of the Maclaurin scries 
about which the motion is supposed to occur. 
Omitting the constant factor \\7r*Gp 2 r*, 
equilibrium configurations will be given by the 
vanishing of the derivative, thus 




Fig. 10. Graph of o> plotted against c 
for Maclaurin spheroids. 



---- 

c j > 

and this equation is satisfied, as it should be, when k has the value given by 
(17), so that the Maclaurin spheroids are the equilibrium forms. To examine 
their stability, we have, omitting positive constant factors, 






It may readily be shown from this expression that -- 2 - vanishes for e 0, but 

immediately thereafter is positive and remains so until that value of e given by 

sm- a e = (l-e 2 ) 1/2 (9e-2e 3 )/(9-8e 3 ) 

at which it changes sign and then remains negative. It can easily be verified 
from (17) that this value of e corresponds precisely to that value at which 

-j-- = 0. Accordingly stability, under the present conditions, disappears, as would 

CtC 

be expected, at the maximum value of the angular velocity. It will be shown 
later, however, that this is not a true instability. In fact, the result means only 
that relative to axes rotating at a rate greater than w max , no spheroidal form 
of equilibrium is possible, and the coordinates in any steady state if measured in 
such a frame would appear to depart more and more from any given initial values. 
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Stability of the spheroids for increasing angular momentum 

If the angular momentum H is adopted as parameter and the deformations 
again restricted to those maintaining the spheroidal form, the Maclaurin series 
can readily be shown to be thoroughly stable, and the preceding result accordingly 
shown to be a spurious instability arising purely from the choice of axes. 

In the present case the appropriate function for discussing secular stability is 



where 



// = Ma 2 aj and / = 



= f Mr*(l - e 2 )- 1 / 3 . 



The displacements at any stage are represented by varying e but with // remain- 
ing constant. We thus have 



77 16 

U = - - 77 



16 

15 



, sin" 1 e 5// 2 



, sin" 1 e 

e 



where K is a factor that remains constant during 
the deformations and is given by 



15 



_ _ 

K ~ ' 

2 



q> 



'2(1- e 2 ) 2 / 3 ' 




Omitting the constant factor $%7r 2 Gp 2 r 3 , the equi- Fig. 11. Graph of // plotted against 
librium forms are, as usual, given by e for Mada ^P^">^- 



dU 
~ de 



_ e 2-5/6 i n -l 



e 2) 



| _ }\ + e -l( 1 _ e 2)-I/3 



_ e 2)-2/3 = Q 



and, as before, this equation is satisfied for the value of K given above. To 
examine the question of stability, we have 



From this it may easily be shown that ,- 2 - vanishes only for e = and that 

thereafter it is always positive. It follows that every member of the Maclaurin 
series is secularly stable for displacements that preserve the spheroidal form 
with axis of symmetry remaining the axis of angular momentum. 

The foregoing results concerning the stability of the spheroids could have 
been deduced at once from the considerations of Chapter IT. Thus, when o> is 
taken as parameter, we see from Fig. 10 that at a certain value of e the derivative 

-v- = 0, so that stability, in the coordinate system implied, is lost. On the other 

U/& 
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hand, if H is taken as parameter, it increases monotonically with e, and there 

is no configuration at which -- = 0. Since the spherical form is stable, the 

etc 

series must therefore be stable everywhere. 

Stability of spheroids and ellipsoids for certain ellipsoidal deformations 

If the displacements are assumed to be such as to preserve always an ellip- 
soidal form with the axis of rotation a principal axis, the question of secular 
stability can again be dealt with quite simply. 

If we plot the values of ajc and b/c for the Maclaurin and Jacobi series, deter- 
mined by Tables I and II, as rectangular coordinates, we get the diagram of 
Fig. 12. In this diagram the point ti corresponds to the spherical form (a = b = r), 

8r 



b/c 



disc 




Last stable ellipsoid C ( cy | inde r a=oo, b- c -0) 

1 

Jacobi ellipsoids I 



B/C 



Fig. 12. Diagram showing representative points of Jacobi and Maclaurin figures. 

D corresponds to the infinite thin disc-spheroid (a = 6 = oo, c = 0), and C corre- 
sponds to the infinitely elongated cylindrical ellipsoid (a =00, fr = e = 0), which 
is the limiting form of the Jacobi figure. The point B is the point of bifurcation 
of the two linear series. 

It is important fo notice that we are considering a system defined by three 
unequal axes and so requiring two coordinates for its specification. If we start 
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with the Maclaurin spheroids we are limited to the case a = 6, which requires 
only a single coordinate, and the existence of the Jacobi ellipsoids does not 
manifest itself. This is why no maximum or critical value of the angular 
momentum, corresponding to the point B, occurs in Table I for the Maclaurin 
forms. On the other hand, the value of H at the point B for the Jacobi series 
has a critical value. Thus the spheroids can be regarded as a special case of 
ellipsoidal forms, so that there are two series of ellipsoidal forms intersecting 
at B t namely, the Maclaurin and Jacobi series. But there is only one series of 
spheroidal forms. 

To discuss their stability by diagrammatic means we can imagine the linear 
series of Jacobi forms described by beginning with H = oo and the angular 
momentum gradually diminished, that is, we can use H as a decreasing para- 
meter. Then, at first, for sufficiently large H, there are three possible equilibrium 
forms (two equivalent ellipsoids and one spheroid), but the parameter eventually 
reaches a value where dH = 0, that is, at B. At this stage the forms with three 
unequal axes disappear, and for smaller H only the spheroidal form is possible. 
The system thus provides an instance of case (i) of Chapter II (p. 10). Since 
the spherical form is stable it follows that the Maclaurin series is stable for the 
portion SB of the series and that it loses its stability and is always unstable 
thereafter. Moreover, since the Jacobi series then comes into existence for 
greater values of //, stability will be transferred to it. Plotted in terms of 
angular momentum as vertical coordinate the Jacobi series turns upwards. 

On the other hand, if the angular velocity were used as parameter and plotted 
as vertical coordinate, it is clear from Tables T and IT that the ellipsoidal series 
would turn downwards after B. If the system were regarded as evolving along 
the Maclaurin series SB, with increasing angular velocity, the ellipsoidal forms 
could not come into existence, and in this sense the Jacobi series is unstable. 



Analytical proof of stability for ellipsoidal displacements 

The foregoing discussion of secular stability for restricted deformations is 
based on arguments given by Jeans. A simple analytical proof, under the same 
restrictions as to the permitted deformations, has been given by Lamb. The 
general lines of his discussion appear to be sound, but it is incomplete in that 
certain points necessary to the argument are left unestablished. These points 
are covered in the proof now to be given. 

For secular stability, as the angular momentum is increased, the function 
V + H 2 /2I = /(a, 6), say, must be an absolute minimum for ellipsoidal deforma- 
tions in which the angular momentum remains constant. We have from (5) 



and since abc = r 3 , this expression, as indicated, may be regarded as a function of 
a and 6 only. 
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Whatever value H may have, this function has the following values at the 
points indicated: 

for a = oo, /(oo, 6) = 0, 

for b = oo, /(a, oo) = 0, 
for a = b = 0, /(O, 0) = oo. 

Moreover, /(a, 6) must always exceed the finite negative value of the first term 
corresponding to the spherical form. 

/* d\ 

Also, if a = 0, so that be = oo, -> 0, 

Jo A 

577 2 
and hence /(0,6) = " .6~ 2 , 

57/2 
and similarly, if b = 0, /(<*, 0) = ; ^r.a~ 2 . 

It is also necessary to the argument to show that as a 2 + 6 2 ->oc, f(a,b) tends 
to zero through negative values provided neither a nor 6 vanishes. To prove 
this, let a = rc~ l , b = rc" 1 ^, so that c = re 2 " 8 , where 0<< 1, so that c->0 and 
a,6->oo, as ->0. Then the gravitational term, omitting the constant factor, 
becomes 

r dx 



Jo {(r 2 e- 2 + A) (r 2 6- 2 + 2 * + A) (rM-* + A)} 1 / 2 

r /v e 2 ~MA 

= "A^JO {(r 2 +6 2 A)(r2 + e 2-2A)(r 2 *-2 + A)} 1 / 2 . 

Now suppose that X = ~*, where fc>0, so that K-+OO as e >0, and that 
A r e 2 ->() and A'e 2 ~ 2s ->0. These require jfc<2 and o<2-2-Jk. The terms in the 
denominator involving e can then be ignored, and the integral becomes 



The second term in f(a,b) is clearly 0(e 2 ). Hence the negative gravitational 
term will predominate at infinity provided 



Since this is certainly true for O 0, it follows that the function approaches the 
infinite circle from below except at the points on the axes a = 0, 6 = 0. The 
secular stability can now be readily established as follows. 

If (i)H< 0-304 #*.MM = 7/ c , say, then Tables I and II show that there is 
only one possible equilibrium form and that this is a Maclaurin spheroid. Hence 
a critical value of /(a, 6) must occur for a certain point such that a = b. If 
therefore a, 6, and / are plotted as rectangular coordinates with the a, b axes 
horizontal and / vertical, the resulting surface (for // < H c ) must be of the form 
shown in Fig. 13. Since /(0,0) = +00 and/(oo,oo) = -0, the critical point M in 
the plane a = b corresponds to a negative value of / and must be an absolute 
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minimum with respect to all variations on the surface, as otherwise there would 
have to be further critical points on the surface. The Maclaurin spheroids for 
H<H C are accordingly secularly stable for displacements of the kind permitted. 




a-b 



(Absolute) Minimum 



Case (i) H<:0304 



Minimum J 1 (Jacobian ellipsoid) 
Case(ii) H^0304 



Fig. 13. Case (i) //)-301. (W (11) //>() 304. The point M in the plane a - - 6, repre- 
senting u Maclaurin form is a minimum for displacements in this plane but a maximum 
for displacements perpendicular to it. 

4(fl=t=&) an absolute minimum, and a second symmetrical point ./ 8 , give together the 
two stable Jacobi forms. 



If (ii) #> 0-304 #Uf M, Tables I and II show that there are three critical 
values of /, one of which always corresponds to a spheroid and is therefore in 
the plane a = b y and the remaining two correspond to equivalent Jacobi ellip- 
soids, with a and b interchanged, and therefore symmetrically placed with regard 
to the plane a = b. Let us denote these critical points by J and J 2 . The plane 
through J v say, and the vertical axis meets the surface in a certain curve, the 
ordinate of which is +00 for a = b = and -0 at the infinite circle a 2 + 6 2 = oo, 
/ = 0. The critical value at J l is therefore negative and finite, and it follows 
that r/j and */ 2 must be points where /is an absolute minimum. Also the point M 
must now be a so-called minimax, since it is a minimum for displacements in 
the plane a = b but a maximum for displacements transverse to this. This 
arrangement of critical points is the only possibility consistent with there being 
three, and only three, stationary points on the surface. 

It can therefore be concluded that the Maclaurin spheroids are secularly 
stable for ellipsoidal displacements of the present kind for e< 0-8 127, though 
this by no means enables us to assert their stability for general deformations. 
On the other hand, their secular instability for e> 0-8 127, which the foregoing 
argument establishes, means that they may be considered physically unstable 
beyond this limit, and that they would not therefore come into existence as a 
result of an evolutionary process that could be represented by increasing angular 
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momentum. It will be shown later that the spheroids are secularly stable for 
all displacements if e < 0-8127. 

Where ordinary stability is concerned, it follows that if the oscillations are 
restricted to forms of the present ellipsoidal type, the spheroids are ordinarily 
stable for e< 0-8127, but for values greater than this nothing can be inferred 
from the present work as to the ordinary stability. It has, however, been shown 
by Cart an that they remain ordinarily stable for all displacements provided 
e < 0-9529. But physically the spheroids are of little interest beyond the point 
where their secular stability ceases. 

For the Jacobi forms the argument has shown that they are always stable 
for these restricted ellipsoidal deformations, but nothing can be inferred from 
this result concerning their stability, or otherwise, for general displacements. 
It will prove in the sequel that the Jacobi forms remain stable for all displace- 
ments up to the configuration for which 

a : b : c : (a&c) 1 / 3 = 1-8858 : 0-8150 : 0-6507 : 1 
and // = 0- 



But the deformation through which instability first enters is only disclosed 
when the displacement is analysed in general form by means of ellipsoidal 
harmonics. 

Evolution along a linear series 

In the foregoing discussions the angular momentum has usually been adopted 
as the gradually changing parameter, while the density of the mass is assumed 
to remain constant. It can be shown, however, that the series of configurations 
obtained is exactly the same, as to its proportions, if the density gradually 
increases while the angular momentum remains fixed. To sec this, since the 
total mass and angular momentum are given by 

M = \irpabc, II = \M(a* + &V, 
there exists the following relation involving only ratios of lengths 



Suppose now that the total mass M , the angular velocity o, and the ratios of 
the axes are fixed, then the value of the right-hand side is fixed, and the relation 
can therefore be satisfied for any assigned p by appropriately choosing //, or 
for any assigned H by appropriately choosing p. It follows that the linear series 
given in Tables I and II obtained by increasing H could equally well be obtained 
by increasing p with H remaining constant. Since, when p is maintained constant, 
we can assume abc = 1 always, this usually makes it more convenient to regard 
the series as described by increasing // with the density remaining constant. 

When the angular momentum is kept constant, the initial spherical member 
of the Maclaurin series will correspond to zero angular velocity, zero density, 
and infinite radius, while the final members of both the Maclaurin and Jacobi 
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series will have infinite density. It is easily shown from Tables I and II that 
when H is kept constant the angular velocity increases always with the density. 
To determine the values of cu and of p for any given form, that is, for given 
ratios of the axes as already tabulated, it is necessary, since p is no longer constant 
along the series, to choose some unit of density. For this it is convenient to 
select that value corresponding to the figure common to the Maclaurin and 
Jacobi series, and to choose the unit of angular momentum so that the constant 
value of H is 0-3035. It may further be supposed that the radius which the bifurca- 
tion figure would possess if changed to spherical form without change of density 
is the unit of length. 

Since the dimensions of H are the same as of G*M*r*, it follows that for any 
other configuration, of given ratios of axes, the value of its spherical radius r 
must be inversely proportional to the square of the tabulated angular momentum, 
and hence that the density must be proportional to the sixth power of this 
quantity. Thus, to take particular examples, for the figure of bifurcation we have 

a : b : c= 1-197 : 1-197 : 0-698, 
w 2 /27rG = 0-1871, H = 0-3035, 
and hence for the critical Jacobi figure for which 

a : b : c = 1-886 : 0-815 : 0-651 and H (tabulated) = 0-3896 
we have a> 2 /27rG = 0-142 x (0-3896/0-3035) 6 = 0-636, (H = 0-3035). 
Again, for the last but one tabulated member of the Jacobi series for which 

a : b : c = 3-129 : 0-588 : 0-543 and H (tabulated) = 0-639 
we have co 2 /27r<? = 0-066 x (0-639/0-3035) 6 = 5-75, (H = 0-3035). 

In the case of uniform density, o> decreases to zero as the angular momentum 
increases to infinity and the kinetic energy of the motion has a maximum value 
after which it decreases to zero for the infinitely elongated form. Since the 
kinetic energy must also be zero in the infinitely dense final form, when the 
density is used as parameter, it must again have a maximum value for some 
density intermediate between and oo. These results appear to have been first 
established by G. H. Darwin. 



Chapter V 
ELLIPSOIDAL HARMONIC ANALYSIS 

In the present chapter we develop the elements of ellipsoidal harmonic analysis 
as far as is necessary in establishing those properties that are required for the 
subsequent chapters. Standard results in the theory of spherical harmonics 
will be assumed. 

Confocal coordinates 

Consider the system of quadrics represented in rectangular cartesian co- 
ordinates x, y, z, by 2 2 



where 6 is a variable parameter, and for which it will be supposed that a>b>c. 
They have the simple property that through any point pass three different 
members of the system. If x,y,z are fixed, the equation can be interpreted as a 
cubic in whose roots are the parameters of these three surfaces. This equation 
may be written 

f(6) = (a 2 + 6) (b* + 6) (c 2 + 6) - x*(b 2 + 6) (c 2 + 0) 

- V 2 (c 2 + 0) (a 2 + 0) - z 2 (a 2 + 0) (6 2 + 6) = 0. (2) 

To locate the roots, we have that /( a 2 ) is negative, /( b 2 ) is positive, /( c 2 ) 
is negative, and /( 4- oo) is positive. Hence the roots of the cubic are all real 
and different. If we denote them by A, /*, v, taken always so that A > ft > v, then 

A> -c 2 >/Lt> -b 2 >v> -a 2 . (3) 



It follows that the A-surfaces are ellipsoids, the /^-surfaces are hyperboloids of 
one sheet, and the v-surfaces hyperboloids of two sheets. It is also easily shown 
that these three families of surfaces are everywhere mutually orthogonal. 
It is evident from symmetry that three given surfaces 

A = constant, p = constant, v constant, 

will intersect in eight points, one in each octant defined by the coordinate planes, 
whose coordinates can be written ( #, y, z). In discussing the properties of 
the surfaces it is therefore sufficient to consider simply the positive octant. 

Returning to the original equation (1), if A, /*, v are the parameters of the 
three surfaces through the point (x, y, z), then we must have as an identity 



Multiplying both sides by a 2 + and then putting = a 2 (and so on in turn) 
we get the following expressions for the rectangular coordinates (x, y, z) in terms 
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of the orthogonal curvilinear coordinates (A, /K, v): 



2= 



(a 2 -6 2 )(a 2 -c 2 ) ' 



Z (a 2 -c 2 )(6 2 -c 2 ) ' 

By direct addition of these, or more simply from consideration of the sum of 
the roots of /(0) = 0, we have 



r* = x* + y* + z* = A + /z + y + a 2 + 6 2 + c 2 . (6) 

Since for given finite values of a, 6, c the quantities /x and v always lie within 
finite limits, by (3), it is evident that at great distance from the origin the para- 
meter A must tend to infinity. That is 



Lam functions and Lame polynomials 

Let us suppose next that f(t) denotes a rational integral polynomial of degree 
m in the variable t, and let us write 



Then the product LMN will be a polynomial symmetrical with regard to 
A, p,, and v. 

If the factors of/(J) are t-8 v t-8 2 , ..., 8 /n say, assuming the coefficient of 
t m in f(t) to be unity, we shall have 



and hence that LMN = Ft {( A - S r ) (p - 8 r ) (v - 8 r )} 



by means of (4), apart from a factor independent of x,y,z. Hence LMN defined 
in this way will always be equal to a product of quadratic factors Qi,Q 2 , >9 m 
each of the form 



If, as will be shown later is possible, the function (7) so resulting, which we 
may denote by V(x, y, 2), is such that it satisfies Laplace's equation 



then the original function L(X) is called a Lami function of the first kind, or of 
the first species. 
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If similarly we commence by putting 

(i) L=/(A)V(A + a*), M=f(n)j(n + a?), N - 
or (ii) L=/(A)V(A + 6), M =/( M ) V(/* + & 2 ), N = 

or (iii) i= 



where in each case /(A) is again a rational integral polynomial whose term of 
highest degree is A, then in exactly the same way it is seen from (5) that the 
product LMN will be a symmetrical polynomial in A, /z, v multiplied in case (i) 
by x, in case (ii) by y, and in case (iii) by z. Moreover, the rational part of the 
product LMN will again be expressible as a product of quadratic factors 
Q r (x,y,z), so that the whole expression LMN when put in terms of x,y,z will 
be of one of the forms 

v = to' 



Tf V resulting in this way satisfies Laplace's equation, the corresponding function 
L(A), which now contains one of the radicals ^/(A + a 2 ), ^/(A + 6 2 ), or ^/(A + c 2 ) as a 
factor, is called a Lame function of the second kind. 
Similarly, if we put 

(i) 

L = /( A) V[( A + & 2 ) (A + c 2 )], M = f(p,) J[(p + 6 2 ) (/* + c 2 )], JV = /(v) V[(v + 6 2 ) (v + c 2 )] ; 

or (ii) 

L = /(A) Vl"( A + c 2 ) ( A + a 2 )], Jf = /(/,) V[(^ + c 2 ) ( M + a 2 )], JV = /(v) Vl(v + c 2 ) (v + a 2 )] ; 

or (iii) 

(v-f 6 2 )]; 



where again /(A) is a rational integral polynomial, we arrive by means of the 
product LMN at a function of one of the forms 



and if V 2 F = 0, then the original function L(\) is called a Lame function of the 
third kind. 
Finally, if 



L = /( A) VK A + 2 ) ( A + 6 2 ) ( A + c 2 )], M = /(/,) ,/[(/. + a 2 ) (/* + 6 2 ) (/x + c 2 )], 

^=/WV[(v4-a 2 )(v 
the resulting function V will be of the form 



54 ELLIPSOIDAL HARMONIC ANALYSIS 

and if V 2 V = 0, then the original function L(X), which now contains 



as a factor, is termed a Lamb function of the fourth kind. 

The several cases may be combined schematically as follows: 




= V(x,y,z). (9) 

When V 2 V = 0, the polynomial V(x, y, z) is called a Lame polynomial or an 
ellipsoidal harmonic, and to each species of Lame functions there corresponds a 
species of ellipsoidal harmonics. 

The possible forms that the Lame functions L(A) can take must depend ulti- 
mately on the requirement that V must satisfy Laplace's equation, and we 
proceed to consider how suitable functions L(A) may be constructed. 

Laplace's equation in confocal coordinates 

Since the surfaces A, p,, v are mutually orthogonal, the square of the line 
element ds corresponding to a displacement dX, dp,, dv will be of the form 



and it is easily found that 



'1 ~~ 4//I.2 



4(a 2 + A)(6 2 + A)(c 2 + A)' 



hi 



3 



If we define positive quantities R, A, B, and C as follows: 



(10) 



then 



7? 



7? 



(11) 



(12) 



(13) 



Expressed in general orthogonal curvilinear coordinates Laplace's equation is 
well known to become 



3 , ia - o 

dX\h 1 ^ + ~ ' 
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and in the present case, since A, B, C depend only on A, /LC, and v respectively, 
this may be written on multiplying through by ABC 



o"-* 

Normal solutions of Laplace's equation 

Suppose that we seek solutions of equation (14) of the form 



wherein, as indicated by the notation, L is a function of A only, M a function of 
fi only, and N a function of v only. Solutions so expressible are sometimes termed 
normal solutions of the equation. For such solutions we have 



and accordingly Laplace's equation becomes 

x C d 



which is of the form 



This equation, if V is a solution of V 2 V = 0, must be true for all A, /z, v. Putting 
p, = v gives at once ^ 2 (/ A ) = ^(/^ an( ^ thus all three functions <f>, < 2 , and < 3 must 
be identical, and such that 



= 0. 
To find <, if A is put zero, this equation can be written 



and since this must hold for all JJL and v, each side must in fact be constant. 
The form of the function < is thus 



where H and K are constants, and the differential equation for L is accordingly 



For reasons that will be apparent in due course it is convenient to replace // by 
a different constant, n, related to it as follows 

J(n +!) = //, 

and it will be seen later that the important cases are when n is an integer. The 
equation for L then becomes 

0. (16) 
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This linear differential equation for L is known as Lame's equation. 

The functions M and N satisfy exactly similar equations, but with /LI and v 
as the independent variable in the respective cases. 

As yet there are no restrictions on the values of H and K, so that whatever 
their values this equation possesses a general solution which may be written 
), and it determines a solution of V 2 F = 0, namely, the product 

F = 



Since Laplace's equation is linear, it follows that any linear sum of solutions of 
this kind is also a solution. 

Polynomial solutions of Laplace's equation 

It will be shown that it is possible to choose the constants II and K so that 
Lamp's equation has solutions of the following forms: 

(i) a polynomial in A; 
or (ii) a polynomial in A multiplied by owe of the radical factors 



factor, 2K 2K 2V W m _. 1 



or (iii) a polynomial in A multiplied by any two of these radical factors; 
or (iv) a polynomial in A multiplied by Vt(A + a 2 ) (A-h& 2 )(A + c 2 )]. 

Such a solution, if it exists, may be written in all cases, apart from a constant 

^ Am _. 1 + j 

say, 

where the indices K lt K 2 and /c 3 arc each independently equal to either or J, 
m is a positive integer, and /(A) a polynomial of degree m. 

If M and N denote respectively the corresponding expressions with ju, and v 
in place of A, the product L(X)M(fji)N(v) will be expressible, by means of (4) and 
(5), as a polynomial in #, y, z in which the terms of highest degree are of order 
2(* 1 + /c 2 + /c 3 + w). This polynomial will not of course be homogeneous. If its 
highest terms are of order n (its identity for integer cases with the n of (16) 
will appear shortly) we shall have always 



so that if n is even all three or one of the /c's must be zero, and if n is odd two 
or none of the /c's must be zero. 

If now this form for L is substituted in equation (15), it must become an 
identity, and this as a first condition requires H to be of a special form. To 
show this we have 



'-* (A + c 2 )*'+* (A + a l \ m ~ l + ...) 
+ K 3 ( A + 2 )* 1+i (A + & 2 )*+* (A + c 2 )"'-* ( A + a x A'"- 1 + . . . ) 
+ m( A + a 2 )**+* (A + & 2 )*^ * ( A + c 2 )*>+* (X m ~ l + ...). 
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Considering at this stage only the highest power of A, its order is 

K i + K z + *3 + i + m = 
and its coefficient is i 
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Hence when L is substituted in the left-hand side of (15), the term of highest 

order occurring will be 
fo 



This must therefore equal the term of highest order on the right arising from 
(HX + K)L, and accordingly we must have as a first condition for a solution 
possessing one of these special forms, and leading eventually to a harmonic 
polynomial of order n y 



where n is a positive integer. 

Before proceeding to show that A" also can be appropriately chosen for solu- 
tions of the present kind to exist, we notice that the polynomial solutions, in 
x, y, z, of V 2 V = can be classified under four different types. Table III shows 
the four cases and the relation between the degree of V the Lame polynomial 
and L the Lame function in the different cases. In this table j^ denotes a poly- 
nomial in x 2 , y 2 , z 2 of the form 

r 2 

F = nl-A-. 



TABLE III 



Species 


Form of solution in 
rectangular coordinates 
(Lame polynomial of 
degree ?) 


Form of corresponding 
Lame function L 


Degree of /(A), 
the rational 
part of L 


No. of values of 
K = no. of 
different L's 


(0 


r=F(,r a ,y*,2 2 ) 
n even 


/(A) 


i" 


KM + 2) 


(ii) 


r = \y *v 2 , ?/ 2 >z 2 ) 

u 

11 odd 


(V(A 1 2 )j 
U'<A+&") /(A) 
U(A}c*)j 


l(ii-l ) 


i(+i) 


(iii) 


M 

V = lzx\F(jr*,y*,z*) 

w 

n even 


x [(A-ffo 2 )(A + ^)j| 
x t(A+^)(A + ^)j /(A) 
N L(A-} 2 )(A+^)jJ 


J(-2) 


in 


(iv) 


V = xyzF(x\y\z*) 
n odd 


Vl(A|a 2 )(A+^ 2 )(Afc 2 )J/(A) 


i(-3) 


Kn-1) 



The last column in the table gives the number of different values of K that can 
be assigned, for given n, in order that the solution shall be of the form shown 
in the third column. We proceed now to demonstrate that K can be suitably 
chosen. There are four cases corresponding to the four species of Table III. 
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Determination of K 

(I) For functions of species (i) the value of n must be even, and the function L(X) 
will then be a polynomial of degree m = \n. It therefore contains \n + 1 co- 
efficients of which that of the highest power has already been taken as unity. 

Thus we have r M MI MO 

L = A^ + OjA* 11 - 1 * a 2 A* n ~ 2 +..., 



where the coefficients a 1? a 2 , ... are to be disposed in such a way that this is a 
solution of Lame's equation. For this, the substitution of this expression in 
the differential equation for L 

A9 d 2 L A dA dL /rr . . _ 
A* - + 4-. = (HX + K)L 
dX 2 dX dX 

must produce an identity. As already seen, the highest power of A occurring 
on the left is \n+ 1, and the resulting expression must contain %n + 2 coefficients 
corresponding to the terms in A* w+1 , A* 11 , ..., A, A. Because of the choice of 
U = \n(n+ 1) the terms of highest degree are the same on both sides, and there 
will remain \n + 1 coefficients that must be equal on the two sides. The expression 
of this yields \n+ 1 equations, and these are evidently linear in K and in the 
coefficients, |n + 1 in number, of the polynomial assumed for L. Elimination of 
these coefficients gives an equation for K that is essentially of degree n+l, 
and accordingly there are \n+ 1 possible values of A' to each of which corresponds 
a polynomial solution. This is the number given in the fourth column of the 
table on p. 57. 

It will be shown later that all the roots of this equation for K are real and 
different and lead in all to %n + 1 different polynomial solutions. 

(II) Consider next functions of species (iii), so that again n is even. 

In this case L will contain two radical factors (which can be selected in three 
ways) so that the degree of the rational part of L(X) is now \n 1. It accordingly 
contains %n coefficients, and by the same procedure as in the previous case their 
elimination leads to \n values of K. Thus with each of the three possible pairs 
of radical factors there are associated \n values of K. 

Again, it will be shown later that these values of K are all real and different 
and lead to different L functions. 

If n is even only functions of the species (i) and (iii) exist and the total number 
of Lame functions of given even order n is accordingly 

= 2n+l. 



(III) Functions of species (ii) exist when n is odd. The one radical factor can 
now be chosen in three ways, and the rational part of L will be of degree \(n- 1). 
It will contain \(n+ 1) coefficients, and their elimination leads to \(n+ 1) values 
of K, which will be shown to be all real and different. 

(IV) Functions of species (iv) also exist when n is odd. L(X) now contains all 
three radical factors, and its rational part is of degree %(n-3). There are 
accordingly \(n-\) different values of A", which can be shown to be all real 
and different, leading to \(n 1) different Lame* functions. 
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If n is odd the Lame* functions are of the species (ii) and (iv) only and their 

total number is , . - x , , . . x _ 

3.(n+l)-f (n-l) = 2n+l. 

Thus, whether n is odd or even, the total number of Lame* functions is 2?i 4- 1 . 
This result is exactly what would be expected, for each of the Lame* functions 
L(X) leads to a Lame* polynomial V(x, y, z) satisfying V 2 V = 0, and the terms of 
highest degree in V give rise to an associated homogeneous polynomial, V h say, 
of order n also satisfying V 2 V h = 0. Now it may easily be proved by counting 
coefficients that there are, quite generally, 2n + 1 independent homogeneous poly- 
nomial solutions of degree n. The terms of highest degree in each of the Lame 
polynomials will evidently be a certain linear sum of such harmonics specially 
adapted to the ellipsoidal coordinates. 

The terms of lower orders in the Lame polynomial V(x,y,z) are determined 
uniquely by the terms of highest order, for we have 




and the homogeneous part given by the terms of highest degree is 




so that if ^ is given the expression for V can be obtained at once. 

It will be shown later that the 2n + 1 homogeneous polynomials arising from 
the Lame polynomials of degree n are themselves linearly independent, and it 
will then follow that any homogeneous harmonic polynomial of degree n is 
uniquely expressible as a linear sum of them. 

The Lam functions and polynomials of order 0, 1, 2, and 3 

By straightforward application of the foregoing theory it is readily found 
that these are as follows: 

The 1 Lame function of order is 1 (or any constant), corresponding to the 
solution V constant of Laplace's equation. 

The 3 Lame functions of order 1 are 



and these correspond to the harmonic functions x, y, and z. 
The 5 Lame* functions of order 2 are 



together with the two functions 

A + (a 2 + 6 2 + c 2 ) J{a 4 + 6 4 + c 4 - 6 2 c 2 - c 2 a 2 - 
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The first three correspond respectively to the second-order harmonics yz, zx, and 
xy, and the last two have for their homogeneous parts functions of the form 
j(# 2 y 2 ) + k(x 2 z 2 ), where j and k are certain constants, of which there are two 
independent types. 

The 7 Lame* functions of order 3 are J[(X + a 2 ) (A + 6 2 ) (A + c 2 )] and the six func- 
tions obtained by cyclic interchange of a, 6, and c in 



These correspond to xyz together with six polynomials containing terms of 
degree 3 and lower whose homogeneous parts constitute six independent linear 
sums of x 3 - Zxy 2 , x 3 - 3#z 2 , ?/ 3 - 3?/z 2 , ?/ 3 - 3x 2 y, z 3 - 3z 2 z, and z 3 - 3?/ 2 z. 

The relation of Lam polynomials to spherical harmonics 

Before proceeding to prove that the A"s of Lame's equation are all real and 
different it is necessary to consider the analogy between the Lame polynomials, 
which are adapted for use in reference to ellipsoids, and the ordinary spherical 
harmonics. The relationship is of special importance in the case of spheroids, 
when a = b or b = r, for then the Lame functions are expressible in terms of 
certain well-known functions, as will be seen later. The results to be obtained, 
in some instances however, involve a limiting process, and cannot be arrived 
at simply by setting a = b in the general Lame functions. 

Suppose variables 0, <f> are introduced replacing //,, v and defined by 

x = J( A + a 2 ) sin cos <, 

y = > /(A + 6 2 )sin^sm^ > I (17) 

z =V( A + c2 ) cos 0- I 

Equations (5) show that 6, <f> depend on /z, v only. It is important to notice that 
these variables are not the 6, <f> of ordinary spherical polar coordinates. They do, 
however, tend to them as A->oo, for then the radical factors in x,y, z in (17) all 
tend to equality with r. It has already been seen that /A an d v are always finite, 
and it is clear also that for A = oo they depend only on the ratios x : y : z. 

A Lame polynomial corresponding to a certain normal solution L(A) M (/u,) N(v) 
will at great distance from the origin be effectively equal to A* w Jf (^) N(v). But 
since A->r 2 this will also be everywhere exactly the value of the homogeneous 
part of the Lame polynomial. But this is a spherical harmonic, and since the 
coordinates 6, < tend to the 6, </> of spherical polars we must have accurately 

MN = F(sin 6 cos <, sin 6 sin <, cos 6) 

where F is a spherical surface harmonic function of 6, </>. Since this is an identity 
it must hold everywhere, but it is only at infinity that 6, <f> become the ordinary 
polar angles; elsewhere they are defined by (17) and have no simple geometrical 
meanings. Thus when expressed in terms of 0,< so defined, the product MN 
always takes the form of a spherical surface harmonic. But it is a known property 
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that any continuous function of 6, </> can be expressed as a linear sum of inde- 
pendent surface harmonics. It follows that any function of \L, v can be expressed 
as a linear sum of the form 



The sum will in general require harmonics of all orders, so that there will be 
2n+l terms involving the 2n+ I harmonics of order n\ that is, 1 constant term, 
3 terms of first order, and so on. 

Oblate spheroids a-b 

To consider the form the functions take when the confocal quadrics are sur- 
faces of revolution round the z-axis, let it be supposed that b->a. Since 

A> -c 2 >^> -b 2 >v> -a 2 , 



it is clear that v can no longer serve as a parameter when a = 6, for it is always 
equal to a 2 . A third coordinate must therefore be introduced to replace v. 
For this purpose we have from (5) 

^ = _(& 2 + A)(& 2 + /*)(& 2 + v) (a 2 -b 2 )(a 2 -c 2 ) 
x 2 ~~ (a 2 +~A)~(a 2 4- fji) (a 2 + v) ' (a 2 - b 2 ) (b 2 - c 2 ) 



v, A 

_ when a b 

a 2 + v 

since the factors (a 2 b 2 ) cancel accurately always, while pairs of factors such as 
6 2 + A and a 2 + A become equal and non-zero and then also cancel. It is therefore 
necessary to retain only b 2 + v and a 2 + r. Hence we have 



where </> is in this case everywhere the ordinary azimuth angle of spherical polar 
coordinates. If the parameter v is replaced by <f>, so that planes through the 
z-axis are now one of the three orthogonal families, we have 

v a 2 sin 2 <f> b 2 cos 2 <, 
and hence 

and we notice that <f> depends only on v, and the function of <j> replacing N in the 
normal solution is to be obtained by a limiting process involving the foregoing 
relations between (/> and v. 

Where 9 is concerned, /* can take values in the range c 2 to 6 2 , and no 
limiting process is required. We can thus put a b and make x, y, z large to find 
the relation between them. So that 



cos* -J^j^), sin^ ^U J (-?$> 



. sin 2 6 cos 2 6 

becomes 2 + -=--- = 0. 
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Hence sin0= &/* and C os0 



Thus when a = b the parameter p depends on only. It is to be remembered 
that 6 is defined by (17) and reduces to the spherical polar angle only at great 
distance. 

Since, however, the product MN is quite generally a spherical surface har- 
monic in 6 and < as here defined, for spheroids it must become such that M 
is a function of 6 only, and N a function of (f> only. But when ordinary spherical 
surface harmonics are put in normal form, the functions depending on B only 
are the Legendre functions P^(cos 6), and those depending on </> only are the 

circular functions c ? s p<. Hence we must have 
sm 



NM = CO %<A (^ = 0,l,...,orrc). (18) 

v ' sin^ r J 

Since no limiting process is involved in the relation between JJL and 6, the form 
of M(p,) arrived at here must be simply the result of writing a = b in the general 
form of M . Also, the function L(X) in the case a = b can be got directly from 
M (/x) simply by writing A for ^. On the other hand it cannot be obtained from 
the present form for N since this has resulted from a limiting process. 

We have, in fact, apart from an arbitrary constant factor which may be 

omitted 

M = pp t (t) where t = cos 6 



where /> = --, and p can take any of the values 0, 1, 2, ...,n. Tf we write 



then M = (l+T 2 )^7>+ M (l+T 2 )' 1 (19) 

and T, so defined, is imaginary for the function M since c 2 //, > 0. 

To find the function L(X) we now simply replace p by A in the foregoing, 

so that 

L = (1 +T 2 )**>nP+ n (l+T 2 ) n , (20) 

c 2 -f A 
where now r 2 = -^ g, (21) 

and T is now a real variable since A + c 2 > 0. 

To each value of p = 1, 2, ... or n, there correspond two harmonic functions, 
LMN, because of the alternative N factor cosp< or sin^. But for p = 0, there 
is only one N factor, namely 1. We thus arrive at 2n+ 1 harmonics in all, as 
would be expected. 
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Prolate spheroids b = c 

If &->c, it is clear, since c 2 >ft> ft 2 , that in this case p, can no longer serve 
as parameter. To deal with this feature we define 8' and by means of 



so that for b = c the coordinate $ is the azimuth angle in the yz-plane, but only 
at great distance does 6' become the polar angle (measured from the #-axis). To 
find the relation between p, and ^r, we have (in a similar way to the case a = b) 



so that /x = b 2 sin 2 if/ c 2 cos 2 if/. 



Hence CQB*= U and sin * = Lt 

and the function M will depend only on i//. 

To find the relation between 6' and v we put b = c and allow a*, y, c to become 
large, whereupon the equation to the surfaces 

2 4 j2 i /v2 

-* + ? TL' = 1 

tt 2 + V C 2 -f l> 

cos 2 0' sin 2 0' ,. 

becomes at great distance - - - + -^ = 
fo a 2 + v c 2 -fv 

and no limiting process is involved in these relations. The function JV now 
depends only on 6'. 

In the same way as before, since MN when expressed in terms of 6', /r must 
have the form of a spherical surface harmonic, we must now have 

cos 



In the present case N is simply the result of putting b = c in the general Lame 
function, whereas M is obtained by a limiting process. We have, in fact, 

N = (I -j8)p#iH-ii(i -t*), (23) 



where t 2 = cos 2 6' = 



Hence to find L(A), we write t 2 = -^ + A 2 , (24) 

ft "" C 

and then since 1 t 2 is now negative we always obtain L in real form by writing 

t*-l) n . (25) 
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The Lam functions that reduce when a = b to cosp(f> and sin/70. 

To examine the question which of the general Lame functions reduce to 
cosp<f> and which to sin^<, it has been seen that M(fji)N(v) when expressed in 
terms of 0, <j> takes the form of a tesseral harmonic, thus 

M(n)N(v) = P(cos 0) V (V = 0, 1 ..... n). 

bill 

If on the right-hand side of this relation we wish to replace 6 and <f> by p, and v, 
then, apart from constant factors irrelevant to the present argument, the appro- 
priate relations are 

(x) sin cos </> oc <][(a 2 + /z) (a 2 + v)], 

(y) sin 6 sin <f> oc J[(b 2 + p) (b 2 + v)J. 
(z) 



Also J?i;(cos 6) oc (sin 0)| -v-J P n (0 where t = cos 0, 

and P n (t) is the ordinary Legendre polynomial of order n. We thus have 

j ( J P n (OJ {sin ^(cos $ + 1 sin (^)}' J 
= (polynomial in cos 0) x 



omitting constant multiplying factors that make no difference to the particular 
points we wish to establish. 

(a) Suppose now that p is even 

The solution involving cosp<f) is the real part of the left-hand side in this 
equation, and hence requires on the right-hand side the real part of (x + iy) p . 
But we have for this (where the coefficients c r denote constants), 



= rational polynomial in x 2 and y 2 . 

Hence a solution of the form V(x 2 , y 2 , z 2 ), for the case a = b will correspond to 
an even value of p and the term cos p<f> for the N function. (See line 1 of Table TV.) 
On the other hand, sinp< evidently requires the imaginary part of (x + iy) p , 
and we have 

3f (x + iy) = px v ~ l y- . . . pxy*>- 1 

= xy (rational polynomial in x 2 and y 2 ). 

Hence a solution originally of the form xyF(x 2 , y 2 , z 2 ) will correspond to p even 
and the function sin^ for N. (This is shown in line 7 of Table TV.) 

(b) Suppose that p is odd 

A solution involving cosp</> again requires the real part of the right-hand side, 
but now 



= x (rational polynomial in x 2 and y 2 ). 
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Hence a solution of the form xF(x 2 , y 2 , z 2 ) will correspond to an odd value of p 
and cosp(f> for N. (Line 2 of Table IV.) 
Similarly 

!&(x + iy) p = y (rational polynomial in x 2 and y 2 ) 

and hence a solution of the form yF{x 2 ,y 2 ,z 2 ) requires p odd and sin#< for the 
N function. 

Finally, since z depends only on 6, the presence or absence of z as a factor of 
the polynomial solution will not affect the foregoing results. That is, a solution 
of the form zF has the same (odd or even) value of p and function for N as does 
a solution of the form F. We can thus set out in a table the form of the iV-function 
and the appropriate p (odd or even) for the eight possible forms of ellipsoidal 
harmonics for the case a = b. 

TABLE IV 



Form of general 
polynomial 



rf 
yF 
zF 

yzF 

zxF 

JcyzF 



Form of general L : 
/(A) rational 



/(A) 
V( A + ')/( A) 




sin p<f> 



Harmonics of orders 1 and 2 for the case a = b 

When a = b, we have (apart from constant factors) 



x oc 



y oc VI (a 2 + A) (a 2 4- / 



= V[( 2 + A) (a 2 + ft)] sin <^ , 



since v == a 2 and no limiting process is involved. 

As illustrations of the foregoing theory we next evaluate the harmonics of 
order 1 and 2 when a = b. 

(i) The three harmonics of order I 
Here n = 1, arid hence we have 

L = (i+ r *)toJDpV(i +r 2 ) (p = 0, 1), 
where r 2 = (c 2 + A)/(a 2 - c 2 ). 



= T = ^(c 2 4- A) omitting constant factors. 

The function LMN then reduces to <]\(c 2 + A) (c 2 + /*)], which corresponds to the 
ellipsoidal harmonic z. 
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Ifp=ll L 



The function LMN can now take either of the forms 



By the preceding table, line 2 shows that the first of these corresponds to the 
harmonic x y and line 3 shows that the second corresponds to y. 

(ii) The five harmonics of order 2 

Here n = 2 and we have for the general form of L 



0: Z, = Z) 2 (l+r 2 ) 2 

= 1 4- 3r 2 = a 2 + 2c 2 + 3 A apart from constant factors. 

Hence the solution LMN is (3A + a 2 + 2c 2 )(3/z + a 2 + 2c 2 ). 

The corresponding harmonic in x,y, z is easily found, since A and p are the 

roots of 9 9 9 

x 2 + y 2 z* 

a 2 + A c 2 + A ' 
so that A/* = a 2 c 2 - c 2 (x 2 + y 2 ) - a 2 z 2 , 



In rectangular coordinates the solution is therefore 

3(a 2 + 2c 2 ) (a: 2 + if + z 2 ) - 9c 2 (x 2 + y 2 ) - 9a 2 z 2 + 9a 2 c 2 - 3(a 2 + 2c 2 ) (a 2 4- c 2 ) + (a 2 + 2c 2 ) 2 . 

The constant portion of this can be omitted and the rest is easily seen to be of 
the form j(x 2 y 2 ) + k(x 2 z 2 ), where j and k are constants. 



Ifp=l: L= (l+T 2 )*7) 3 (l+ r 2 ) 2 = rV(l+r 2 ) = V[( 

cos 

The LMN solutions are therefore ^[(a 2 + A) (c 2 + A) (a 2 + /it) (c 2 4- /*)] . <f>. Line 6 of 

sin 

Table IV shows that the solution with factor cos^ gives the second-order har- 
monic xz, and line 5 shows that the one with sin <f> as factor gives yz. 

Ifp = 2: L= (l + r 2 )D 4 (l + T 2 ) 2 = l4-r 2 = a 2 + A. 

The LMN solutions are therefore (a 2 + A) (a 2 -h /z) . 2</>, which correspond respec- 

sin 

tively to x 2 y 2 (by line 1 of Table IV) and xy (by line 7). 

These results may be compared with those already given on page 59. The 
solutions when expressed in terms of x, y, z must obviously reduce to ordinary 
spherical harmonics, but the importance of the present analysis is that it enables 
us to find the particular linear sums of these appropriate to an ellipsoid in the 
general case, and to certain spheroids in the cases a = 6, and b = c. 
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The values of K are all real and different 

Having established the form of the Lame* functions in the case of spheroids, 
we are now in a position to return to the question of the values that K must 
have in the general case (a > b > c) in order to give a solution of Lam6 ) s equation 
in the form of a polynomial multiplied by 0, 1, 2 or 3 of the radical factors. 

Supposing first that a = b, we have seen that v must be replaced by <f> related 

to it by 

v = a 2 sin 2 </> 6 2 cos 2 <. 

The equation for the function N, in its general form, is 



Now C 2 = (a 2 + v) (b 2 + v) (c 2 + v) = - (a 2 - b 2 ) 2 sin 2 $ cos 2 <f> (c 2 + v) , 

so that C = i(a 2 b 2 ) sin < cos < ^(c 2 + v). 

Also dv 2(a 2 b 2 ) sin </> cos <f>d<t> 

and hence C = - \i J(c 2 + v) -jj 

dv d<p 

and the equation for N when <f> is the independent variable is 

Idv dN 



If now we let 6 -> a, then in the limit v = a 2 and - , = 0, and the equation for 
N reduces to 



From the discussion of the preceding sections it is necessary that this shall be 
satisfied in turn by 

N = . p<f> (p = 0,1,2, ...,n). 

Accordingly we must have 

n(n + I ) a 2 + 4K = p 2 (a 2 - c 2 ) 
or K = p 2 (a 2 - c 2 ) -n(n+l) a 2 . 

It is thus seen that for a = 6, the equation giving the values of K must have 
these n+ 1 roots corresponding to p = 0, 1,2, ...,n, and they are obviously all 
real and different, and the roots corresponding to odd values of p separate those 
corresponding to even values. 

It may now readily be proved that in the general case when a =4= 6 the preceding 
result still holds in the form that for Lame* functions of any one type of given 
order n the values of K are all real and different and separate those corresponding 
to any other type of the same order. 
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To show this, let us consider first the case when n is even. When a = 6 it 
has been seen there are n + 1 values of K. From Table IV it is seen that for p 
even (0, 2, ...) the general solutions (a 4= b) are of the form 

/(A) and VL(A + a 2 )(A + 6 2 )]/(A), 
while for p odd (1,3,...) the general solutions are of the form 

Vl(A + 6 2 )(A + c 2 )]/(A) and J[(X + c*)(X + a*)]f(X), 

where as usual /(A) denotes a polynomial in A. Also, in the case a = 6, the roots 
in K for even values of p separate those for odd values. Supposing now that 6 
is allowed to vary slowly away from a so that a =j= 6, the foregoing property could 
cease to hold only if at some stage a root K corresponding to a function of the 
first group /(A) Qr V[(A + tt2)(A + 62)]/(A) 

became equal to a root K corresponding to a function of the second group 

or 



It will now be shown that this cannot possibly occur. For if it did, let the two 
corresponding L-functions having the same values of K be L and L 2 . Then if 

u is defined by 7 , 

J d\ = * 



Lame's equation for them may be written 



Hence in the customary way, dashes denoting differentiation with respect to u t 



which integrates to T , r 7 , T , . ~ 

& //i/yg Z/g/^ constant = (7, say, 

or ( 

which in terms of A becomes 

-2V[(A + a 2 )(A + &*)(A + c*)] ( ^(L J /L 8 ) = C\L\. (26) 

It will be shown that unless (7 = this equation cannot be satisfied. 

First it is seen that whatever the form of L 2 its square, Z/|, is rational. Now 
a possible form for L^L 2 , with an obvious notation, is 



and differentiation of this with respect to A cannot introduce a radical factor 
^/(A + a 2 ) anywhere. But the left-hand side of (26) contains ^/(A + a 2 ), and hence 
the equation cannot hold unless C = 0, in which case it reduces to 



= 0, or 
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An exactly similar argument holds in each of the remaining possibilities for 
L 1 /L 2 . But the two L-functions are essentially of different species and hence 
cannot bear a constant ratio to each other. Hence the initial assumption is 
false and their K values cannot be equal; therefore the K values corresponding 
to one species must always separate those of any other species but of the same 
order. It follows also that only real values of K are possible, since, as b changes, 
two real values would have to become equal before complex or imaginary roots 
could come into existence. 

The whole of the preceding argument applies in exactly the same way for n 
odd but with the species of functions now changed in accordance with the forms 
indicated by Table IV. 

We may therefore conclude that for functions of a given species and order, 
the general equation for K when a>b>c has all its roots real and different. 

The linear independence of the 2/i+l Lame functions of a given order 

The Lame functions of order n satisfy 

A dAdL , 



and denoting by A the whole of the operation affecting L on the left-hand side, 
this equation can be written shortly* 

AL = KL. 

Supposing now that a linear relation holds between the 2w + 1 possible Lame* 
functions, of the form 

S c,4 = (8 = l,2,...,2w+l), 
i 

wherein the c s are constants not all zero, it is plain that if we break up the rela- 
tion into parts in each of which the radical factors are the same, then the several 
linear sums so obtained must each vanish identically. For example, if n were 
even, there would be \n+ \ rational Lame functions (of the first species) and 
three sets each of In Lame functions containing a pair of radical factors. Evi- 
dently in an identically zero linear sum of these 2n + 1 functions, the rational 
part, the part containing Vl(^ + ^ 2 ) (^ + r - 2 )l the P ar ^ containing VlX^ + f2 ) (^ + a2 )J 
and the part containing N /[(A + 2 )(A + 6 2 )| must all vanish separately. The same 
must hold when n is odd, and the parts containing J(h + a 2 ), ^(^ + ^ 2 )> V(^ + c ' 2 )> 
and ^[(A + a 2 ) (A + 6 2 ) (A + c 2 )] must vanish separately. 

If we select any one of these parts it will be a linear sum 

S<VL.sO (27) 

s=l 

in which the number of L functions, j, is equal to the number of different possible 
values of K, as given in the last column of Table III. If K^K^...,Kj denote 

* There will be no confusion with the earlier use of A in Chapter IV. 
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these different values of K we shall have for each L a 

Hence, operating on (27) with A we get 

j 



s-i 
A second operation with A gives 

i 



K.C.L. = 0. 



and so on. Operating successively j - 1 times, we obtain in all j linear relations, 
and if the j quantities c 8 L 8 are eliminated we obtain 

1 

= 0. 



But this equation is exactly equivalent to 

i 
i 

and such a result cannot be true, since all the K H are different. Accordingly 
there cannot subsist a linear relation between the Lame functions of a given 
order, and their linear independence is thus established. 

The zeros of the Lame functions 

In the present section we will denote by L r the rational part of any selected 
Lame* function, that is, when radical factors, if any, are omitted. Radical factors, 
if present in a Lame* function, imply zeros at A =?= a 2 , A = 6 2 , or A = c 2 , so 
that if L is regarded as a polynomial in A -fa 2 , as is found convenient in the 
sequel, the presence of radical factors will mean zeros of L at A + a 2 = 0, or 
a 2 6 2 , or a 2 c 2 . We notice also that 

0<a 2 -6 2 <a 2 -c 2 . 

We proceed now to establish the following theorem. 

All the roots of L r regarded as an equation in A + a 2 lie between and a 2 c 2 . 

Suppose first that a h. Then it has been seen that the function L reduces 
to the form L = (!+,*) D*( 



where r 2 = (c 2 + A)/(a 2 -c 2 ), and p has one of the values 0, 1, 2, ...,n. 

The equation L = considered as an equation in r 2 has \p or \(p+ 1) coincident 
roots at r 2 = - 1, according as p is even or odd, together with the roots arising 
from the polynomial factor D p+n (l+T 2 ) n . To establish the distribution of the 
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roots of this polynomial we can begin with the equation (l-fT 2 ) n = 0, and 
consider the chain of polynomials obtained by successive differentiations. This 
equation itself, if r 2 is treated as the variable, has n roots at r 2 = 1. Again, 

the derived equation 

D(l -f r 2 )" = r(l + r 2 )"- 1 = 

has n I roots at r 2 = 1 and one root at r 0. Similarly we get a second 
derived equation 

Z) 2 (l+T 2 ) n = {l + (2rc-l) r 2 }(l+T 2 ) n - 2 = 0. 



Since n>l (for the second differentiation to arise), we have that 2n 1^3, 
and this equation, in r 2 , has n 2 roots at r 2 = 1 and one root intermediate 
between 1 and 0. 

Continuing the process it is seen that each differentiation produces a further 
root between 1 and 0, and all these intermediate roots are separate and inter- 
laced with those of the preceding equation of the chain, since they arise one at 
a time from successive critical values. 

The factor roc^A + c 2 ), when present, gives a root at a 2 c 2 , but since radical 
factors are being supposed omitted from L such a root will not occur. Since 
1-fr 2 = (a 2 + A)/(a 2 c 2 ), roots r 2 = 1 correspond to roots AH- a 2 = 0. Hence, 
when a = b, all the roots of the rational equation 

L r (A + a 2 ) = 
must satisfy < roots in A + a 2 < a 2 - c 2 , 

and moreover the intermediate roots are necessarily all different. 

Before proceeding to show that the same holds good for the roots of the 
general equation L r = 0, we require to establish the following two simple conse- 
quences of Lame's equation for all a, 6, and c. 

(i) Apart from radical factors, L(X + a 2 ) = cannot have a root at 0, a 2 6 2 , 
or a 2 c 2 . 

This means that L r cannot have A -fa 2 , A-f & 2 , or A + c 2 as a factor apart from 
the possible occurrence of these as radicals. To establish this, we have that the 
general function L satisfies Lamp's equation, which may bo written 



Suppose now that L (A-f a 2 )/(A), where/ is a polynomial that does not vanish 
when A + a 2 = 0. Then we have 



)/'=/ for 

,/2 r 

* = 2 /' + ( A + a2 )/" = 2 /' for A + a 2 = 0. 
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Inserting these in the equation and putting A + a 2 = gives 

)/=0 for 



and hence / must contain A + a 2 as a factor. But this is contrary to hypothesis, 
and hence L itself cannot contain A + a 2 as factor. In exactly the same way it 
follows that neither A + b 2 nor A + c 2 can be a factor of L. 

Similarly, if we suppose L = (A + a 2 ) 3/2 /(A), or suppose L = (X + a 2 ) J f(X) where 
j^ 1, and differentiate Lame's equation an appropriate number of times before 
substituting for the derivatives of L and then putting A-fa 2 = 0, we arrive at a 
similar contradiction. 

(ii) The equation L = cannot have a double, or multiple root. 
For, supposing first a double root, then at this value of A we should have 
L = and .- = 0, and hence by Lame's equation 



But, by (i), such a root certainly cannot occur at a 2 , -6 2 , or -c 2 , and hence 
at the root - . 2 = 0. This would mean that it must be a triple root, or higher*, 

and is contrary to hypothesis. Hence no root can be double. 

If we assume a triple, or higher order root, then by differentiating Lame's 
equation an appropriate number of times first we can again arrive at a contra- 
diction. 

We can now return to the equation L r (X + a 2 ) = for the general case a^b 
and prove that all the roots are real and different and lie between and 
a 2 -c 2 . 

For suppose that the equation has one root, or more, actually greater than 
a 2 c 2 . If then we allow 6 to change continuously so that it approaches eqiiality 
with a, this root would at some stage have to pass through the value A + a 2 = a 2 c 2 
since it has been shown that for a b all the roots are less than a 2 - c 2 . But by 
(i) no such root can exist, and accordingly a root greater than a 2 c 2 cannot 
disappear in this way. Similarly there can be no root less than 0. 

There is alternatively the possibility, if the equation L r had two roots 
greater than a 2 c 2 , that these might disappear and become unreal at some 
stage as b approaches a. But, by continuity, two such roots would first have to 
become coincident before becoming complex or imaginary, and this again is 
not possible, by (ii). 

Since when a = b all the roots of L r = are real and different it follows that 
the same must be true when a 4= b. Thus we may finally conclude that all the 
roots of 



are real and different, and are intermediate between and a 2 c 2 . 
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ALTERNATIVE METHOD OF CONSTRUCTION 
OP ELLIPSOIDAL HARMONICS 

The arrangement of the zeros of the several Lame functions of a given kind 
can be specified even more closely than has so far been demonstrated, and it is 
necessary in connexion with the question of the stability of the Jacobi ellipsoids 
for this more detailed location of the roots to be effected. For this purpose we 
shall establish an important theorem due to Stieltjes, but the proof of this 
requires as a preliminary the consideration of the direct construction of ellipsoidal 
harmonics in rectangular coordinates. 

Construction of harmonics of the first species 

Harmonics of the first kind can be put in the form 



For example, to find the second-order harmonics of this type, 

(r 



_ / J_ __1 1 \ 

\ a * + X + l ) 2 + \ + c 2 + A/' 
and hence Q will be a spherical harmonic provided A is a root of the equation 

(A + 6 2 )(A + c 2 ) + (A + c 2 )(A + a 2 ) + (A + a 2 )(A + 6 2 ) = . 
This is a quadratic equation with roots A 1? A 2 , say, such that 

-a 2 <A x < -6 2 <A 2 < -c 2 
or 0< 



There are thus two ellipsoidal harmonics of this kind corresponding to these two 
different possible values of A. 

Next let us consider the conditions for 



to be harmonic. By direct differentiation we have 
VII dll [ 2x 



BII 2_ d 2 II Sx* 

W ~ iW m *~+ 
Hence 



pH* +*,)/' 
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But by partial fractions 

y 

a, 6, c 

and also 



Q P -Q q 



an 



summation over 5 not being implied. Using these relations to remove the second 
derivatives of II we obtain 



2 2 



V' 



Q 



where 2' means, as usual, that the term for which p = q is omitted from the 
sum. It follows therefore that II will be harmonic if A lf A 2 , ..., A m are chosen so 
that each of the separate m quantities in this sum is itself zero. That is, if 
A!, A 2 , ..., A, rt are given by a solution of the set of m equations 



- S' >- > =0, 

9=1 ^l~~ "q 



A,-; 



= 0, 
= 0, 



(28) 



'^- ,- = o. 



These equations, not being linear, will in general have a number of different 
solutions, and to each solution will correspond an ellipsoidal harmonic of the 
assumed form. 

Let us consider now the polynomial L(A) of degree m in A defined by 



9=1 

which we shall show to be the Lame polynomial associated with the foregoing 
ellipsoidal harmonic. Denoting by dashes differentiation with respect to A, we 

have L'( A) = 2 products of A - A 1? A - A 2 , . . . taken m 1 together, 

and L"( A) = 22 products of A - X l9 A - A 2 , . . . taken m 2 together. 

Hence, putting A = A ; ,, it is seen that L"(X p )/L'(X p ) is equal to twice the sum of 
the reciprocals of A^- X v X p A 2 , ..., A p A m , the term A^-A^ not occurring, of 
course. That is r///\ x , 



Hence if A has any of the special values A lf A 2 , ... A p given by equations (28) that 
make II harmonic, the expression 

1 . 1 . 1 , o^'(A) 
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vanishes. Multiplying through by the factors involved in the denominators, 
this means, since L itself is a polynomial, that the expression 

(a 2 + A)(6 2 + A)(c 2 + A)Z/(A) + if 2 (*>* + A)(c 2 + A))z/(A) 

la, b, c ) 

is a polynomial in A that vanishes whenever A has any of the special values 
AI, A 2 , . . . , A m . Accordingly it must possess A X lt A A 2 , . . . , A A m as factors, and 
therefore has the polynomial L itself as a factor. Moreover, it is of degree m+l 
in A, and the coefficient of the highest term is clearly 



ra(ra 
Hence the remaining factor must be of the form 



or, if m = \n %n(n+l)X + K, 

so that L is simply a polynomial solution of the differential equation 



which is precisely Lame*'s equation. 

Harmonics of other species 

An exactly similar procedure of direct differentiation can be applied for func- 
tions of the three remaining species with the same final result, but with slight 
changes in the equations corresponding to (27). 

The Lame function of general species can be written 

m 

(A + a 2 )" (A + b*Y> (A + c 2 )- II ( A - A,), (29) 

P~I 

where the indices /c 1? /c 2 , /c 3 are always either or J. Also, it has already been 
shown that A lt A 2 , . . . , A w are all real and unequal, and different from a 2 , ft 2 , c 2 . 
The degree, n, of the corresponding harmonic is given by 



and when K ly /c 2 , /c 3 are specified, the number of Lame functions of the particular 
form is settled and always equal to ra-h 1. (Table III.) 

If we follow out the same procedure of differentiation to find the necessary 
conditions for 

rx yz 
1 u zx x\ 



xy 

to be a solution of V 2 V = 0, the equations corresponding to (28) in all cases are 
found to take the form 

*3 + ; 
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wherein K V /c 2 , /c 3 have the values or as occurring in (29), and the appro- 
priate sets of Aj/s are given by solutions of these equations. 

We are now able to establish Stieltjes' theorem on the positions and distribu- 
tions of the zeros of the Lam6 functions of a given order and species. 

Stieltjes 9 theorem 

For given /q, K 2 , * 3 , the ra+1 Lame functions L(X) of a given order can be 
arranged in a sequence in such a way that the rational part of the r-th function of the 
set fias rl of its zeros between -a 2 and -b 2 and the remaining m-r+1 zeros 
between b 2 and c 2 . 

To prove this, ]et ^!,^ 2 > "-^m ^ e rca ^ variables, m in number, such that 



-& for p=1 f 2,...,r-l,] 

(31) 
and -& 2 ^<^,^-c 2 for p = r,r+ 1, ...,mj 



and let us consider the following product 



The function P then has these properties. Tt is zero when any one or more of 
the variables <f> p have their least values, and also when they have their greatest 
values. It is essentially positive when the variables ^ p are all unequal and also 
unequal to any of a' 2 , 6 2 , or c 2 . Furthermore P is a continuous and bounded 
function in the domain specified by (31), since each of its vseparate factors is 
bounded and continuous over the ranges concerned. It follows that there is at 
least one set of values of the variables, within the domain, at which P attains 
its upper bound; that is, there exists a set of values within the domain for which 
P, and therefore logP, has a maximum value. These values of (/> V ^ 2J ---^ 
be given by the equations 



dlogP _ dlogP _ _ _ 

<l 302 "" ty m 

But these equations, when written out, are precisely equations (30), with <f> p 
occurring instead of X p . Hence each solution of (30) necessarily lies within the 
domain (31), and not on its boundary. 

Thus if r has any of the values l,2,...,w+l, there exists a corresponding 
Lame function whose rational part is such that r 1 of its zeros lie between 
a 2 and b 2 and the remaining m r+1 lie between b 2 and c 2 . Moreover, 
since when K V K 2 , /c 3 are specified, there are precisely m + 1 Lame functions of 
the given species, it follows that they are all obtained in turn when r is given 
successively the values l,2,...,m+l. 

We notice finally that if r = 1 there are no roots of L(\) between -a 2 and 
-6 2 , while if r = m+ 1 there are no roots between -ft 2 and -c 2 . 
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FURTHER PROPERTIES OF LAM& FUNCTIONS 
AND APPLICATION TO GRAVITATION 

The ellipsoid, E, ?! + j! + *] = i ( 1 } 

can be regarded as the surface A = 0, and any point of its surface (in the positive 
quadrant) can be specified by coordinates (ju,,v). The perpendicular distance 
from the origin to the tangent plane at the point (ju,, v) is easily found to be 

abc 



It will be convenient to denote this perpendicular distance by -arabc so that 
we have 



and is necessarily always finite and positive. This quantity is of special importance 
in establishing the following so-called orthogonal property of the surface har- 
monics associated with an ellipsoid. 

// MN and A^N^ are two different surface harmonics, then 

dS = 0, (2) 



where dS is an element of the surface (1), and the double integral is extended over 
the whole surface of the ellipsoid (E). 

This result can readily be established by means of Green's theorem. For 
consider the two solid harmonics related to these surface harmonics, namely 

V = LMN and ^ = L^M^. 
Then V 2 V = and V 2 T = everywhere. But by Green's theorem 



where the volume integral is taken through the interior of the ellipsoid and the 
surface integral over its surface, and denotes differentiation along the outward 

normal. Since V and V : are harmonic, the left-hand side is identically zero and 
the surface integral on the right must accordingly vanish. 

Now for an infinitesimal step along the normal only the coordinate A changes, 
and hence aT7 ^ T JT %\ 

dL A 
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But by equations (13), Chapter V, p. 54, 



= nr dX on the surface A = 
2abc 



= o T'- 
2-nraoc 

The surface integral thus becomes 



But the factor -- L x depends only on A, and therefore remains constant 
ft A dA 

over (E), or over any A-surface, and may thus be taken outside the integral. 
Also, this factor cannot be zero, for if it were we would have 



^ _^ 

LI dX ~ L dX ' 

or LI = kL 9 

where k is a constant, which would mean that L and L l were not essentially 
different Lame functions. The orthogonality property (2) is accordingly proved. 
It is important to remember that the result holds not only when M N and 
M^NI are of different orders (n), but also for any two different harmonics of the 
same order. 

Development of a function as a series of ellipsoidal surface harmonics 

It has been seen that any continuous function of position (p,, v) on the surface 
of an ellipsoid can be expressed as a linear sum of the surface harmonics M N, 
and the foregoing orthogonality property now enables the coefficients in the 
expansion to be formally obtained. 

Thus for any function O(/A, v) we must have an expansion of the form 



where the A k are certain constants. Multiplying through by TxrM k N k and integrat- 
ing over the surface of the ellipsoid gives 



IT vr$>(tJL, 



and all the products of different harmonics on the right vanish when integrated 
except that for which i = k. Whence 



f f <bM k N k dS = A k f f -urMlNldS, (3) 
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and this relation formally determines the constants A k of the ellipsoidal (surface) 
harmonic expansion of O. It will be convenient to term them the expansion 
constants of the function O(ft,v). It is to be noticed that the factor multiplying 
A k on the right in (3) is an essentially positive quantity. 

The linear independence of the 2n+l ellipsoidal harmonics of a given 
order n 

This can now be easily established as a simple consequence of the ortho- 
gonality. For suppose L S M 8 N 8 (s = 1, 2, ..., 2a+ 1) are the ellipsoidal harmonics of 
a given order n, and let it be assumed that 2n+ 1 constants C A exist, not all zero, 
such that 2w 4 1 

2 c.L.M.N.^Q. 

1 

If we multiply this by mM k N k where M k N k is any selected one of the 2n+l 
ellipsoidal surface harmonics of order n, and integrate over the surface of the 
ellipsoid, we get // 

c k L k \ \-arM 2 N 2 dS = 0, 

all other terms disappearing because of the orthogonal property. But the integral 
is essentially positive, so that we must have c k for all k, and this is contrary 
to the original supposition. Accordingly no such linear relation can exist, and 
the ellipsoidal harmonics of any given order are therefore linearly independent. 

The function 5(X) associated with L 

Since Lam6's equation is a second order differential equation there must exist 
two independent solutions for each value of K. For given appropriate K the 
function L can be regarded as one of these, and there must then exist another 
independent solution which will be denoted by S. 

To determine this second solution, suppose the independent variable A is 
replaced in Lame's equation by u where 

d A = - 2 ^[(a 2 + A) (6 2 + A ) (c 2 + A)] du = - 2Adu. 



-Liien ^j- TX i ^x 7X 

aA 

and Lamp's equation may be written 



du 2 
The function S will satisfy the same equation, so that 



Multiplying the first of these equations by 8 and the second by L, and sub- 
tracting, gives d * s d*L 

JU ,'~o O ~^ ii~ == 0. 

du 2 du 2 

r dS ndL 
or, after integration, L-J S-^- = constant. 



80 FURTHER PROPERTIES OF LAM^ FUNCTIONS 

If this constant of integration were taken as zero, a further integration would 
give S = kL, as seen above. This must, of course, be a possible solution, but it 
is obviously not a new function. If, however, the constant is taken to be other 
than zero, tf then turns out to be a new function altogether. The actual value 
taken for the constant, so long as it is not zero, is immaterial, since changing it 
can only multiply S by a constant factor. It turns out to be convenient to choose 
it to be 2n + 1 where n is the order of the harmonic LMN. We thus have 



- 
L 2 \ du du 



^ ^- 4. 4. 4. - & f" 2?i+l 

which integrates to give T = - T<r - du, 

^ Jo L* 



or, in tenns of A, S W = 

the upper limit of integration again being chosen with a view to later con- 
venience. 

The function S(X) found in this way is therefore such that 



also satisfies Laplace's equation in curvilinear coordinates A, ^, v. 

It will be seen that #->0 as A-^oo. For A large, we have L = aA i/l = ar n , 
where a is any constant factor, and also A -> A 3/2 . Hence, for A large 



f 

1 



a*\ n '2 A 3 / 2 



Thus, in the analogy with spherical harmonics, the function 8MN corresponds 

to functions of the type r~ n ~ l P*(cos 6) . p</>. 

sin 

It thus appears that the general normal solution of Laplace's equation in 
terms of A, //,, v is of the form 



corresponding to the function 

jMcos 0) + ft^(cos 6)} (, cos^ + ft 



of ordinary spherical harmonics. 
It may be noted also that as A->oo 



independently of n and a. 
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The Dirichlet problem for an ellipsoid 

Suppose we are given the values that a function takes on the surface of the 
ellipsoid (1) and we wish to find a harmonic function V defined in all space that 
takes these values on the ellipsoid. 

Let O(/x, v) denote the given function, to which V must reduce on the ellipsoid. 
It must be expressible as a linear sum of surface harmonics in the form 

o = ZA^N!. 

We can replace each constant A l in this by a slightly different constant by putting 



where L t is the Lam6 function in A associated with M t N i9 and A% is the correspond- 
ing function. We thus have, on the ellipsoid 



where ^(0) and $ t -(0) mean the values of these functions on A = 0. 
Now consider the solid harmonic function defined as follows: 

V = Z(*ti(0)LMN inside the ellipsoid,] 
V = X<xL(0)SMN outside. J 



This function evidently satisfies V 2 V = everywhere, and takes the required 
values O on the surface of the ellipsoid, and accordingly represents the solution 
of the problem. 

The close analogy with the comparable result in spherical harmonic analysis 
is obvious. 

The gravitational potential of a surface layer on the ellipsoid 

The above function V is evidently continuous across the surface of the ellipsoid, 

8V 

but its normal gradient is not. To see this, we calculate -_-- on the two sides 

en 

of the surface. At points just outside we have, by outward differentiation, 

f?.g.f* 
du dX dn 



on the surface A = 0, 

where dashes denote differentiation with respect to u. Similarly at points just 
inside we have, by inward differentiation, 
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If we regard V as the gravitational potential of a surface layer of material of 
given surface density <r(jLt,v), then V must be continuous across the surface, but 
its gradient need not in general be so. Indeed, by Gauss's theorem, considering 
in the usual way the flux out of a flat elementary cylinder containing the surface 
element of mass adS and with its end faces parallel to the surface, 



assuming the constant of gravitation to be unity, or absorbed in a. Instead of a 
surface distribution of zero thickness and (variable) surface density a, it can 
equally be supposed that V is the potential of a distribution of unit density 
everywhere but of variable infinitesimal normal depth f , and we shall then 
have a corresponding equation 

8V\ tdV\ 



Accordingly, if we are given a potential function V of the above form, the thick- 
ness of the surface distribution giving rise to it will be determined by 

477-1 = ^*{L(Q)S\Q)- 
and since L and S are so defined that 



where n is the order of the L function, we obtain finally 

47rf = TD-Sa(2?i+ l)MN. 

Reversing the process, if instead we are given a surface distribution of unit 
density and depth of the form 



then the potential at all points of space can be written down at once as 
F(outside) = S-^L L^S 

F(inside) = X ^ 
and the value of the potential on the surface is 



Liouville's formula 

We can now establish a certain important relation of which we shall have 
need later. 

If we are given a surface distribution 



LIOUVILLE'S FORMULA 
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its potential at a general point P of space can be written down directly. In 
fact, if Q denotes the distance of the element dS from P, then 



and this must be equal to one or other of expressions (6) according as P is outside 
or inside the surface. Since, however, the fa may be regarded as arbitrary, we 
can suppose them all zero except a particular one, and the equality of the two 
forms for the potential becomes 



-- - 

Zn + 1 



L(0)SMN 



if P is outside the ellipsoid, and 



= - 2 *~-~S(0)LMN 

if P is inside. 

If now P is allowed to become a point P of the surface A = 0, the relation 
reduces to Liouville's formula, namely 



where Q now denotes the distance of dS from the point P of the surface itself. 

Use of constant angular velocity for a freely rotating mass in stability 
considerations 

We return now to the question raised in Chapter II (p. 27), and show that 
for certain disturbances or deformations of a freely rotating ellipsoidal liquid 
mass, it is permissible to examine the stability as if to remained constant. That 
is, the equations can be set up by using axes rotating with constant w equal 
to the value in the equilibrium configuration. 

For, consider a surface distribution of (small) depth f measured normal to 
the surface of the ellipsoid (f is of course a negative quantity at some parts of 
the surface). In the case that will concern us 
most, f/to- is a third-order surface harmonic. 
The change in moment of inertia of the system, 
81, can be found simply by calculating the 
moment of inertia of the surface layer, and this is 



Undisturbed surface 



81= f [\$+hco*j)*dhdS 

JEJO 




Fi 14 



where g is the distance of the surface element 

dS of the ellipsoid (E) from the axis of rotation 

0z, h is the elevation above the surface of an 

element of volume dhdS, and </> is the angle between the radial direction, which 

is that of , and the normal at dS. 
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Thus 87 



f (\$ 2 
JK Jo 



= f 
JE 



Now g 2 = x 2 . + y 2 , and on the surface of the ellipsoid this is expressible as a linear 
sum of second-order surface harmonics. For, in the solid expression 

( r 2 7y 2 2 2 \ 
Sr. + * + H- 

6 can be chosen so that the function is harmonic, namely, so that 



The quantity g 2 is then expressible in terms of second-order solid ellipsoidal har- 
monics, that is, as a sum TA i L i M i N i for each term of which n = 2. On the surface 

it therefore reduces to ^LoL i M i N l since A is constant. Accordingly, $ 2 d8 = 0, 

JK 

since the integrand is the product of second-order and third-order harmonics 
together with the factor m. Thus 87 contains terms involving only f 2 and higher 
powers, and is therefore of second order of smallness, which may be indicated Jby 

Writing / = /,+ 8 2 /. 

But since the total angular momentum is invariable 8(Ioj) = and hence 



and the change in angular velocity is therefore also of second order. 

Now in considering secular stability we are concerned with the second-order 
variation of W + M 2 /2I' t the first-order change is already zero since the deforma- 
tions are from an equilibrium form. We thus have 

IT) = 8a w ~ ^p Sa/ to second order 



correct to the second order of smallness. Hence the conditions will be exactly 
the same as for the function W - \u?I with o> treated as constant. 

The same will evidently also be true in discussing ordinary stability if the 
deformation involved in the oscillations is expressible by third-order harmonics 
or higher. For the equations of small motion consist of terms of first order of 
smallness, and hence any additional terms that would arise through SOD, which 
has been seen is necessarily second order, can be neglected as of the second order 
of smallness. Thus constant angular velocity may again be assumed. 

We may notice that these results will also apply to certain second-order 
harmonic deformations provided they are such as not to affect the moment of 
inertia to the first order of small quantities. 
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Calculation of surface gravity on the ellipsoid 

We next establish the important result that on the surface of an equilibrium 
ellipsoidal figure the value of apparent gravity, g, at any point is inversely 
proportional to the perpendicular distance from the centre to the tangent plane 
at the point. That is, g-m is constant over the surface. 

To prove this, we have that at points within the ellipsoid the total mechanical 
potential, that is, gravitational plus centrifugal, may be written 

U = - 7r(oz 2 + fy 2 + yz 2 ) + i^ 2 (x 2 + y 2 ) + constant, 

where a,)3,y are defined by equations (4) of Chapter IV, and we have taken 
Gp = 1 for convenience, as is permissible. The components of apparent gravity 
within and on the surface, measured parallel to the coordinate axes, are 

Bll_ 8U dU 
'fa 9 Hy' dz' 

But in equilibrium g is along the inward normal to the ellipsoid, whose direction 

cosines are 

TUX my TJTZ 

~a 2 ' " 6*"' ~c 2 ' 
Hence considering the third component we have 

-mz 8U 



Accordingly the product gw is constant and is given by 



The result can be expressed more briefly, and also in a form specially suited 
to our purpose, in terms of Lame* functions. For, the first-order (n 1) Lame 
function associated with the third axis c of the ellipsoid, namely, 



= L 1 , say, 
has corresponding ^-function given by 



These two functions L^ and 8^ are of course functions of A, but on the surface 
A = we have at once that 



but it will be convenient to write this result for brevity simply as 



wherein it is to be understood that the Lame functions are to be evaluated on 
the surface of the ellipsoid A = 0, and the product of the axes has been taken 
as unity, that is, abc = 1. 
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Calculation of the coefficients of stability of an ellipsoidal configuration 

. With the results of this chapter available we can now obtain the change of 
potential energy (gravitational and centrifugal) of a rotating ellipsoidal mass of 
liquid when its free surface is subjected to a general small deformation. 

The undisturbed form of the surface has equation A = in ellipsoidal co- 
ordinates, and any point of it has coordinates (/x, v). The motion is a pure rota- 
tion round Oz. We now suppose the surface to undergo an infinitesimal continuous 
displacement without change of total volume of figure, and we denote by the 
distance from a point of the original ellipsoidal surface to the deformed surface 
measured along the normal to the ellipsoid. As before, we denote by g the 
resultant force at (/x,v) due to the gravitational attraction and the centrifugal 
force. Also, we denote by dS an element of surface area surrounding the point (/x, v). 
The total force acting on a particle of matter in the deformed configuration 
may be regarded as made up of three contributions, thus: 
(i) the attraction of the original equilibrium ellipsoid, 

(ii) the attraction of the surface layer formed by the difference, plus and 
minus, between the deformed surface and the ellipsoid, 

(iii) the centrifugal force. 

If dm K denotes an element of mass of the ellipsoid, and Q its distance from a 
given arbitrary point P } the total mechanical potential at P due to (i) and (iii) is 

(8) 

where r is the distance of P from the axis of rotation. The product of the constant 
of gravitation and the density, Gp, may be taken as unity. 

Since the ellipsoid is, by hypothesis, an equilibrium form, V will be constant 
over its surface and equal to ^, say. Also, the gradient of V will be the total 
force that would act on a unit mass at P at rest in the rotating frame due to 
gravitation plus rotation. If P is at a small distance h along the normal to the 
equilibrium surface at a point P lt then to the first order of small quantities 

and here g may be calculated anywhere within first-order distance of P or P l 
without introducing first-order error in V(P). 

If now dm and dm' are any two elements of mass at distance A apart, after 
the deformation has been applied, the potential energy due to gravitation will 

1)6 -iffA^dmdw' 

and the potential energy due to centrifugal force will be 



the integrations being extended to every element of mass of the deformed body, 
and the sign conventions taken in both cases so that the force is minus the 
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gradient of the potential. Thus, the total potential 7, say, is given by 
2U 



= - f JA-Mmctei'- a>* (r*dm 



If we denote by U the value of U in the undisturbed configuration, then stability 
depends on the function U-U . To evaluate this in suitable form, let us denote 
by dm E ,dm' K elements of mass of the original ellipsoid, and by dm L ,dm' L elements 
of mass of the surface layer. Since the whole mass, and volume, are unchanged, 

we have always / / 

dm L = = dm' L , 



\&- l dmdm' \ 



and 

Hence with this subdivision of the total mass we have 



-2C7 = \&- l dm JS dm' K + \&- l dm E dm' L + \br l dm' K dm L + \br l dm L dm' L 

+ oj 2 I r 2 dm E + ou 2 I r 2 dm L 

and 2C7 = kr l dm K dm' K + cu 2 r 2 dm E , 

where in each integral A denotes the distance between the elements of mass con- 
cerned. But obviously 



and hence 



ff -1 ' ff ff 

J J A ldmK<im L = J J ldm * dm L = J J A l dm E d?n L , 

f f rr r * 

-2(tf-y o ) = 2JJ A-irfm^^ + JJ A idm^m^ + ai'Jr-d^ 

r rr 

= 2 Frfm L + br l dm L dtri L 

in virtue of equation (8). 

To reduce this further, consider an arbitrary element of area dS of the surface 
of the ellipsoid. The normals at its boundary generate a cylinder with dS as 
base, and the volume element dm L may be defined by taking a thin slice of this 
cylinder bounded by planes parallel to dS at distances h and h + dh above the 
surface. We may then write 

dm L = dSdh, dm' L = dS'dh', 

and h and h' lie between and f (which takes positive and negative values), 
where is itself a small quantity. We thus have 

r lm r 



Vdm L = V Q dm L + ghd 

= V Q \dm L + \gdS\hdh 
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correct to the second order of small quantities. 

- l dhdSdh'dS', 



Again ((Q- l dm L dm f L = H(fQ- 



and since f and ' are first-order small quantities, Q may be taken as the distance 
between the surface elements dS, dS' of the ellipsoid without affecting the accuracy 
as far as second order. Thus 



and we have, correct to the second order, 

(10) 



Let us now assume that /TCT, which is a function of position on the surface 
of the ellipsoid, consists of a convergent infinite series of the form 



wherein the A t 's are (infinitesimal) constant coefficients. These A i may be 
regarded as the coordinates of the dynamical system, but now infinite in number. 
Thus if (/z, v) is a point within dS, and (/*' ', v) a point within dS', we have 



where the dashes indicate that the functions are calculated at (//, v'). 
Thus p = S Al<m 2 M*N* + 2 S 

i i 

and hence 



fsr^dS = X-4J f^m 2 J/?#?d,Sf + 2 2 S -4^A: ( 

= tTrLj^X^? L-Jf J^f diSf + j-w/^ S S -4^A; 

J ifk J 



making use of the fact that gw = ^nL l S l and is constant over the surface, so 
that it may be taken outside the sign of integration. But every term in the 
second double sum is zero on account of the orthogonality property of the surface 
harmonics, so that finally 



Inserting the values of f and f in the second integral on the right of (10) gives 
2 S A^ 

i.k 
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and here the summations are over all i and k, equal or not. But by Liouville's 
formula 

4ir -- L k (0)S k (0) Jf fc i 



where n is the order of the Lame function L k . The double integral accordingly 
becomes 






In this, in any term for which i^lc the integral vanishes, and the expression 
reduces to the single summation 



Accordingly we arrive at the final expression for the change of potential energy, 
namely, 

U-U = 2L l S l --L i SA\MlN\d8. (11) 



The form of the displaced free surface may be regarded as defined by the 
(small) coefficients A iy which correspond to the q t of ordinary dynamical systems 
involving a finite number of degrees of freedom. Here, however, the A^s are 
infinite in number, there being one corresponding to each Lame function, so 
that there are 2n+l such coordinates associated, one to one, with the 2n+\ 
surface harmonics M i N i of a given order n. We have accordingly obtained the 
expression of U as far as the second order in the coordinates. First-order terms 
do not occur, of course, since the displacement is from an equilibrium configura- 
tion. To each coordinate A i there corresponds a coefficient of stability 



and in this coefficient the factor 4?r \'mM*N\d8 is essentially positive whatever 

the ratios of the axes a : b : c. Hence in order that a given spheroidal or ellip- 
soidal configuration shall be a form of bifurcation it is necessary that one or 
more of these coefficients shall vanish for the ratios of axes involved. Thus 
secular stability is seen to depend on the equation 



which, since A = 0, constitutes a (transcendental) relation between the ratios of 
the axes a : b : c. The discussion of its possible roots, first when a = 6, and 
second for the Jacobi ellipsoids, forms the subject of the two following chapters. 
In the subsequent work it will be convenient to ignore the essentially positive 
factors and call simply J 1 ^ the coefficient of stability. 
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Surface displacements given by first-order harmonics 

It will be recalled that the Lame" function L t occurring in the coefficients of 
stability is ^/(A + c 2 ) (before A is put equal to zero). Accordingly the third of the 
three coefficients of stability corresponding to n = I is always zero. The reason 
for this is because the equilibrium is neutral for any small translation of the 
ellipsoid as a whole parallel to the axis of rotation. To see this, suppose the 
ellipsoid displaced by Sc parallel to 0z, then the equations of the free surface 
before and after when expressed in the same coordinate system are 

x 2 y 2 z 2 = l 

x 2 y 2 (z-8c) 2 , 

and + ' +!_' = i. 

a 2 b 2 c 2 

If, as usual, we denote by the normal displacement of a point P r of the second 
surface from a point P(x, y, z) of the original surface, the coordinates of P 1 are 

(vrX .. TCT?/ . TET< 
x -}- , y ~\ - , z -\ 

and this lies on the second surface if 



Retaining only the first-order terms gives 



But ai + J ^ + c4 - -2 (assuming abc = 1 ) , 

and hence = 2 -arz. 

c 

But z is one of the three ellipsoidal harmonics of order 1, and on the surface 

of the ellipsoid 

z 

whence we have 

where the coordinate A is a small quantity proportional to Sc. 

In the same way a small displacement parallel to Ox or Oy also leads to a 
value of f/iir proportional to the corresponding first-order harmonic. Such 
deformations, however, involve a displacement of the centre of gravity of the 
whole mass off the axis of rotation and are accordingly ruled out since this will 
not affect the question of the relative stability of figure of the mass as a whole 
considered quite apart from any possible translatory motion. For this purpose 
the centre of gravity can always be taken at rest on the axis of rotation. 
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THE SECULAR STABILITY OF THE 
MACLAURIN SPHEROIDS 

The determination of the secular stability for all small displacements involves 
finding the points of bifurcation on the Maclaurin series. We may regard this 
series as starting at the spherical form, which has already been shown to be 
thoroughly stable, and developing in the direction of increasing angular mo- 
mentum. The vanishing of a coefficient of stability then requires 

/^ = (A = 0), (1) 

and for the spheroids we have, since a = 6, 



where r 2 = (c 2 + A)/(a 2 - c 2 ). 

Now equation (1) involves only a and c (in effect their ratio) since A is put 
zero, but if regarded as an equation for a 2 it is clear that A can be retained and 
the equation equally well be thought of as one in A -fa 2 . In the following dis- 
cussion it is found more convenient to proceed in this way; that is, A is not put 
zero but instead A -fa 2 is treated as the unknown in the equation. 

It is to be noticed that when a = 6, although the total number of ellipsoidal 
harmonics of order n is still 2n+ 1, as in the general case, there are only n+ 1 
different L-functions (and M -functions). To each of these, however, there 
correspond two different ^-functions, namely cosp<, siii^, except for p = 0, 
when there is only one, namely N = 1. This makes 2n+ 1 harmonics in all, but 
only n+l different coefficients of stability, since these depend only on the L- 
and $-f unctions. Thus of the 2n + I coefficients of stability, 2n are equal in 
pairs, and if any one of these (for which p4=0) should vanish, it means two 
coefficients vanishing simultaneously. 

The spheroid common to the Maclaurin and Jacobi series 

We first show that, when a = 6, the equation 

^2 = ^1^1-^2^2 = (A = 0), (2) 

wherein the L 2 concerned (of the three such second-order Lame* functions) corre- 
sponds to p = 2, n 2, has a unique root and that it corresponds to the point on 
the Maclaurin series at which the Jacobi series branches off. That is, it corre- 
sponds to the point of bifurcation on these series already established when 
certain restricted second-order displacements only are permitted (pp. 45 et seq.). 
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In the present case we have 

2 =(1+T 2 )Z> 4 (1+T 2 ) 2 

= a 2 + A, apart from positive constant factors. 
Whence S 2 =(a 2 + 

Also we have 

S 1 = V(c 2 + A) 

Substituting these expressions in the coefficient of stability and putting A = 0, 
equation (2), written at length, is 

of 00 dX f dX 

C 2 1 . a 4 ! =0, 

and if we now put a = b in this it becomes 

= 0, (3) 



and this is identical with the condition for a Jacobi ellipsoid with a = b put in it 
(p. 40). This establishes therefore not only that. (2) has a root but that it is the 
only root, when a b, for there is only one spheroidal member of the Jacobi sefies. 

The equation JP/ = 

(It is particularly to be noticed that the arguments of the present section refer 
to the general case a>b>c and therefore apply equally to the Maclaurin and 
Jacobi figures.) 

If r 2 is used as the variable, the spherical form 

a = b = c corresponds to 2 2 = oo, that is, ^ 

a c 

to r 2 = oo. Again, the infinitely flattened disc 
spheroid a = 6 = oo, c = corresponds (since A = 0) 
to A + c 2 = 0, that is, to r 2 0. Hence if we regard 
r 2 as gradually decreasing from oo to 0, the corre- 
sponding equilibrium configuration will describe 
the Maclaurin series from the spherical form to 
the infinite circular disc-spheroid of vanishing thickness. 

Consider then the graph of the function F 2 (r 2 ) over this range. Near r 2 = oo, 
its value is (p. 80) approximately ( )A~* and the function is accordingly 
positive and decreasing for large A. Let us denote by r 2 k 2 the unique finite 
root of the equation. Then the graph of F 2 must be of the form shown in Fig. 15. 
That is, it must be positive at r 2 = oo, begin by increasing as r 2 decreases, and 
eventually cross the axis at r 2 = k 2 . (There might, of course, be more than one 
maximum between oo and k 2 , but the point is not relevant to the argument 
to be given.) 
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Now it will certainly establish that F 2 is the first coefficient of stability to 
vanish if it can be shown that for T 2 ^k 2 any other coefficient of stability F t 
exceeds F 2 . (This result, which can in fact be established, is clearly a sufficient 
condition though not a necessary one.) That is, it will suffice if it can be shown 
that il 

for 



that is, if L z Si- - 4-S Y f >0 for 

- 1 



It will be proved that this inequality holds strongly for all r 2 > 0, that is, for all 
A -f a 2 > a 2 c 2 , and this ensures that it holds for r 2 ^ k 2 since this value corresponds 
to the first Jacobi figure and must therefore be a positive quantity. 

To establish this we begin by considering the more general equation in A + a 2 , 

L t S t = 0, (4) 

l l v ' 



where m denotes the order of the Lame function L k . Let us examine the question 
of what roots, if any, this may have for r 2 > 0. 

We notice first that since LS->X~* as A->oo for any order n, then F-+Q as 
T 2 ->oo. Also, it has been shown (p. 70) that all the roots of the rational part 
of the equation L(A + a*) = 

lie between and a 2 c 2 , and that there is a zero at a 2 c 2 only if L actually 
contains the radical ^/(A-hc 2 ) as factor. Hence, if we consider only values of 
A -ha 2 essentially greater than a 2 c 2 (that is, exclude A + a 2 = a 2 c 2 ), the roots 
of F = are the same as of FjK^ = 0; that is, the same as of the equation 

1 (5) 

' V ' 



2n+i'L t 'L 2 k 

If now this equation is differentiated with regard to A, we shall get an equation 
whose roots separate those of (5). The derived equation is 

i t*di s t \_ i atdim 

' ' V } 



dX\2m+l'L k) 
But by the definition of the ^-functions 



- 



, , d 1 S 1 1 

and hence TT L 7-^r = T^-^ > 

dA(2n+l L) L 2 2^ 

Using this result, the first two terms in (6) disappear, and the derived equation 
is simply 



and its roots separate those of (5) or (4). 
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Since L { and therefore also 8 i are essentially positive over the range 

a 2 c 2 < A-f 2 <oo, 

the factor $ t ./L t . cannot vanish, and (7) simply expresses the condition for 

to have a critical value in this range. Hence the values at which this occurs 

separate the roots of (4). 

Differentiation with regard to A may be replaced by differentiation with regard 

dlL 

to u (p. 79) since -^ cannot vanish for A-f 2 >a 2 c 2 , so that (7) becomes 
d\ 

dL t dL k 

Ltjf. -- J-'; ; = U. (o) 

k du l du v ' 

A further differentiation with regard to u will now provide an equation whose 
roots separate those of (8), namely 



d*L k 
"-du* ~ L < du" = - 

But Lame's equation, with u as independent variable, for L t and L k gives 



1 ) A - K k } L k , 






and inserting these in (9) and omitting the factor L t L k , which cannot vanish for 
A + a 2 > a 2 c 2 , we arrive finally at the simple linear equation 



= K t -K k . (10) 

Since this is of the first degree it obviously can have at most one root exceeding 
a 2 c 2 . Hence (8) can have at most two such roots, and (4) can have at most 
three such roots, that is, values of A + a 2 exceeding a 2 c 2 . 

We notice also that if oo turns out to be a root of (4), as in fact will be shown 
to be so, it would count as one of the three at most whose possible existence is 
demonstrated by the above argument. This would then mean that (4) could 
have at most two finite roots > a 2 - c 2 . We make important use of this result in 
case (ii) below. 

Having established this for the general equation F = we return to the 
original form of the equation (1), viz. F t = 0. The discussion requires considera- 
tion of three cases, as follows: 



(i) Suppose L t contains ^(A + c 2 ) as factor. 
We may then write L t = V( A + c2 )<W/i 

where <(A) = 1, ( A + a2 ) ( A + & 2 ) or 



and/^(A) is a rational polynomial in A in which the highest power occurring has 
positive coefficient. 
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Now it has been shown that f i = 0, as an equation in A + a 2 , has all its roots 

real and lying between and a 2 c 2 . Hence -~ = also has all its roots real and 

df 

in this same interval, and hence for A + a 2 increasing from a 2 -c 2 to oo, -~ has 

dA, 

constant sign, and therefore f t is always increasing, since /^ (oo)= +00. The same 

also holds for both the radical factors ^/(A + a 2 ) and <J(X + b 2 ), and hence in the 

wide sense for any possible < 

Hence in the present case 



LJI* = LJ1* 
is always increasing with A in the range a 2 c 2 < A + a 2 < oo, which means that 



and hence, by (5) and (6), that 



Accordingly FJLl is always decreasing t * A increases, and since it -> + as A->oo, 
it must be always positive for finite A -I- a 2 exceeding a 2 c 2 . We thus have the 
following result (for r 2 > 0) : 

// L t contains <J(X + c 2 ) as factor the equation F t = has no root. 

We notice that in the case a = 6, the radical ^(A + c 2 ) being a factor means 
that T is a factor of L when expressed in the form 

L= (l+ T 2)iP/)lM(l+ T 2)ii. 

But clearly L has r as factor whenever p + n is odd, and hence in the case of 
spheroids for there to be a possibility of F t = possessing a root, the L-function 
must be such that p + n is even. 



(ii) Suppose L t does not contain ^(X + c 2 ) as factor. 
Assuming n ^ 2, we have for A large 



ll -^ ii --- for A<oo 

3 l l 2n+l l i \3 2n+l/ 

and = for A = oo. 

For A + a 2 = a 2 c 2 , since L l = J(X + c 2 ) we have L 1 S l = 0, whereas, on the 
present hypothesis, L^ does not vanish for A -he 2 = 0. But since L t = has no 
root ^a 2 c 2 , the product L^ has the same sign for all values of A-fa 2 >a 2 c 2 , 
and near infinity this sign is plus. Hence for A + a 2 = a 2 c 2 we have 
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The function F i therefore changes sign for some finite A-f-a 2 exceeding a 2 c 2 , 
and the number of roots between a 2 c 2 and oo (not counting oo as a root) must 
therefore be odd. Infinity itself is also a root. But it was shown above that this 
equation can have at most two finite roots, and since now it is necessary that 
the number must be odd, it follows that the actual number of finite roots must 
be just one. 

Hence every L t for spheroids for which p + n is even gives rise to an equation 
F t = possessing a single root A + a 2 > a 2 - c 2 . That is, for some positive finite 
value of A + c 2 (before the infinite disc spheroid is reached) the coefficient of 
stability corresponding to every L i for which p + n is even will vanish and change 
sign. Thus secular instability will enter at some stage for each of the corresponding 
harmonic displacements. 

It will be shown later that the particular L t for which the spheroid first becomes 
unstable is L 2 = 1 +r 2 ; that is, the second order harmonic corresponding to n = 2, 
p = 2. But before proceeding to show this we consider the remaining possible 
form for L+. 



(iii) Suppose n = 1 so that L i = J(\ + a 2 ) or ^/(A + 6 2 ). 

It has already been shown (p. 90) that the displacements corresponding to first- 
order harmonics consist of (small) trans' itions of the ellipsoid parallel to the co- 
ordinate axes. For the third of these, vis. ^/(A + c 2 ), the coefficient of stability is 
zero identically, and the equilibrium therefore neutral for the corresponding 
displacement, as already explained. 

Where the first two are concerned, it is clear since a>b>c that both the ratios 



are always decreasing as A -fa 2 increases from a 2 c 2 to oo, and their final values 
are both unity. Hence they are always positive, and equation (1) accordingly 
has no finite root >a 2 c 2 , and it follows that the corresponding coefficients of 
stability are always negative. 

But, as already mentioned in Chapter VI, such deformations are ruled out 
since they involve a displacement of the centre of gravity off the axis of rotation, 
which would indicate a continued general translatory motion of the whole mass 
within the adopted frame of reference, whereas it is necessary to consider only 
motion relative to axes fixed at the centre of gravity. 

As far as the Maclaurin spheroids are concerned, (i), (ii), and (iii) above cover 
all possible cases. Before proceeding it may be convenient to summarize the 
results of these foregoing sections on the number of roots of the equation 
= lying between a 2 -c 2 and oo. 



(i) // L t is divisible by ^(A + c 2 ). root; L^ - - - - L^ > always. 

~~ 



(ii) // LI is not divisible by ^(A-f c 2 ) and n^2. 1 finite root; 
L 1 S 1 
finite r 2 . 



-L 1 S 1 - - - L i 8 i at first (r 2 = oo) positive, but becoming negative beyond some 
o ~~ 
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(iii) // L t is not divisible by ^/(A + c 2 ) and n = 1. root; and the coefficient of 
stability always negative. But this type of deformation is ruled out as it displaces 
the centre of gravity off the axis of rotation. 

It may again be mentioned that these results hold generally for all a > b > c, and 
therefore apply equally to the discussions of the Maclaurin series and the Jacobi 
series. We proceed now to establish the following result for the Maclaurin series. 

The first coefficient of stability to vanish is that corresponding to L 2 (p = 2). 

As the mass gradually flattens from the spherical form and describes the 
Maclaurin series, the first coefficient of stability to vanish will be shown to be 
that corresponding to 



i.e. 2> = 2, n = 2 
= 1 -f r 2 , apart from a constant factor. 

As far as L L is concerned, it is only necessary to consider p + n even, since it 
has been proved that there can be no root otherwise. We wish, therefore, to show 
that the root r 2 = 1c 2 (corresponding to the spheroidal member of the Jacobi 
series) of the equation 

O (p-2) 



is greater than (since r 2 decreases from oo as the spheroid flattens) the root 

T = *, Of 



For this purpose, as has been seen, it will certainly suffice if it can be proved that 

3 LlSl ~ >2n + 1 L(8i > 3 
or that Ji 2 S 2 

and this has already been seen to be equivalent to proving that L t /L 2 is always 
increasing with r 2 . 

To prove this it is necessary to consider three separate forms for L t covering 
all cases. It is important to notice that we are again able to establish sufficient 
conditions; that is, the above inequality in fact holds strongly and we are able 
to take advantage of this in the proof. 

Case (a). Suppose p^ 2. Then 

LI/LI = (1 +T 2 )*<"- 2) /)"+"( 1 +r 2 )", p + n is even. 



In this expression, the first factor (l+ T 2 ) io >~ 2) either reduces to 1 or else is 
always increasing with r 2 , according as p = 2 or p>2. Again, since (l+r 2 ) n is 
a polynomial in r 2 with all its coefficients positive, so also must be Z) pfn (l + r 2 ) n , 
and it is therefore always an increasing function of r 2 . Hence L t /L 2 is always 
increasing, and accordingly it follows that for L L of the present type the coefficient 
of stability corresponding to L 2 must vanish first. 
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Case (b). Suppose p = 0; then n = 2j, where j is an integer ^ 1 (since p + n 
is necessarily even). 
We now have 



(This form includes the remaining Lame* function L = 1 + 3r 2 of order 2, given 
by 7> 2 (l+r 2 ) 2 . The second function L = r^l + r 2 ) contains T as factor and its 
corresponding coefficient of stability cannot vanish.) 

In this case Z^/ 2 ls always positive since both the numerator and denominator 
are always positive. Tf therefore it can be shown that its differential coefficient 
is always positive it would follow that it is always increasing with T. The 
logarithmic differential coefficient is 



V r 2 +l' 
and what we wish to show is equivalent to proving that 

E(r) = (r 2 + 1 ) Z^'+V + 1 ) 2j - 2r D^(r 2 + I ) > 0. 

Now (T2+1)2 ; = s, !( |^ y] r- <t-0,l,...,2fl, 

so that 

2?"i *>M 

/)2j/_2 , 1 \2; = y _____ __ *J -*" _ -2k-2j (1. _ ,/ ,/ , 1 f),}} 

( } tt!(2j-t)!(2t-2j)! (K-J,J + L,...,-J) 

= 2 -V 2 *- 2 ', say, 

k 

and it is seen that all the coefficients A k are essentially positive. We thus have 
E(r) = 



and it will be shown that this polynomial in T has all its coefficients positive. 
The coefficient of r 2 *- 2 ^ 1 is found to be 



Since A k > and k ^ j, the only one of these coefficients that might be negative 
is that for which k = j, and this one is 2A } + 2A j+l , since for k >j all the factors 
multiplying A k and A k+l are positive. But this one coefficient will also be positive 
if Aj +1 > Aj. That is if 



that is, if (2j + 2)(2j+l) > 2(j+l)lj, 

or if j(2j+l)>l, 

and since j^ 1, this is certainly satisfied. 

Thus all the coefficients in E(r) are positive, and therefore LJL 2 is positive 
and always increasing with r 2 . Hence the coefficient of stability corresponding 
to L 2 = 1 + r 2 vanishes before any of those corresponding to L t of the present type 
as the mass gradually evolves along the spheroidal series from the spherical form. 
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Case (c). Lastly, we have p = 1, n = 2j+ 1 (since ^? + w is even). It is to be 
noticed that j > 1 since j = would correspond to a first-order harmonic, and 
these are already excluded. 

In the present case 

fy = ( 1 + r 2 )* /) 2 '+ 2 ( 1 + r 2 

and hence L /L 2 



This will be always a positive increasing function of r if its logarithmic differential 
coefficient is always positive, that is, if 



is always positive, and this will be so if 

E( r ) = (r 2 + 1 ) Z> 2 '+ 3 (r 2 + 1 ) 2i+l - TZ) 2 >'+ 2 (T 2 + 1 J 2 ^ 1 > 0. 
Proceeding in a similar way as in case (6), we have 

t 2 , l\2}+l = y _______ ( 2 ilil^A!_____ T 2k-2j-2 

( ' T 



, say, 

k 

and here all the coefficients B k are positive. Thus we have 

E(r) = (r 2 +l)X(2A;-2j- 
and in this the coefficient of r 2 *-^'" 1 is 



This might possibly be negative if 2k < 2J + 3, that is, if k = j+ 1, in which event 
it becomes 2/? ; . 42 J3 j+1 . But even this will be positive if 



or if 2j(2j + 3) > 1, 

and this is certainly satisfied even for j 1. 

Accordingly it follows that the coefficient of stability corresponding to 
Z/ 2 = 1 H- r 2 vanishes before that corresponding to any L t of the present form. 

Taking these three cases together we may therefore conclude that as r 2 decreases 
from oo, the first coefficient of stability to vanish is that associated with the 
second-order Lame* function L 2 = a 2 + X. Thus secular stability of the Maclaurin 
series is lost at this point. It is to be remembered, however, that the present 
result has been reached from consideration of the function TT ico 2 /, and this 
corresponds to stability relative to a frame of reference rotating with constant 
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angular velocity. Since the harmonics associated with this L 2 are the second- 
order polynomials x 2 y 2 and xy, it is possible that o> is changed to the first 
order of small quantities as a result of the corresponding surface deformations. 
Accordingly it is necessary to establish that the instability arrived at here is a 
true instability (see p. 26). 

As far as xy is concerned, it is readily shown that a deformation f /m involving 
the surface harmonic MN to which xy reduces on the ellipsoid will not affect 
the moment of inertia / to the first order. For it has been seen (p. 83) that the 
first-order change in / is 

!/= | (x 2 + y 2 )dS 
JE 

and on the surface x 2 + y 2 can be expressed in terms of solid harmonics in the form 



where h, k, and I are constants. Accordingly, since xy is a harmonic independent 
of this, the surface integral giving S l I vanishes. Thus / is changed only to second 
order and the instability indicated by the vanishing of the coefficient of stability 
associated with xy must be a true one. 

On the other hand the surface deformation corresponding to x 2 y 2 will produce 
in Sj 7 a term of the form 

(uI^MldS 

and the moment of inertia is therefore changed to the first order. But the sur- 
face harmonics corresponding to x 2 y 2 and y 2 z 2 are precisely those involved 
in the restricted ellipsoidal displacements considered in Chapter IV (p. 45). 
It is easily verified that the normal displacement f required to deform the 
standard ellipsoid 222 

-- + ?- + -= 1 
a 2 + & 2 c 2 

into a slightly different one with common principal axes, say 

y =1 with 



(a+Sa) 2 

when made harmonic, is given by an expression of the form 

f /tnr = p(x 2 - y 2 ) + q(y 2 - z 2 ) + r, 

where p,q,r are small constants vanishing if Sa, S6, 8c all vanish. Accordingly 
the argument given in Chapter IV, which is based on the function W + H 2 /2I 
appropriate to systems with changing angular velocity, can be adduced to show 
that the instability corresponding to x 2 y 2 is also a true instability. 

Thus it may finally be concluded that the Maclaurin spheroid is secularly 
stable, and therefore also ordinarily stable, for all deformations (representable 
by ellipsoidal surface harmonic expansions) provided its eccentricity of meridian 
section is less than that of the spheroid of bifurcation, that is, for e< 0-8 127. 
For values greater than this the spheroid is secularly unstable. 
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The question of its ordinary stability, which cannot of course be settled from 
the present considerations, has been studied by Cartan, who has shown that the 
spheroidal form remains ordinarily stable for all possible small oscillations (repre- 
sentable by ellipsoidal surface harmonics) provided the eccentricity of section is 
less than 0-9529. 

From a practical standpoint, however, the secularly unstable spheroids are of 
little interest, for if we are regarding the mass as undergoing evolution owing 
to a gradual increase of angular momentum, then when the value exceeds that 
corresponding to the first point of bifurcation, the mass would proceed to evolve 
along the Jacobi series, which (as will be shown in the next chapter) begins by 
being thoroughly stable. The departure from axial symmetry thereby involved 
may be regarded as brought about by some unspecified small external disturbance 
which can always be regarded as operative in actual physical systems. 

The Maclaurin spheroids beyond the first form of bifurcation 

If the series of spheroids is followed beyond the point at which the Jacobi 
series bifurcates, and regarded as subjected at any stage only to harmonic 
deformations of a given order n, then it can be proved that the first coefficient 
of stability to vanish corresponds always to the highest value of the index p t 
that is to p n. In other words, the first coefficient of stability to vanish is 
that corresponding to the Lame function 



It may further be proved that the coefficients of stability corresponding to 
p = n 2, p = n 4, . . . vanish successively as the series continues to be described 
in the direction of increasing eccentricity and angular momentum. But it does 
not hold that all such instabilities of a given order n are reached before the first 
of those entering from harmonics of the higher orders n + 1 , n -4- 2, 

Since, however, all such configurations would already have become secularly 
unstable for displacements corresponding to L(n = 2, p 2) they are of no 
physical application and could not come into existence through any natural 
evolution of a liquid mass. If a system possessed an amount of angular momentum 
appropriate to any such spheroidal form, it would in fact settle down through 
the action of internal friction to the corresponding equilibrium configuration on 
the Jacobi series provided such configuration (with this angular momentum) is 
itself secularly stable. Accordingly we proceed next to consider the secular 
stability of the ellipsoidal forms. 



Chapter VHI 

THE SECULAR STABILITY OF THE 
JACOBI ELLIPSOIDS 

The object of the discussion of the present chapter is to discover at what stages 
on the Jacobi series forms of bifurcation are reached. The main result from a 
practical standpoint is that the first coefficient of stability to vanish is that 
corresponding to a certain third-order harmonic. 

We have now the general case a>b>c. As usual we put 

A + c 2 2 ,, . A + a 2 2 f 

-_- = r 2 , so that - = T 2 + 1. 
a 2 - c 2 a 2 - c 2 

It is found convenient to represent description of the Jacobi series, in the direc- 
tion of increasing angular momentum, by making a 2 6 2 vary from to a 2 c 2 . 
The value a 2 6 2 = clearly corresponds to the initial spheroidal form, while 
a 2 6 2 = a 2 c 2 , that is 6 = c, corresponds to the infinitely elongated final Jacobi 
form a==oo,& = c = (Table II, p. 41). To each value of a 2 6 2 in this range 
there will correspond a certain Jacobi figure, and hence a certain value of r 

a 2 b 2 
(since A is always zero). That is, for given 2 ^ between and 1 there is a 

tt ~~ C 

unique Jacobi ellipsoid, and therefore definite values of the ratios a : b : c, or of 
a, b and c themselves if we suppose abc = 1 always. In terms of r the axes of 
this ellipsoid will be 



so that a:b:c = J(r*+ 1) : 



The condition for a Jacobi ellipsoid 

This condition has already been found from first principles in Chapter IV 
(p. 38) but it can be derived quite simply in a form more suitable for our purpose 
in the following way. 

Suppose we have an equilibrium Jacobi figure E, and it is given a rigid body 
rotation through an infinitesimal angle 6 round the z-axis, so that referred to 
the axes of E it is a new ellipsoid E' of the same size and shape as E. If the 
normal distance of points of E' from E is denoted by f , it is easily found that 

= *(i- 
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where M 2 N 2 is that one of the five second-order surface harmonics given by 

M 2 = VKfi + a) (/* + 6 2 )], N 2 = V[(" + 2 ) (" + 6 1 )]- 

For a displacement of the present kind it is evident that equilibrium must be 
neutral, for E' is simply a Jacobi form referred to different axes. Hence the 
corresponding coefficient of stability must vanish. That is 

iL^-iV^O (A = 0), (1) 

and this must accordingly be the condition for a Jacobi ellipsoid. It is easily 
verified that it is equivalent to the condition already found in Chapter IV. 

a 2 6 2 

It then follows from Table II that for a given value of - 2 - 2 between and 1, 

a c 

the equation has just one real solution r = j 2 , say. 

Further, if b->c for a configuration on the Jacobi series the axes become in 

a 2 -b 2 
the ratios a : b : c = 1 : : 0, and hence r = 0. Accordingly as 2 -- jj"^' ^ e 

tt ~~ C 

root j 2 -> 0. 

Condition for a point of bifurcation 

For this it is necessary that another coefficient of stability vanishes. That is, 
in addition to (1), we must have 



Hence for a point of bifurcation on a Jacobi figure we require simultaneously 

L& (A = 0), (3) 



where t = V(A + c 2 ), and L 2 

Now it has already been proved (p. 94) that equation (2) cannot have any real 
root if L t is divisible by ^/(A + c 2 ), and that if it is not so divisible then it has 
one real finite root. Let us suppose that when L t is not divisible by ^(A + c 2 ) 
the root of (2) is denoted by r = /^. Then for both of equations (3) to be satis- 

fied simultaneously requires . _ . 

J2 K i- 

It will now be shown that for a giveri order n, there must be at least one L { 
for which a solution of (3) exists. To do this let us consider first the situation 
at the two ends of the range, viz. a = b, and b = c. 

(i) a = b. This condition alone means that the ellipsoid is always on the 
Maclaurin series; r oo corresponds to the sphere and r to the infinite disc. 
The root r J 2 of (1), taken in conjunction with a = 6, gives the first member 
of the Jacobi series, and the vanishing of any other coefficient of stability (on 
the spheroidal series) has been shown to occur for a value of r smaller than j 2 , 
that is, at a more flattened form. Accordingly, at this end of the range we have 

j 2 > k { for a = b. 
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(ii) b = c. It has been seen (p. 63) that in this case the function L i is given by 

I* = (-!) jD*+(* a -l) (p = 0,l,2,...,n), 
where t 2 = (a 2 + A)/(a 2 c 2 ). Also we have 



Hence if #H=0, equation (2) is satisfied for t = 1, or T = 0, since then both Z^ 
and Li vanish separately. Consequently in such cases the root Jc t is zero, and 

therefore we have 

J2 ki, (i = 2) for b = c. 

If, however, # = 0, we have 



and t = 1 is no longer a solution of (2) since Z t - does not vanish. Accordingly 
the solution now is not T = but must be some finite positive value. Hence 
in this case we must have 



Hence for that particular function L i of given order n that reduces when b c 
to the .L-function for which p = 0, the equation j 2 = k t is satisfied at least once, 
and possibly an odd number of times, for some finite value, or values, of T. Also, 
when b = c there is only one Z-fimction (for each given order n) for which p = 0; 
that is, only one of the 2n + 1 general jL-functions of given order reduces for 
6 = c to that particular L t for which p = 0. 

It is thus established, by combining the results of these two cases, that the 
coefficient of stability corresponding to this particular general .L-function certainly 
vanishes for some intermediate a : b : c. It will be shown in due course that 
it is in fact the only one that can do so. To distinguish it from the remaining 
2n coefficients of stability, which turn out to be always positive, it is termed 'the 
characteristic coefficient of stability of order n\ 

The equation ^Z^ - ^ - \LjSj = cannot have a root if L/ has j/(X+i 2 ) 
as a factor 

Let us now take the equation 

' < 4 > 



and in order to find the forms of L t for which it cannot have a root let us consider 
in what cases LJL 2 is always increasing. 

Omitting Z-functions divisible by ^(A + c 2 ) (for which it has already been seen 
the equation can have no root), the functions L i can take four different possible 



THE EQUATION FORMS OF L t 105 

forms, and written in factorized form these are as follows: 

(i) . = (A + a 2 - ai )(A + a 2 - a2 )...(A + a 2 - a J \ for 

(ii) L t = J[(X + a 2 )(A + 6 2 )].(A + a 2 - ai )... (A + a 2 -^^)) n = 2m 
and 

(iii) ^ = V(* + a 2 )-(A + a 2 -i)-MA + a 2 - a , n ) for 

(iv) L i 



(It will be understood that in these expressions the roots are denoted convention- 
ally by !, a 2 , a 3 , ... in each of the forms without implying that the roots are the 
same in each expression.) 

In each of the four expressions we may suppose the factors arranged in 
decreasing order of the a's, that is, so that 



Now it has been shown (p. 70) that for any /^ the roots in A + a 2 are all real and 
positive and less than a 2 c 2 . That is, 

< a < a 2 c 2 for every a. 
Hence A + a 2 - a > A + c 2 ^ 0, 

and therefore all the factors A + a 2 a will be positive increasing quantities as 
A + a 2 increases from a 2 c 2 to oo. 

We have also that for A + a 2 increasing 



increasing quantity decreasing quantity 
increasing quantity for every a. 



Also 

and this will certainly be always increasing provided a > a 2 b 2 . 

Now for the Jacobi ellipsoids L 2 = J[(X + a 2 ) (A + 6 2 )], and hence the ratio L t /L 2 
can take four possible forms corresponding to the above, and these can be written: 

7 2 (Y \ ( \ _L n ;2 



(ii) j- = (A + a 2 - ai ) ... (A + a 2 -,_!) 

and 



odd. 




106 SECULAR STABILITY OF JACOBI ELLIPSOIDS 

In these four products all factors involved, except for the first factor {} in each 
of (i) and (iii), have already been shown to be always positive and increasing. 

Hence in cases (ii) and (iv) the ratio LJL 2 is always increasing, and it is seen 
that these are precisely the cases in which Z^ is divisible by ^/(A + 6 2 ). 

In cases (i) and (iii) LJL 2 will be always increasing provided <x l >a 2 ft 2 . That 
is, (4) can only be satisfied if the greatest of the roots ot lt a 2 , ... is less than a 2 6 2 , 
and this means that all the roots of the rational part of L t must be less than a 2 6 2 . 

Hence we can conclude that if equation (4) is to have a root, L t must not be 
divisible by ^(A + 6 2 ) and also all its zeros (of the rational part) must be less 
than a 2 6 2 . (It is not yet established that it does actually have a root, but 
only that it certainly cannot do so if these conditions are not satisfied by L t .) 

Only one coefficient of stability of given order n is capable of vanishing 

To establish this it is necessary to consider the cases of n odd or even separately. 

(i) Suppose n even. The JO-functions can in the first instance take only one 
of the following four forms : 

4, V[(*+& 2 )(A+c*)]/, r , VP+c 2 )(A+ 2 )]4. Vr(A+ 2 )(A+& 2 )]4> 

where L r denotes a rational polynomial in A. But it was shown in Chapter VII 
that no coefficient of stability can vanish corresponding to an L t that is divisible 
by ^(A + c 2 ), and it has been shown in the foregoing section that it cannot vanish 
if ^ is divisible by ^/(A + 6 2 ). Accordingly, for n even, the last three of the above 
forms are ruled out and only a coefficient of stability corresponding to a wholly 
rational L t can possibly vanish. For these rational functions, /c x = /c 2 = /c 3 = 
(see p. 57), and there are %(n + 2) = m+l such functions. But by Stieltjes' 
theorem, only one of these has all its zeros less than a 2 6 2 , and this is that one 
for which r = m+l = ln+l (in the notation of the proof of this theorem). 

(ii) Suppose n odd. To begin with, the possible forms for L t are now 



where L r denotes a rational polynomial, and again the last three forms are ruled 
out because of the presence of ^/(A + 6 2 ), or ^(A + c 2 ), or both, as a factor. There 
remains only the typo ^(A + a 2 )L r , and since now /c x = i, /c 2 = * 3 = 0, there are 
%(n+l) = m+l of these in all. Again, however, only one of these has all its 
zeros less than a 2 6 2 , namely that for which r = m+ 1 = \(n+ 1). 

Hence in either event, n odd or even, only one coefficient of stability (of the 
2n + 1 in all of a given order n) can vanish. This one has been termed the 
characteristic coefficient of stability of order n. 

This does not of itself establish that the characteristic coefficient does vanish 
for some a : b : c, but the result now follows from the conclusion reached on p. 104, 
above. For it was shown that there is certainly one function L t (of the 2n+ 1 in 
all) for which j 2 = lc i is satisfied, that is, for which the corresponding coefficient of 
stability vanishes. This, together with the result of the present section, finally 
establishes that for any given order only the characteristic coefficient of stability 
can vanish and that it does do so for some finite Jacobi figure. 
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Again, the argument of p. 93 et seq. is available in the present case to show 
that the equation 

3^-2^1^ = < 2 > 

has only a single finite root. Now the initial (spheroidal) Jacobi form has 
a : b : c = 1-197 : 1-197 : 0-698 and therefore r 2 = 0-515. The final Jacobi form 
has a : b : c = oo : : and hence r 2 0. For intermediate values the ratio a : b 
is uniquely settled if the ratio a : c is assigned (Table II, p. 41) so that equation 
(2) is simply an equation for a : c. But the argument of Chapter VII (p. 95), 
applies generally to show that the number of finite roots of (2) is one. Accordingly, 
as the Jacobi series is described in the direction of increasing angular momentum, 
each characteristic coefficient of stability changes sign but once. If displace- 
ments of a given order n only are contemplated, the mass eventually becomes 
secularly unstable for some finite value of a : b : c, and thereafter remains 
secularly unstable. 

It can be established (see below) that instability first enters through the 
characteristic third-order harmonic deformation, and that as the series continues 
to be described, further instabilities enter successively through the characteristic 
harmonic of the fourth order, then of the fifth order, next of the sixth order, 
and so on. If for a given order n we denote the characteristic coefficient of 
stability by C n , then the following result can be proved, establishing that these 
coefficients vanish successively. 

// n > p, then C H > C p . 
If C n = we have 



where L n denotes the Lam6 function corresponding to the characteristic coefficient 
of stability of order n. To establish the result it will be sufficient if it can be 
shown that LJL p , where n>p, is always increasing as A + a 2 increases from 
a 2 c 2 to oo. For if this is so, then 



and hence (from Chapter VII, p. 93) 



Or, in terms of C p and C n , ~ {(C H - C P )/L*} < 0. 
In this, 7/|, is essentially positive, and since 

^ if n>p as A ~ >00> 



this means that C n C p is always positive, and hence that C p must vanish first. 
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To show that LJL p is always increasing we consider three cases as follows. 

(i) Suppose a = b. 

Since the //-function giving the characteristic coefficient of stability has all its 
roots in A + a 2 between and a 2 6 2 , when a becomes equal to b all these roots 
must come into coincidence at A + a 2 = 0, and hence we have 

n even L n = (a 2 + A)* /l , 

n odd L n = V(a 2 + A)...(a a + A)K*- 1 > = (a 2 + A)*" again. 
Hence LJL p == (a 2 + A)<-">, 

and since n>p, this is always increasing as a 2 + A increases from a 2 c 2 to oo. 

(ii) Suppose b = c. 

In this case it has been seen that L n = D n (t 2 l) n 9 and hence 

A, />(*- 1) 
L p j)p(t*-l)p' 

and as A increases, t 2 = (a 2 + A)/(a 2 -c 2 ) increases from 1 to oo. To prove that 
L n /L p increases, if we put t = 1 + 2x, then apart from positive factors, 



LJL p = D[x(x + l) n }ID{xi>(x+ I)**}. (6) 

If x = this has the value n\/pl, and it is easily shown that 
p I D n {x n (x + 1 )"} - n \ D*>{x p (x + 1 ) p } 

contains only terms with positive coefficients if n>p. It is therefore always 
increasing with x ( > 0), and hence, so is LJL jr 

(iii) Suppose b arbitrary and between a and c. 

When u (p. 79) is taken as independent variable instead of A the Lame equations 
for L n and L Jt become 

d du* = {n(n + ] ) ( A + a 2 ) - K' n }L n , (7) 

~du* = 

(the constants K' n , K' p have been so written, with dashes, to allow for A + a 2 being 
taken as variable on the right instead of simply A). Now in order that the ratio 
L n /L p shall not be always increasing it is necessary (as has been shown in Chapter 
VII, p. 94) that the equation in A + a 2 



shall have a root greater than a 2 c 2 , and this can be shown to be impossible 
by the following argument. 
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We have (from p. 79), 



Hence the expression L n -= n - L n - 7 - - 

* du n du 

vanishes for A + a 2 = a 2 -b 2 and for A + a 2 = a 2 -c 2 , since A + 6 2 = and A + c 2 = 
make the radical factor vanish, and L n itself does not contain ^/(A + 6 2 ) or ^(h + c 2 ) 
as a factor. But we have, by means of the differential equations for L n and L p9 



<"> 



and hence i/jL n L p must vanish at least once in the interval 



But L n L p cannot vanish, since all the zeros of both L n and L p are less than 
a 2 b 2 . Hence i/j itself must vanish for some A + a 2 in the interval, and since \fs 
is linear it cannot also vanish elsewhere. 

Accordingly it is established that in all cases L n /L p is always increasing, so 
that C n > C p if n >p, and it follows that the coefficients of stability corresponding 
to harmonics of order 3, 4, 5, ... vanish successively as the Jacobi scries is des- 
cribed in the direction of increasing angular momentum. 

The L-function corresponding to the characteristic coefficient 

To find the characteristic coefficient of any given order n, use can be made of 
the property that when b c the corresponding L-function reduces to D n (t 2 l) n . 

This result can be established independently from the fact that in the general 
case L must be such that all its roots are less than a 2 b 2 . For if we put b = c, 
roots <a 2 b 2 will remain so. But then L reduces to 



L = (t 2 - l)**>7>+*(* 2 - i), (11) 

and unless p = there would be a root at t 2 = 1, i.e. A + a 2 = a 2 c 2 >a 2 b 2 
(p = 1 is excluded since it would give ^/(A + c 2 ) as factor). 

For example, for n = 3, we require to find the general third-order Lame* func- 
tion that reduces when b = c to 

D*(t 2 - 1 ) 3 = t (M 2 - 3) = V( 2 + A) ( A + fa 2 + f c 2 ), 

apart from constant factors. The third-order Lame functions were found on 
page 60, and it is easily verified that the appropriate one of these so reducing is 
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The Jacob! ellipsoid at which the third-order characteristic coefficient of 
stability vanishes has been determined by Darwin,* who finds 

a : b : c = 1-8858 : 0-8150 : 0-6507 
and a> 2 /27rGp = 0-1420, H = 0-3896. (See Table II, p. 41.) 

We may conclude therefore that all ellipsoidal forms less elongated than this are 
thoroughly stable, that is, both secularly and ordinarily stable, and all ellpisoids 
more elongated are secularly unstable at least for third-order harmonic displac- 
ments. 



The pear-shaped figure 

According to the theory of linear series of equilibrium configurations, the 
occurrence of a point of bifurcation on the ellipsoidal series indicates that another 
linear series of equilibrium configurations crosses it at this point. So long as 
the free surface of the liquid mass is restricted to ellipsoidal form, that is, so 
long as l/Ttr depends only on second-order harmonics, no point of bifurcation 
can reveal itself, just as the Maclaurin series is found to be always secularly stable 
if spheroidal forms only are permitted. Moreover, the initial form of the new series 
can be obtained from the ellipsoidal series by applying a small normal deformation 
f/OT depending on the particular surface har- 
monic J/a^ through which instability enters, 
just as in the spheroidal case the initial members 
of the Jacobi series can be obtained by means 
of a small deformation proportional to that 
M 2 N 2 for which the instability enters. 

If /m is taken small but finite, indications 
of the shape of the members of the new series 
can be calculated. In his original Memoir on this 
subject, Poincare gave a sketch of the resulting 
surface, and although the figure is not of course 

axially symmetrical round the longest axis (of the original ellipsoid, for which 
in any case 6=j=c) the resemblance to a pear was so striking that it seemed appro- 
priate to call the new form 'the pear-shaped figure' and the series 'the pear- 
shaped series'. Since then it has also been referred to as 'the piriform figure' 
and sometimes simply 'the pear'. It is of interest, however, that Darwin's 
detailed calculations, made later, both of the axes of the critical ellipsoid and 
of the normal displacements corresponding to the crucial third-order harmonic 
lead to a form of the figure distinctly more elongated than Poincare had supposed 
and bearing rather less resemblance to a pear. Poincar^'s sketch and the form 
as found by Darwin are shown in Figs. 16 and 17, which make clear the extent 
of the resemblance to a pear. But a name is useful in any case for the new series 
and there has been general acceptance and usage of the term pear-shaped. 




Fig. 1G. (After Poincare', 'Figures 
(TEquilibre\ p. 101.) 



* Collected Works, in, p. 312. 
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The stability of the pear-shaped figure 

The problem of the stability of the pear-shaped figures has been studied by 
Liapounoff, Darwin, and Jeans, and has proved to be one of the utmost intricacy. 
The chief difficulty lies in calculating with sufficient accuracy the gravitational 
potential of the pear-shaped form, the value of which is necessary for the expres- 
sion of the constancy of the total mechanical potential (gravitational plus centri- 
fugal) over its boundary surface. The pear-shaped figure is not representable 




M 



O M 

Fig. 17. (After Darwin, 'Scientific Papers', in, p. 314.) 



by any simple closed analytical form, but it has been found possible to examine 
its properties in the immediate neighbourhood of the critical Jacobi form at 
which the piriform series bifurcates, for this can be achieved by regarding the pear- 
shaped series as produced by the addition to the critical ellipsoid of a small 
deformation measured by an infinitesimal parameter in powers of which series 
expansions can be introduced. Thus the equation of the pear-shaped figure can 
be shown to be capable of being represented to sufficient accuracy by an 
expression of the form 



(12) 



where the first term (in brackets) equated to zero is the ellipsoid of bifurcation, 
e is the small parameter representing the degree of (small) departure from the 
point of bifurcation, and P , Q Q , R Q are certain polynomials in x, ?/, z in which 
the terms of highest degree are respectively of order 3, 4, and 5. The problem 
is then to calculate the coefficients in these polynomials appropriate to the 
pear-shaped equilibrium form by means of the condition for relative hydrostatic 



112 SECULAR STABILITY OF JACOBI ELLIPSOIDS 

equilibrium. In the present case this condition is found to take the form 



(13) 

Here V g is the gravitational potential of the mass bounded by (12), and it turns 
out to be necessary to calculate V g as far as terms of the third order in e, that 
is, to e 3 . It is this requirement that renders the problem of such great length 
and complexity. The term e 2 6(o> 2 ) represents the adjustment to the square of 
the angular velocity, o> 2 , associated with the critical Jacobi form, and its value 
is found in the process of calculating the coefficients in P , # , R Q . The factor 6 
is simply a constant of proportionality included here since it is only necessary 
that the total potential shall be constant over the free surface. 



Pear-shaped 




Fig. 18. 

When the coefficients involved in P ,Q ,R Q have been found, it is then com- 
paratively easy to calculate the moment of inertia of the corresponding pear- 
shaped mass as far as terms in e 2 , which is all that is finally necessary. Then 
from the calculated value of the angular velocity co 2 + e 2 8(o> 2 ), the value of the 
angular momentum appropriate to the equilibrium pear-shaped forms in the 
neighbourhood of the point of bifurcation follows at once. According to Jeans's 
calculations the values are as follows: 

o> 2 + e 2 S(o> 2 ) = 2nGpx 0-142(1 + 0-05227e 2 ) (14) 

for the angular velocity of the equilibrium pear-shaped figure, and 

H = // (l-0-06765e 2 ) (15) 

for its angular momentum, where H = 0-3896(6r*Jtf*r*), r 3 = abc. 
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If therefore, as usual, we regard the system as evolving slowly in the direction 
of increasing angular momentum, and represent the Jacobi series by a vertical 
line in a diagram, and the parameter e, representing departure from it, by a 
horizontal coordinate, as in Fig. 18, then it is evident from (15) that the pear- 
shaped series of equilibrium configurations are such that their representative 
curve initially turns downwards from the point of bifurcation (7, as indicated 
in the diagram. Accordingly, by the considerations of Chapter IT (p. 10), the 
pear-shaped series must be secularly unstable initially. This conclusion has also 
been reached independently by LiapounofF and, according to Jeans,* can also 
in fact be established from Darwin's researches once a certain technical error is 
put right. 

What the result means physically is that if the mass were set in motion in 
precisely the state consistent with a pear-shaped figure (in the neighbourhood 
of (7), it would be in equilibrium, and if entirely undisturbed would continue in 
rigid-body rotation in that form. Also, its total angular momentum would be 
slightly less than that associated with the critical Jacobi figure, but its angular 
velocity slightly faster. If, however, any slight disturbance were given to the 
system, the presence of internal frictional forces would cause at least some of 
the amplitudes associated with the disturbance to increase, and the system 
would gradually depart more and more from the pear-shaped form to settle 
down finally in that Jacobi form possessing the same angular momentum as the 
original pear-shaped configuration. It is therefore clear that by no natural 
course of evolution can any of the initial members of the pear-shaped series 
ever come into existence. As H gradually increases the mass will evolve along 
the Jacobi series to C, and if it should increase further the only equilibrium form 
available for the system will be the appropriate Jacobi ellipsoid, which itself 
will at this stage have already become secularly unstable. 

The nature of the further development of the mass will then depend on whether 
these ellipsoids are ordinarily stable or not. If they were ordinarily stable, it 
would mean that any slight disturbance would increase only at a rate dependent 
on the amount of friction present, and this might not involve a rapid departure. 
An example of such a case is provided by the lunar orbit which, though secularly 
unstable, is ordinarily stable, and the rate of departure from the present con- 
figuration is extremely slow because of the very small amount of tidal friction 
available. On the other hand if the Jacobi ellipsoids beyond the point of bifurca- 
tion are ordinarily unstable, as in fact they will be shown to be, the exponential 
factors indicating the instability are independent of friction and do not vanish 
with it. They accordingly may involve dynamical changes from the Jacobi 
form proceeding at an acceleration comparable with any other acceleration in 
the system. The examination of the ordinary stability or otherwise of the Jacobi 
ellipsoids accordingly forms the subject of the following chapter. 

* Problems of Cosmogony, p. 92. 



Chapter IX 

THE ORDINARY STABILITY OF THE 
JACOBI ELLIPSOIDS 

To investigate ordinary stability it is necessary to set up the equations of small 
motion of the liquid and consider the conditions under which their solutions 
are purely periodic or otherwise. 

General equations of small motion of a rotating liquid 

We adopt rectangular axes Ox,Oy,Oz and suppose them rotating round Oz 
with angular velocity a), which may be assumed constant since we intend 
to examine small motion about an equilibrium configuration having angular 
velocity w. The velocity components u, v,w relative to these axes then suffice 
to deal with any small departure. 

The component accelerations, relative to fixed axes coinciding with the 
instantaneous position of the moving axes, are 

Du Dv Dw 

-2M-*x, Di + 2vu-ty. -ft, 

D o o d d 

where -f- = f ^ + u---i-v- \-w . 

Dt c)t ox dy oz 

In small motion all squares and products of u, v, w and their derivatives are 
neglected, and the component accelerations accordingly simplify to 

9 ov ^ . Bw 



The equations of small motion are therefore 

Bu 2 ty I dp 

-- -2a>v-a> 2 x = /- /, 

ct ox p ox 



ow _ty_ldp 

~di ~"dz~~~pdz 9 

where p denotes the pressure and </> is the gravitational potential of the liquid 
mass, no external forces other than gravitation being supposed present. Also, 
uniform density is assumed throughout, so writing 
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these equations take the form 



(1) 



-----, 
8v 



In addition we have for uniform density the equation of continuity 



In any configuration of equilibrium relative to the rotating frame, u,v t w 
would all be zero, and hence 

?X = dx = fy = 
dx dy dz 

Thus x^plp-t- ^(tf + y*) = ^ (3) 



where, in the first instance, ^ is independent of space coordinates. But XQ nmst 
also be independent of time, since in a natural equilibrium configuration p and 
</> depend only on the space coordinates x, y, z. 

Supposing now that there exists some configuration of equilibrium of the 
system, it is proposed to consider the possibility of motions iri which the maxi- 
mum displacement of any particular element of the liquid from its mean position 
is a quantity of the first order of smallness at most, since the object is to arrive 
.at valid linear equations. The resulting non-equilibrium configuration at any 
instant will be referred to as the disturbed configuration, and will differ in 
dimensions from the equilibrium configuration only by small amounts of the 
first order. In the usual manner of treating such motions for the purpose of 
discussing their stability, it is supposed in building up the relevant equations 
that each element of the liquid undergoes an oscillation about a mean position, 
and the reality or otherwise of the period of this motion, when found, deter 
mines whether or not the postulated motion is capable of remaining small. 
Thus, if X Q , y Q) z is the mean position of the clement of liquid that at time t has 
position x, y, z, then we may suppose 



(4) 



where & denotes the real part of the functions following it. (i = ^l.) The 
period 2?r/A connected with the motion is the same for every particle of the 
liquid in accordance with the usual procedure of finding small oscillations. 
(There will be no confusion of the present A with that of the confocal para- 
meters.) The quantities ,77, will be functions of position, and as X-X Q , etc. 
are small, so must be the magnitudes of these functions. It is therefore immaterial 
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to the present order of accuracy whether they are regarded as functions of 
x o> 2/o z o or f x > y> z - As they stand in (4) they will in general involve complex 
coefficients owing to the factor e* A ', and it will be supposed without continuing 
to write the sign i& that the real parts give the appropriate displacements. 

Whatever form these quantities f , 77, take, it is to be expected that as in 
ordinary dynamical examples of small oscillations, they will be undetermined 
to within one or more (infinitesimal) factors that measure the amplitudes, as it 
were, of the oscillations. If this factor (or factors) is reduced to zero, no oscilla- 
tion would take place, since all amplitudes involved are linear multiples of it. 
Thus (# >2/o> z o) must in fact be the position of the element in the equilibrium 
configuration, and accordingly in the motion considered every element oscillates 
harmonically (if A is real) about its position in the assumed equilibrium state. 
Thus x Qt 2/ , 2 refer to the equilibrium configuration. 

For such a motion equations (4) require to the first order that 

_c)x _ 

U = -- lAg< 
vt 



and = - 



with similar expressions for v, ; w and . Inserting these in (1), a small 
motion of this kind requires 



(5') 



8z * ' 

These equations show that if x contains t in any other form than e' A * it cannot 
do so in any terms involving x, ?/, or z. But any terms not containing x, y, z can 
be omitted, since only space differentiations of x enter the equations of motion. 
Thus we may write lU/f ^ ^V ( } (4'\ 



where, again, may have complex coefficients, and the real part of the right- 
hand side gives the x appropriate to the motion. Thus if the oscillation is reduced 
to zero, x must become ^ , and hence /r must reduce to zero. 

Equations (5'), together with the equation of continuity, now take the form 



and 



8y 8z 



(5) 



GENERAL EQUATIONS 117 

The small oscillations of a general rotating liquid are thus seen to depend on 
these four partial differential equations. The functional forms of f , 77, , and 0, 
as is usual in such circumstances, will depend on the boundary conditions obtain- 
ing in the particular problem to which the equations are to be applied. 

As already mentioned, the expressions giving the motion must, as for systems 
of a finite number of degrees of freedom, contain one or more (independent) 
arbitrary infinitesimal constants determining the general amplitude of the 
motion, in terms of which all other constants concerned, for each particular A, 
must be capable of being linearly expressed, and if this one constant is made 
to vanish so must the others. The equations are in fact seen to be satisfied by 
f = 77 = = = 0, and if ,77, , i/j is a solution, so also is k,krj,k,,ki// for any 
small constant k. The equation determining the periods, or A, is to be obtained, 
as for finite dynamical systems, by eliminating all the constants associated with 
the various amplitudes of the oscillation. 

The surface condition for an ellipsoidal configuration 

If the liquid is bounded by a free surface, as in the problem to be considered, 
then p = at this surface whatever its form. From (3) the surface value of x 
is therefore 




2 (X 2 + y 2 ), Free surface of 

. equilibrium configuration 

where now <f> is the gravitational potential of 

the disturbed configuration. To calculate <f> at 
a general point, the disturbed configuration 
may be regarded as consisting of the equili- 
brium configuration together with a layer of 
liquid spread on its surface and of depth 
PN (Fig. 19) given by Fig. 19. 

PN = l(x - x ) + m(y - y Q ) + n(z - z ), 

where Z, m, n are the direction cosines of the outward normal to the equilibrium 
surface. (At places where the original STirface is depressed the density of the 
layer must be considered negative, its whole volume always being precisely 
zero.) Since PP is an infinitesimal of the first order, the contribution to the 
potential (at a general point of space) by the column PN may be calculated 
as if it were all situated at P with a resulting error of second order only, while 
the self-gravitational energy of the layer is also of the second order of smallness. 
The surface layer may therefore be dealt with, to the first order, by regarding 

a mass of amount 

mrj + n) e tA< per unit area 



as added at the point (x ,y Qt z ). Let us denote the potential due to this layer 



While the foregoing simplification is possible in calculating the potential of 
an infinitesimal column surrounding PN, it is seen that a similar simplification 
cannot be made in formulating the condition expressing the vanishing of the 
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pressure at the free surface. For when an oscillation is present, this free surface 
is no longer at (x ,y ,z Q ) but a ^ some point (x,y, z). Hence the contribution to 
the total mechanical potential by the gravitation of the undisturbed configura- 
tion plus the centrifugal force must be evaluated at (x,y,z) in order to obtain 
the value of p correct to first-order accuracy. The gravitational potential of the 
surface layer can, on the other hand, be evaluated at (x , y , z ). 

If XE * s * ne total mechanical potential at a point (# , y Q , z ) of the free sur- 
face in the equilibrium configuration, the potential at (x,y, z) in the disturbed 
configuration will be 



so that x(Xj y, z, t) = XK 

But, as explained above, xts * s constant over the equilibrium free surface (and 
indeed throughout its interior) and equal to x . Hence, since x ~ ' 
we have 



Now -~ E - , ^, ~rT E are simply the components of total force (gravitational plus 

" x o Wo "*<> 
centrifugal) at the surface of the equilibrium configuration and are therefore the 

components of (apparent) surface gravity. Thus if g(x Q , y& z ) denotes the value 
of gravity 

, } _ / . . . i 

?2/ 0) ~ 



and hence t/f(x, y, z) = g(l + mj] + n > )- <f) L . (6) 

Assuming now that the undisturbed free surface is the surface of the ellipsoid 



then 



It will be recalled that p = mabc, where w is the quantity defined on page 77. 
Writing now 



then we have, rewriting (6), \jj = gpcr-<l> L . (8) 

(From this point onwards the suffix can be dropped from X ,y 0f z .) 
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But it was shown (p. 85) that the product g-ar, and hence gp, is constant over 
the surface of the ellipsoid and given by 

gp = 
wherein L r = ^(A + c 2 ), 

a - 

1 ~ 

and A = is of course the ellipsoid under consideration. The factor Gp may 
conveniently be taken as unity. 

To arrive at the surface condition in an appropriate form it is only necessary 
now to express a in terms of surface harmonics. The quantities f , 77, f are func- 
tions of position on the surface of the ellipsoid, and hence so is a. It can there- 
fore in general be developed in a series of Lame* functions (involving of course 
only M ,v) thus 



(9) 

k 

where the ^'s are (small) constants. The potential <f> L due to a surface layer of 
depth pa spread on the ellipsoid may now be written down at once from page 82, 
and is 

M k N k , (10) 



where n is the order of the Lame functions L k> M k , N k . Hence on the surface of 
the ellipsoid we have from (8) 



11 - ^ L k S k M k N k . (II) 

\3 l l 2n+l k */ A=0 

For brevity we write 

and 74 



Then the boundary condition to be satisfied by the functions f, 17, , is seen to 
be equivalent to the following condition: 

If the constants of the ellipsoidal harmonic expansion of the function 



or = 
are denoted by A k , then the expansion constants of i/f on the ellipsoid must be 



Expression of the boundary condition in terms of fy and its derivatives 

Solving the first three of equations (5) for ,77, gives, provided A(A 2 - 4o> 2 ) =j= 0, 



(12) 
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and inserting these in the equation of continuity shows that iff, in space, must 
satisfy the second-order partial differential equation 

-* 



This equation is usually referred to as Poincare's equation. 
We have also, from the definition of cr, that 



_ , lM 

a 2 d* & 2 V \ A 

= D A (0) say. (14) 

It may be remarked here, from the above definition of -D A (0), that if tfj is a poly- 
nomial of degree n so is D x (t/j). The problem therefore reduces to the following 
requirement: 

To find a solution tf/ t of PoincarPs equation, defined everywhere in and on the 
ellipsoid and such that on the surface of the ellipsoid /r and /? A (0) have the constants 
of their ellipsoidal harmonic expansions respectively of the forms 

Zn(JI -II k )A k and 

Properties of Poincare's equation 

Before proceeding to this problem it is necessary to establish certain properties 
of Poincare's equation. These are stated under Theorems I-1V below. 

We shall suppose that A can take any values, real, complex, or pure imaginary, 
but with the exception of real values in the range 2o> ^ A ^ 2co. That is, the 
only real values of A permitted are such that A 2 > 4aj 2 or 1 4o> 2 /A 2 > 0. Arguments 
based on this restriction can be applied directly to discussing stability, since this 
may be dealt with by showing that imaginary values of A are impossible. That 
is, we can assume imaginary values for A and show that this leads to a contra- 
diction. 

Theorem I. Poincare's equation possesses no solution that vanishes on the, surface 
of the ellipsoid, other than the identical one t/r = 0. 

4o 2 
This result would obviously not be generally true for l--y-<0, for then 

x 2 i/ 2 z 2 

~2 + 72 + ~2~~ l would be a solution vanishing on the ellipsoid if A were such 



that 1- ~ =- 

To prove this we make use of Green's Lemma, or the so-called divergence 
theorem, 



in which the integral on the left will be applied to the surface of the ellipsoid 
and the volume integral to its interior. Denoting by tft the conjugate complex 
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function to a solution of Poincare's equation, let us take for 7, V, W the 
expressions - , . 



Inserting these in Green's formula gives 
f M, ty tofi L 4o>V(/r) 10 ff(W ty &A &A /, 4ajV<A ftA) , , , 

Jrr^aM 1 -^]^ 

the remaining terms, multiplying and arising on the right, disappearing since 
j/r satisfies Poincare's equation. 

If now it is assumed that ifj vanishes on the surface, then so also must ^, and 
hence the left-hand side of the above is zero. Thus in that event the integral 

f^i/r tiili ?>ilf Bib 

on the right must also vanish. But , .7,- is real and positive, as also is , .- r . 

c)x fix dy dy 

Hence if A 2 is complex, or pure imaginary, the vanishing of the imaginary part 
of the integral on the right gives 



But r -.~-^() essentially everywhere, and therefore, in fact, -, = 0. Hence we 

cz dz cz 

toti ?d 

must also have J J = 0. 
dx cy 

In the same way, if A 2 is real and 1 4w 2 /A 2 > 0, the same result holds, since 
all three terms on the right are now essentially non-negative. Hence // is constant 
within the ellipsoid, and since it vanishes on the surface it must be zero everywhere. 

It follows immediately that Poincare's equation can have only one solution 
taking given values on the surface of the ellipsoid. For if fa and 2 were two 
essentially different solutions taking the same given values at the surface, then 
ifj l i/J 2 would be a solution taking zero value on the surface. Accordingly fa 1// 2 
would be zero everywhere within the ellipsoid, that is t/j^i/j^. 

Theorem II. There always exists a polynomial P(x,y,z) of degree n that satisfies 
Poincare's equation and takes on the ellipsoid the same values as an arbitrarily 
given polynomial Q(x, ?/, z) of the same degree. 

This will certainly be true if it can be shown that there exists a polynomial 
R(x, y, z) of degree n 2 such that the polynomial defined by 



is a solution of the equation, i.e. Vf P 0. (The polynomials involved here are 
not being supposed homogeneous.) Now R, being of degree n 2, contains 
^n(n 2 1) coefficients which are as yet at our disposal. Writing down the condi- 
tion that Vf P = gives 



( r 2 ,,,2 2 2 \ fift 

- 2 + ; 2 + 2 -l VfU + 4 (const.)*?, + (const.) B = 0. 
a u U C J ? t V,o v% 
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The expression on the left-hand side of this equation also has $n(n 2 1) coefficients, 
since it also is of degree n 2. Also all the coefficients of R occur only linearly, 
and hence the conditions for it to vanish identically provide just enough (linear) 
equations to enable the $n(n 2 1) coefficients of R to be uniquely found in general. 
Now in order that the process of solution shall be possible, it is necessary that 
the determinant of the equations for the unknown coefficients in R (that is, of 
the left-hand sides of the resulting equations for the coefficients of R) shall not 
be zero. It can easily be shown that this requirement is fulfilled. For if we 
make the equations homogeneous by putting Q = temporarily, the resulting 
determinant cannot vanish, for if it did, it would mean that each of the coeffici- 
ents in R could be expressed as a multiple of a selected one of them. That is 
\n(n 2 1) 1 of them could be expressed as multiples of the remaining one, 
simply by omitting one of the equations and solving the remainder. Denoting 
the selected coefficient by a r this would then give 

R = a r (polynomial of degree n-2 not identically zero). 

(.2 yZ ^2 \ 

~2+/2+ a ~~ M ^ wou ld be a polynomial satisfying Vf = and Vanish- 
tt U C ] 

ing on the ellipsoid. Since by Theorem I this is not possible, it follows that the 
determinant of the coefficients cannot vanish. Hence given Q, an R can always 
be found, and so P can be found. 

Definition. If Q is chosen to be a Lame* polynomial of order n, viz. of the form 




where V 2 Q = 0, the corresponding polynomial P will be called a Poincare poly- 
nomial. 

Since, in space, Q will be of the form L k M k N k , it follows that both polynomials 
will reduce on the surface of the ellipsoid to a Lame surface function M k N k . 



Theorem III. If P is a Poincare polynomial of degree n, the function 
reduces on the ellipsoid to a sum of Lame surface functions M k N k of order n only. 

To prove this, let it be supposed that $ and i// are any two different solutions 
of Poincare J s equation. Then we have, p denoting as usual the perpendicular 
from the centre to the tangent plane at an arbitrary point of the ellipsoid, 



c*~8z A a 2 By b 2 dx 



and similarly, changing A to A, 
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where Z, m, n are the direction cosines of the outward normal to the ellipsoid. We 
then have that the surface integral over the ellipsoid 



,8>fi 



, 8 <t> 2io>.8>/> 2io> ,8<f>\ 

arir^--T*^ 

8J> 2io> . 8<li 2ia> . 8<k\ /, 4o> 2 \ 8 / . 8Ji , 8d>\\ , . , 
+ ^^^ 



since <f> and satisfy Poincare*'s equation. 

Let it be supposed now that P is a Poincare polynomial of degree n, and 
reducing therefore on the ellipsoid to a product of Lame functions MN of order 
w, and let us suppose $ P. Also let us suppose <f> = P' where P' is any other 
Poincare polynomial of degree m<n. (The w, n here are of course integers and not 
the conventional symbols for direction cosines in the surface integrals.) Then by 
the preceding result 



But the right-hand side of this vanishes, for 7)_ A (P') is simply a polynomial of 
degree m<n, and therefore expressible on the ellipsoid as a sum of Lame func- 
tions MN of orders m and less, and therefore of orders all less than n. That is 

7)_ A (P') = ^A k M k N k the order of every term being less than n. 
Also P = MN wherein the order is n. 

But by the orthogonality property 

k N k pdS = to. 



[[ 



Hence p'D^pdS = 0. 



That is, D^(P) on the ellipsoid is orthogonal to all Lame functions of order less 
than n. But in space D X (P) is a polynomial of degree n, and hence because of 
this result it must reduce on the ellipsoid to a sum of Lame functions of order n 
only. (If any other order were involved the integral could not vanish for every 
P' of order less than n.) 
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If therefore the 2n + 1 independent Lame polynomials L M N of order n are 
considered, and P l9 P 2 ,.--> ^n+i denote the corresponding Poincare* polynomials, 
then on the ellipsoid there must exist relations of the form 



wherein the constants ^ fc depend rationally on A since D A (P) contains A only 
rationally. 

Similarly D_ A (P) reduces on the ellipsoid to a sum of Lame functions of degree 
n and the coefficients are again rational in A. 

Theorem IV. If iff is any function satisfying Poincar&s equation, the expansion 
constants of order n of U\($) considered as a function of position on the surface of 
the ellipsoid, are linear combinations withjijced coefficients of the expansion constants 
of order n only of 0. 

To establish this we make use of the result arrived at in proving Theorem III, 
namely, .... .... 

UDM)pd8 = W 

/ / .' V 

In this let us write c/> = P, a Poincare polynomial of order n, and therefore reduoing 
to a Lame product MN of order n on the ellipsoid. In space P will correspond 
to a certain Lame polynomial LMN also of order n. The corresponding expansion 
constant in D^jt) will be 



and this is equal to \i/j/)_i(P)pd8. 

But, by Theorem III, D_ A (P) reduces on the ellipsoid to a sum of Lame functions 
of order n only. Hence this integral is equal to a linear sum of the expansion 
constants of of order n only, which establishes the result. 

From this it follows that if in the development of 0, considered on the surface, 
into a sum of Lame products lf\, say, = MN, and we consider tho terms of order 



then to each F i will correspond a Poincare' polynomial P, and we can consider 
i/r in space to be given by 



Then by Theorem IV 

2n + l 

# A (<A) = 2 Zc,, |fc p * in s P ace 



on *h e ellipsoid. 
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Hence the complete value of Z> A (^) on the surface will be 



Now the problem to be solved is to find a solution ifj of Poincare's equation 
such that, on the ellipsoid, and D A (0) have the constants of their ellipsoidal 
harmonic expansions respectively of the forms 2Tr(H II k )A k and (A 2 4oj 2 ) A k . 
In order to satisfy this requirement we thus obtain the following conditions on A 



2n X h ik (H -H i )A i (fc= l,2,...,2n+l), (15) 

i=i ' 

and these conditions evidently constitute 2n + 1 homogeneous linear relations 
between the A k in which the coefficients h { k depend rationally on A. 

The expansion constants of o-, i.e. in effect of Z> A (0), are therefore given by an 
infinity of such linear and homogeneous equations. But this system breaks up 
into an infinity of partial systems each consisting of a finite number, 2n + 1 , in 
this same number of unknowns. That is, the constants corresponding to a given 
order n are involved in *2n + 1 homogeneous equations that contain none of the 
constants corresponding to any other order. The equations for the permissible 
values of A are obtained by equating to zero the determinants of each of these 
partial systems, and accordingly lead for each separate order n to a rational 
algebraic equation for A of a finite degree. 

It is shown later that the degree of the equation in A corresponding to order n 
is always 7i 2 + 4?&+ 1. The roots of this equation then give the frequencies, A/277, 
of the possible (small) free oscillations, and then for any selected root, A, equations 
(15) give the corresponding ratios of coefficients A l : A 2 : ... : A 2n , x associated with 
the harmonics in the oscillation. It has been further shown by Cartan that the 
procedure explained above in fact applies for all frequencies including also any 
within the range 2co to 2w. 

Ordinary stability 

It will not be necessary to our purpose to obtain in detail expressions for the 
small oscillations of the Jacobi ellipsoids when these arc secularly stable and 
therefore also ordinarily stable, but to consider only the question of whether 
they remain ordinarily stable beyond the configuration at which secular stability 
first ceases. It is assumed throughout, as its definition, that ordinary stability 
means that all the roots in A must be real, for if there were an imaginary or 
complex root there would be a motion in which the displacement (measured to 
the first order) increases indefinitely and therefore, at least, could not be regarded 
as remaining small. Moreover, it is not possible in the present case for a motion 
depending on a term such as e~ kt to occur without there also being a term e kt . 
For to each solution of equations (5) there exists a corresponding solution sym- 
metrical with respect to the #z-plane but with the sign of A changed, since the 
equations are unaltered if the signs of ?/, 77 and A are changed. This shows that 
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the roots of the final equation in A derived from (5) are equal and opposite in 
pairs, as might have been expected but cannot be asserted from the general 
result for finite rotating systems. 

If the ellipsoid is secularly stable for displacements of order n, the potential 
energy will be an absolute minimum and therefore the system will automatically 
be ordinarily stable, and so the equation in A will necessarily have all its roots 
real. Accordingly, as the system evolves along the Jacobian series in the direc- 
tion of increasing angular momentum, the question of ordinary instability first 
arises for third-order harmonic deformations. It will next be shown that ordinary 
stability ceases at the same time as the secular stability for these displacements. 
From a physical standpoint, this property amounts to the Jacobian ellipsoids 
being ordinarily unstable at all stages beyond the critical form at which the 
pear-shaped series bifurcates. 

To establish this it is necessary to consider in further detail the actual equation 
for the possible values of A, but it turns out that the result depends only the 
highest and lowest terms of the equation, and we proceed to consider these. 

The degree of the equation for X for oscillations of order n 

It has been seen that the small oscillations are to be obtained by integrating 
the equations 



dx 



ex 



In these, is a polynomial of degree n in x,y,z\ and f, 77, are polynomials of 
degree n\ in x, y, z. These polynomials are not homogeneous, but the number 
of coefficients in the terms of highest degree in /r is \(n -f 1) (n + 2), and in each of 
,17, is ln(n+\). 
Now it has been seen that \fj and a must be of the related forms 



and a = f + ij + p t = *A k M k N k . 

Hence if R k = L k M k N k is any Lame polynomial of order n, we have 



and 



THE DEGREE OF THE EQUATION FOR A 127 

Hence, combining these results, we have the following 2n + 1 conditions on 
0> f > *? arising from the relationship between if/ and a 

+ ,, + { pdS, (17) 



all the integrals being taken over the surface of the ellipsoid. 

As has been seen, the equation for A, giving the frequencies of the oscillations, 
is to be obtained by eliminating the coefficients of the polynomials 0, , 17, . 
To obtain the frequencies of order n, it is necessary to effect a process of elimina- 
tion on the terms of highest degree in x,y, z in the polynomials ^,,17, , for it 
is these that contain (only) the constants associated with the MN's of order n. 
The lower terms will contain constants associated with surface harmonics of 
lower orders. It will be convenient, however, to retain the same letters (i/r, f , 17, ) 
to denote the relevant parts of these polynomials comprising only their highest 
order terms. 

To consider the elimination in further detail, let us write = A(, then equa- 
tions (5) or (16) become 



& - " s - 

together with 

(19) 



Suppose the elimination is begun by dealing first with the coefficients of the 
polynomial <, which refers now to simply the highest terms of the original <f> = ifj/X. 
If we consider the three equations (18a) involving <f>, the coefficient of a term 
like x r , or if, or z r in <f> will appear once only; the coefficient of a term like x r y s , 
or y r z 8 , or z r x s will appear in two places; while the coefficient of a term like 
x r y*sf will appear in three places. For example, suppose the polynomials are 
denoted by: 
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r\ i 

The equation - - = Af + 2ia>rj must then give rise to the identity 

GX 

^ 



Equating coefficients of the various terms gives the following necessary relations 
between the original coefficients: 

x n ~ l : na n = Xb n _- l + 2iajc n _- L " 

r-i s _ o- > (20) 

The coefficients in </> associated with terms such as x r that appear once only in 
these equations obviously cannot be eliminated at this stage. On the other 
hand the coefficients in (f> that qualify terms like x r y" appear twice and therefore 
can be eliminated. The number of such coefficients is 3(?i 1). Since the <-terms 
are not multiplied by A in the three equations (18a), the elimination of these 
coefficients leads to 3(n 1) relations linear in A and in the coefficients of the 
polynomials , 77, . 

Considering next the coefficients in <f> that qualify terms like x r y ft z i , these appear 
three times each in equations (18a), and the total number of such coefficients is 
\(n l)(n 2) [viz. 7\(?i-f 1) (n + 2) 3(?& 1) 3]. Each of these can be eliminated 
in two independent ways, arid provide two equations linear in A and the coefficients 
associated with ,77, ; that is, (n 1) (n 2) equations in all. 

Thus the elimination in this way of those coefficients of </) that can be eliminated 
by equations (I8a) leads in all to 



equations, in each of which A and polynomial coefficients occur only linearly. 

It will be noticed also that none of the coefficients of <f> appears multiplied 
by A in any of the equations derived from (I8a), and hence if the values of these 
coefficients in terms of those of , rj, (and A) are substituted on the left-hand 
side of (19), then A will occur both as A 2 and A. 

The total number of coefficients remaining now to be eliminated is simply the 
number of coefficients of f, rj and taken together, that is |n(n+l). For this 
purpose there are already available n 2 1 equations linear in the coefficients of 
f , 77, f and containing A to the first power. There are also the 2^+1 equations 
provided by (19) which contain A both to the first and second powers. And 
finally there are the equations arising from div(,7/, ) = which give \n(nl) 
relations not containing A at all. The total number available is thus 
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which is precisely the requisite number. The elimination of the coefficients of 
, 17, can therefore be effected by means of a determinant A M (A) in which 

the first n 2 1 rows contain elements of order A, 

the next 2n+l rows contain elements of order A 2 , 

the remaining \n(n 1) rows contain elements independent of A. 

Hence the degree in A of the eliminant is at once 



Since the roots in A of A H (A) = are equal and opposite in pairs, this number 
7fc 2 + 47i+ 1 when even gives twice the number of independent harmonic motions 
possible, and when odd is twice this number plus one. Thus for n = 3, the total 
number of independent harmonic vibrations of a stable Jacobi ellipsoid is eleven. 

The degree /t 2 +4/i+l is always attained 

The next step in the argument is to show that this degree is always reached. 
Tf in any circumstances it were not attained, this would mean that if in the set 
of equations finally giving rise to the determinant A,, (A) only the terms of highest 
degree in A in each element were retained, the resulting set of equations would 
be linearly dependent. 

It is to be remembered that those coefficients of A that can be eliminated are 
first removed to provide 3(n 1) relations between the coefficients of the , 77, f 
polynomials, and these are some of the equations in which now only terms in A 
are being retained. Writing for the moment A^ = <, then first of all the coefficients 
of x are eliminated. Thus, the linear dependence, after retaining only the highest 
terms in A, is equivalent to asserting the consistency of the following set of 
equations 



and R k xpdU = 0. 



~dx ?> 'By r]i dz *" 

2+5+2-* '[ 
J> 

But if this set of equations held, it would mean that x is a harmonic polynomial 
of degree n orthogonal to all the Lame polynomials of degree n, and this cannot 
be so. 

Hence it follows that the degree of the highest power of A occurring in A M (A) 
is essentially n z -f 4/i + 1. 

Form of the constant term in An(X) 

This term will be obtained simply by putting A = in every element of the 
determinant. Now if in the 2n+ I rows that arise from (17) we put A = 0, the 
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remaining terms have a common factor of H H k . For putting A = 0, (17) 
becomes 





and the right-hand side is linear in the coefficients of , 77, and also independent 
of A. Hence in A W (A) the constant term contains each of the coefficients of stability 
H H k as factor, and thus 2n+ 1 such factors in all. We must therefore have in 
the first place 2ri+1 



where M is some multiplying factor, to be considered next. This form is what 
would be expected by analogy with systems of a finite number of degrees of 
freedom (see p. 20). 

We proceed to consider now whether M vanishes in any particular circum- 
stances. If M were zero, it would imply the consistency of the equations got 
by putting A = in (16) and (17). That is, the following set of equations would 
be self-consistent: 



8, 



Assuming this to be so, it is simple to find what the equations imply. For, since 

r\ I 

-- = 0, must be independent of z. Hence f , 77 are also independent of z, since 

dz 

they are proportional to derivatives of </>. Also we then have 

o 



- _ 

^ fa By 2ii*> \dxdy 
Hence is independent of z. It follows then that the expression 

*^ 

a*dy 

is a polynomial of degree n containing z only to the first power, where it occurs 
explicitly in the last term. But then the last of equations (22) means that this 
polynomial (linear in z) reduces on the surface of the ellipsoid to a sum of Lame* 

y ij z 

polynomials of order less than n. But the degree o^-^f + ^^ + ^S is necessarily 

a o c 

n, since ,77, are each of degree n-l. Hence the situation must be that 
5 + f^ + Ji= (Polynomial of degree w-2) g + g + ?* - 

+ polynomial of degree n-l or less. 
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But f , rj, now represent homogeneous polynomials (since we are retaining only 

terms of order n in finding A n (A)), and hence ^f + y4>? + 4is also homogeneous 

doc 

of degree n. If therefore we rewrite the above in the form 



(y^ ifi z 2 \ 
(homogeneous polynomial of degree ^-2)1-2 + + 2) 

+ polynomial of degree n-l or less, 



this must be homogeneous of degree n, and hence the last term, in italics, 

X tJ Z X 2 I/ 2 Z 2 

must be zero. Thus - O f+ri ?? + - must be divisible by + + - which is 
a 2 6 2 ' c 2 a 2 b 2 c 2 

impossible (since it contains z only to the first power) unless it is identically 0. 
This requires first = 0, since the term in z must vanish identically, and second 

x d<f> _ y d</> 
~ b*~dx 



x 2 u 2 
This second requirement means that ^ depends solely on ~2 + ^ 2 , and since it is 

homogeneous it must be proportional simply to some power of this; moreover, 
since it is of degree n, we must have 



and this is only a rational polynomial if n is even. 

It may be concluded, therefore, that the equations (22) are self-consistent and 
that M vanishes when n is even. 

While if n is odd, M is essentially non-zero. 

It would be expected that M = for n even, since then the total number of 
roots n 2 + 4n + 1 is odd, and since the roots are equal and opposite in pairs A = 
must be a root. The case of special importance to our purpose is n = 3, and the 
constant term is then 

2n+l 

A n (0) = M (non-zero) H (# -/4). 
i 

(Note: When in the equations resulting from (18) and (19) we retain only the 
terms in A, or later only the terms not containing A, the function <f> then appearing 
denotes of course only that part of the original <f> balancing the retained terms 
on the right.) 

Ordinary stability ceases simultaneously with secular stability 

We can now establish the important result that the Jacobi ellipsoid becomes 
ordinarily unstable simultaneously with its becoming secularly unstable through 
the changing of sign of the characteristic coefficient of stability of order three. 
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For if we denote C n , the characteristic coefficient of stability of order n, by 
H H 2n+lt then for (7 rt >0, while the ellipsoid is thoroughly stable, we have by 
the foregoing result 

(23) 



where we may suppose A x < A 2 < ... < A m are the positive roots of A W (A) = 0, and 
are m in number where 2m = n 2 + 4n+ l. Also 6 is the product of M with the 
remaining non-vanishing coefficients of stability, and must therefore always 
remain positive since it is positive initially. On the other hand, for a Jacobi 
figure just beyond the critical ellipsoid, we must have 

A?At...A? tt =0(# -// 2n+1 )<0. (24) 

Now, by continuity, only the smallest of the roots, viz. Af, can have changed 
sign, and hence now X\ must be negative, A t therefore imaginary, and the system 
ordinarily unstable. 

It is not necessary to consider specially the possibility of the coincidence of the 
smallest roots of Af = A|, for then we should have simply AJ negative, and so A x 
would be complex, and exponential terms with real indices would enter just the 
same. 

Finally, it may be noted (cf. p. S4) that for a third-order harmonic deformation 
the angular velocity is not changed to the first order. Hence the instability Is a 
true one not capable of being removed by selecting a frame of reference with a 
slightly different rate of rotation. 

This completes the proof of the ordinary instability of the Jacobi ellipsoids 
beyond the point at which the pear-shaped series bifurcates. 



Chapter X 
COSMOGONICAL IMPLICATIONS 

Impossibility of the fission hypothesis 

According to the fission theory, the origin of binary stars is regarded as due to 
the break-up of a single mass by its rotation. There are a number of grave 
astrophysical difficulties in this theory such as the question of the source of the 
necessary angular momentum, whether the density distribution in stars and its 
evolution are consistent with fission al disruption, whether the distribution of 
mass-ratios and separations are explicable on the theory, and the whole general 
problem of stellar evolution itself. We will not attempt to discuss these here 
but will concern ourselves solely with examining the fission process itself, as 
proposed by its advocates, on purely dynamical grounds. A general secular 
increase of angular momentum will be assumed, with the density remaining 
uniform, since this is equivalent to constant angular momentum with density 
increasing gradually. These are the assumptions on which both Darwin and 
Jeans based their ideas on fission. 

Now, as already explained in the introduction, had Darwin's conclusion, that 
the pear-shaped figure is secularly stable, been correct, it might then have been 
fairly plausible to suppose that the deepening of the furrow with evolution along 
the series was some indication that the mass would eventually divide into two 
parts in orbital motion about each other (though information on the initial 
stages of this process by no means necessarily secured that the pear-shaped 
series itself did not bifurcate later to some new form). But when Jeans's studies, 
in agreement with those of LiapounofF, contradicted Darwin's conclusion, the 
sole grounds on which Darwin based his description of the fissional process were 
completely removed. Yet Jeans nevertheless finally maintained exactly the same 
outcome of the process as Darwin, namely, fission into two detached masses mov- 
ing in almost circular orbits about each other. Jeans's view obviously amounted 
to asserting that the evolution of the mass would be quite independent of whether 
the pear-shaped series happened to be secularly stable or unstable, and hence that 
the whole of the investigations establishing the incorrectness of Darwin's con- 
clusions, or indeed any studies of this problem at all, were valueless from the 
point of view of cosmogony, since the outcome of the process in either case 
would be the same. As the fission theory will always be associated with Jeans, 
and as its acceptance is still one of the main errors besetting even modern writ- 
ings on stellar evolution, it seems worth while dealing in some detail with the 
many points at which Jeans fell into error. 

The first of these concerns the nature of the difference between secular and 
ordinary stability of which at many parts of his work Jeans seems to have had 
an inadequate conception. Thus at one point* he says: 

* Astronomy and Cosmogony, p. 199. 
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... as the physical conditions of a system gradually change secular stability neces- 
sarily sets in before ordinary stability. Thus, for problems of cosmogony it is secular 
instability alone which is of interest. A system never attains to a configuration in 
which ordinary instability comes into operation since secular instability must always 
have previously intervened. 

Now this statement is plainly incorrect, since, although a system cannot become 
ordinarily unstable before becoming secularly unstable, the two kinds of stability 
can vanish simultaneously. Any finite system for which a single coefficient of 
stability changes sign, for example, will become both ordinarily and secularly 
unstable at the same stage, since the product of the squares of the frequencies 
of oscillations has always the same sign as the product of the coefficients of 
stability. A simple instance of this was discussed in Chapter II, where it was 
seen that a particle in a circular orbit under a central force \i\r n becomes simul- 
taneously both secularly and ordinarily unstable when n rises above 3. Also we 
have seen that the same actually holds for the Jacobian ellipsoids, which become 
both secularly and ordinarily unstable for third-order harmonic displacements at 
the same stage. Naturally enough Jeans, holding the view he stated, failed to 
see the possible importance of ordinary instability. But as it happens, it can be 
shown on general grounds, from the property that the ellipsoids become ordinarily 
unstable, that the fission process as described by Jeans is dynamically impossible. 
For the essence of ordinary instability is that it is independent of friction, and 
therefore that any motion ensuing from it must be of a strictly reversible char- 
acter; that is, if t is written for t in the equations of motion they must remain 
unaltered. But if in the complete absence of friction we imagine a close double- 
system as the product of rotational instability of a single ellipsoidal mass, a 
reversal of the motion of the bodies simply leads to them describing their orbits 
in the reverse direction and not to their re-uniting into the original mass. Thus 
at a single blow the fact that the system becomes ordinarily unstable disposes 
of the fission hypothesis for the formation of double systems. 

But Jeans seems even to have misunderstood the nature of secular instability 
too. It is clear, from the fact that the initial members of the pear-shaped series 
have less angular momentum than the last stable ellipsoid, that the system cannot 
evolve along the pear-shaped series, since this would necessitate a continually 
decreasing angular momentum for the body. Yet Jeans supposed it to be possible 
to find the initial details of the process of break-up by fission by calculating to a 
fairly high order of accuracy the shapes of the members of the pear-shaped 
series in the corresponding two-dimensional cylindrical problem which has 
analogous properties to the three-dimensional problem but is far simpler of 
treatment. The resulting figures (Problems of Cosmogony, p. 116; Astronomy and 
Cosmogony, pp. 220, 221) in fact can establish nothing of the precise course of 
development when the critical Jacobi form is passed. For the ordinary instability 
at this stage means that the precise further development of the system depends 
on the detailed initial conditions of disturbance, and these from their nature are 
essentially unknowable and different from one system to another, though it may 
nevertheless happen that the general outcome of the instability is independent 
of the exact nature of the disturbance. Thus, a rod standing upright on a fixed 
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pivot at its lower end and acted on solely by gravity is ordinarily unstable, but 
the plane in which it falls away from this position and proceeds to move like a 
pendulum depends entirely on the initial disturbance given to it, however small 
this may be. Again, a particle describing an unstable circular orbit under a 
central force, when disturbed may either eventually fall into the centre of force 
or depart an infinite distance from it. There are only the two possible general 
results, but the paths by which the particle may reach either one of them are 
infinite in number and depend on the precise conditions of disturbance. Thus the 
final outcome of the instability of the Jacobi form could conceivably be fairly 
definite though the intervening motion might be subject to complete uncertainty. 
In some parts of his writings on the fission problem, Jeans seems waveringly 
to have held various opinions quite inconsistent with others expressed elsewhere 
in his work. Thus he says:* 

The fact that the pear-shaped series is initially unstable shews that a rotating mass 
cannot evolve by slow secular changes through a series of pear-shaped figures. This 
somewhat diminishes the interest of the pear-shaped series in the problem of cosmogony, 
but nevertheless it remains important to obtain as clear an idea as we can of the nature 
of the series. 

Yet a few pages later,t after having calculated the two-dimensional forms, as 
explained above, he says: 

Thus we may with fair confidence assert that the two-dimensional [pear-shaped] 
series ends by fission into two detached masses and ... it seems highly probable that 
the three-dimensional series also will end by a similar fission. 

At a later point in the same workj Jeans appears to have still another view of 
the implications of secular instability, for he there says: 

The pear-shaped figure is unstable, so that as soon as it is formed dynamical motion 
ensues arid fission results. The masses are at first projected away from one another 
with considerable velocity, but seem likely to settle down finally to describe steady 
orbits about one another. 

In his researches on this problem, Darwin supposed that the initial orbits would be 
strictly circular, but this was because ho believed the process of fission to be a statical 
process and not a dynamical process, as we have seen it to be. According to Darwin's 
view, the changes in the star while fission was taking place wore, initially at least, of 
a purely secular nature, and it was natural to suppose that the final result would be 
two masses rotating in actual contact and at rest relative to one another. 

There being no longer any theoretical justification for supposing that the initial 
orbits will be strictly circular, we have to consider the possibility of the masses being 
thrown apart with appreciable radial velocities, and describing elliptic orbits about 
one another. 

Consider for simplicity the case in which the original star is supposed to divide into 
equal masses, and suppose that fission occurs when the centre of each mans is at a 
distance r from the common centre of gravity. Let each star be supposed to have a 
radial velocity v in addition to the tangential velocity cor in space resulting from 
rotation. Each mass will describe approximately an elliptic orbit in space so that 
after the orbits are partially described the masses will again each be at a distance 
r from their common centre of gravity, but are now approaching each other with a 
radial velocity v. A collision of some kind must occur, and since the masses will not 

* Problems of Cosmogony, p. 102. 
f Ibid., p. 115. 
} Ibid., p. 252. 
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be perfectly elastic, their velocity of recession after collision will be some velocity v 
less than v, while the [transverse]* velocity cor must, from the conservation of angular 
momentum, be the same as before. It follows that the new orbit will bo of less eccentricity 
than the old, and the eccentricity will further diminish at each subsequent collision. 
We cannot argue that the eccentricity will be finally reduced to xero; a limiting value 
will be reached such that the masses just graze one another at periastron. 

This account of the matter shows some appreciation of the likely consequences 
of instability apart from overlooking the possibility that v might be great enough for 
escape to occur, but even so there are a number of other strong objections to the 
suggested course of development. To begin with, the assumption of equal masses 
after break-up can easily be shown to be invalid since such a system would 
necessarily require more angular momentum than the last stable Jacobi form, 
but is of small importance compared with the ideas expressed concerning the 
collisions. If the masses are projected away from each other with considerable 
velocity, as may reasonably be assumed, they will perforce come together with 
practically the same velocity, as stated, but the two masses would re-unite 
colliding liquids do not retain their identities arid separate like billiard balls 
and there seems no possibility of them retaining their individual identity in the 
way Jeans supposes. But even if this difficulty is passed over, there is the objec- 
tion that the angular momentum of the critical Jacobi form is appreciably 
smaller than that of two equal or comparable masses in stable orbital motion 
about each other (see Table VI, p. 142), and the discrepancy would be increased 
by the introduction of any eccentricity for the purpose of keeping the smaller 
mass outside the limit within which it is tidally unstable. 

Although these ideas of Jeans are somewhat confused it may nevertheless be 
the case that when the Jacobi form becomes ordinarily unstable a slight disturb- 
ance will cause it eventually to divide into two masses. For there is little doubt 
that some sort of break-up must occur for the reason that the single system 
possesses more angular momentum than it can store as a stable mass in any 
known form. By some means the mass must rid itself of at least part of the 
angular momentum, and its transference to that of orbital motion of separate 
pieces is obviously a possible way of effecting this. No other way out of the 
difficulty has ever been proposed apart from ejection of a part of the mass and 
no other way seems physically possible. But even so it remains necessary to 
avoid the difficulty arrived at above concerning the collision and re-uniting of 
the pieces. If, for simplicity of discussion, we assume break-up into two masses, 
the only way of avoiding a subsequent coalescence of the pieces is if the initial 
rate of separation is sufficiently great to carry the pieces to infinity, that is, for 
them to have hyperbolic relative speed of separation. If this did not happen at 
first, the collision and re-uniting of the masses would necessarily involve dissipa- 
tion of energy and a rendering of the system more unstable still. (The ultimate 
effect of dissipation owing to a collision would be an increase of density and this 
would be equivalent to an increase of angular momentum without change of 
density.) A further break-up, but this time with greater violence would follow, 

* Jeans has 'radial*. 
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and so on, until a disruption was reached that succeeded in sending the pieces 
apart with hyperbolic speed. It appears that in no other manner can the system 
find its way to a new state of stable steady motion. Each piece will then become 
a separate stable rotating system and the relative velocity will eventually decrease 
to a constant value. 

It does not seem to be generally known, and it is certainly not made clear 
in either of Jeans's books (Problems of Cosmogony, 191 9; Astronomy and Cosmogony, 
1929), that the conclusion originally arrived at by him concerning the stability 
of the two-dimensional pear-shaped figures was incorrect. Thus in 1902* Jeans 
announced that his analysis established that the cylindrical pear-shaped series 
was secularly stable, a conclusion that seemed to lend strong support to Darwin's 
corresponding conclusion for the three-dimensional series. It was then natural 
to suppose that further evolution would take place, at least to begin with, along 
the pear-shaped series with a gradual deepening of the furrow, suggesting eventual 
division into two masses, and in the two-dimensional case the forms could be 
calculated to a fairly high order of accuracy. But many years later in 191 Of 
Jeans briefly mentions that his former value for the moment of inertia of the 
two-dimensional pear was in fact wrong owing to a simple numerical error and 
that the corrected value implied instability. This signal discovery appears to 
pass unmentioned in both the 1919 and 1929 books, though the unstable nature 
of the forms is concluded there, but curiously enough Jeans retains the remaining 
conclusions of his 1902 paper in detail and again claims, as then, that the mass 
will still evolve along the unstable series. 

To sum up on the question of the fission theory of the origin of binary stars: 
even if all its extremely doubtful initial assumptions are permitted to give the 
theory its most favourable conditions for success, we are led purely by dynamical 
evidence to the conclusion that the process of instability cannot result in the 
formation of a double system consisting of two comparable masses moving in 
close orbital motion about each other. When it is remembered that added to this 
the particular theory of stellar evolution of which fission forms a crucial part itself 
depends to a large extent on conjectural speculations on the internal structure 
and source of energy of the stars, long since known to be invalid, there appears 
to be no remaining evidence in favour of the fission theory. 

The origin of satellites within the solar system 

It does not of course follow that the process of break-up with ejection of part 
of the mass to infinity has played no part at all in cosmogonical changes in any 
celestial systems, but if the present dynamical theory is to be applicable to even 
a moderately close degree, it can only be for systems having some fair approxi- 
mation to uniform density in the greater portion of their mass. There is reason 
to believe that this holds for the main bodies of the planets of the solar system, 
though for reasons that will become clear we shall not attempt to apply it directly 
to any of the planets in their present forms but instead only to what may be 

* Phil. Trans. 200 A, p. 96. 
t Ibid. 217 A, p. 28. 
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termed primitive planets, by which is meant the forms that these bodies may 
have taken at an earlier stage in their development. There are cogent physical 
arguments for supposing that in its early stages the solar system developed 
through a state that consisted of a number of bodies having planetary distances 
from the sun and arriving at a condition of rotational instability either through 
increase of density or through the addition at their surfaces of matter endowed 
with the necessary angular momentum, or through possibly a combination of 
both processes. 

According to certain hypotheses of the origin of the solar system, the planetary 
material was initially detached from that of a companion star to the sun moving 
at such a separation for its material to have angular momentum per unit mass 
comparable with that at present associated with the great outer planets. It is 
not necessary for our present purpose to make any decision between the various 
hypotheses of this kind, for in all of them the result is arrived at that only a 
few large planets, comparable in mass at least with those of the present great 
planets, could be formed initially. The problem then remains how the numerous 
satellite systems have arisen. The satellites are of such small masses that there 
is no possibility of them forming directly by condensation from material at 
stellar temperatures, nor could they form and remain in orbital motion round 
a planet itself forming by condensation and contraction from a disc of material 
revolving round the sun. Yet another difficulty is the existence of the four much 
smaller terrestrial planets which exhibit surprising differences of composition 
among themselves. Of these planets Mercury in particular, and Mars to some 
extent, are more nearly comparable in mass with the largest satellites than with 
other planets. 

The problem can be approached from a slightly different point of view by 
considering the general question of the contraction through gravitational causes 
of a mass of material endowed with angular momentum. With continually 
increasing density a state of rotational instability must always arise sooner or 
later unless some change of physical state enters the system to arrest further 
increase of density. In most stars, for instance, energy generation sets in in the 
central regions long before the density increase becomes consistent with rotational 
instability of any kind. In a planet what evidently finally protects the system 
from rotational instability is its reaching a liquid or solid state such that further 
contraction ceases. But there is no reason to expect that the angular momentum 
per unit mass associated with the material from which any primitive planet 
formed would be insufficient to cause instability before such a state was reached. 
If the rotational momentum were great enough, a primitive planet would become 
unstable before this and result in a fissional break-up with the ejection to infinity 
of a part of the mass. For uniform density it can be shown in the following way 
that the portion so lost to the main body must be less than about one-eighth of 
the original mass. 

Suppose the original mass M breaks into two separate pieces of individual 
masses mM and (1 m)M, where m<%, without change of density, and which 
eventually just escape from each other to infinity. The speeds of rotation of these 
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masses when they are first formed would probably not differ much from that of 
the critical Jacobi ellipsoid, but their shapes will not be critical ellipsoids. They 
could, in fact, possess this particular angular velocity (given by w*l%TrQp = 0-142) 
while having the form of a thoroughly stable Maclaurin spheroid of eccentricity 
of section about 0' 71, as is shown by Table I, p. 39. The masses could, however, 
retain most energy when in the critical Jacobi form, and if we attempt to impose 
this upper limit of internal energy on the dividing pieces, the condition they can 
just succeed in escaping from each other is, at once 

- 0-745 GM^pV* > - 0-745 G(mM)p - 0-745 G(l - m) 5 / 3 Jlf 5 /V /3 , 
or 



and this cannot be satisfied for any non-zero ra in the range 0<ra<J. This 
establishes, therefore, that if the pieces into which the original mass breaks 
eventually separate to infinity, they will do so as stable bodies (excluding, of 
course, the possibility of further internal contraction). 

If at the other extreme it is assumed that the two resulting components 
eventually settle down into bodies of negligible rotation, so that they retain the 
least possible amount of energy, the internal energies of the pieces will be given 
by the appropriate expression for the spherical form (Chapter IV, equation 6), 
and the energy condition for complete separation of the components is now, 
omitting common factors, 

- 0-745 ^ - 0-967 w 5 / 3 - 0-967 (1 - m) 5/3 , 
or 

and this requires ra< 0-19, giving a limiting mass-ratio of rather more than 4:1. 
But in any actual case the mass-ratio must exceed this value, for the components 
must inevitably possess a certain amount of rotation. 

If the reasonable supposition is made that the angular velocity of the com- 
ponents is the same as that of the critical ellipsoid, but that the separate pieces 
have the form of Maclaurin spheroids, then it is readily found by interpolation 
from Table I (p. 39) that when e = 0-71, the potential and kinetic energies for a 
mass mM are given by 

V = -0-9571 G(mM)*'*p l < B 

T = 0-0868 G(mM)*'*pV* 

in the units of the present discussion. Accordingly, the total energy of 
the component mM when separated to infinite distance from (\-m)M is 
0-8703 G(mM) 5/3 p l/3 t and the energy condition in this case is 



-0-745^ - 
or 0- 856 ^ w 5/3 + ( 1 - w) 5/3 . 

This form of the condition requires w<0-ll about, and the limiting mass-ratio 
is therefore about 8:1. 
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Turning now to the question of the conservation of angular momentum during 
the disruption, this must similarly place a restriction on m. But a slight difficulty 
is met with here in that the orbital angular momentum of the separating pieces 
is determined solely by the transverse component of relative velocity of the 
individual masses, whereas the criterion for escape involves only the total relative 
velocity and its direction is immaterial. If it is assumed that the transverse 
relative velocity of the smaller mass itself involves sufficient energy for escape, 
it is readily found that the parabolic orbital angular momentum of mM and 
(1 m)M, supposed to possess spherical form in making this part of the calcula- 
tion, is given by 



Also, the angular momentum of the critical ellipsoid is 



If now, as in considering the energy balance, it is supposed that mM and (1 m)M 
reach great separation in the form of Maclaurin spheroids with u) 2 /27rGp = 0-142, 
then it is readily found from Table I (p. 39), that in the present units the angular 
momentum of mM is given by 



with an exactly similar expression for (1 m)M. Accordingly the conservation 
of angular momentum requires 



or 

and this requires a value of m at most slightly less than 0-14, in other words, 
a mass-ratio rather greater than 7:1. 

In view of the unavoidable element of approximation involved in the fore- 
going calculations, the general agreement between the two limits, arising separately 
from the energy and angular momentum requirements, gives strong indication 
that a mass-ratio of at least 7 or S to 1 must be the outcome of the instability. 

Consideration of double systems in circular orbits 

The validity or otherwise of the fission mechanism for the formation of double 
systems can be investigated from another standpoint by considering the dynamical 
properties of stable double systems, and it is found that this provides further 
evidence on the question of the mass-ratio and the ultimate result of fission. 

The appropriate dynamical problem has been studied by Darwin, who has 
shown that if two masses of equal uniform density are moving in circular orbits 
about each other so that the whole motion is effectively one of rigid-body rota- 
tion, the forms of the free surfaces of the components are to a high degree of 
approximation representable by ellipsoids. If we adopt Darwin's notation and 
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denote the masses by m and M , then A = m/M will be the mass-ratio. The unit 
of length may also be chosen in such a way that 

4 A 4 1 



so that the total combined mass, without change of density, would occupy a 
sphere of unit radius. The distance between the centres of the two masses is 
denoted by R, and the axes of m by a,6,c, and those of M by A,B,C. Then 
Darwin's calculations show that in relative stable equilibrium the two ellip- 
soids are so oriented that the longest axes, denoted now by c and (7, are in the 
line joining the centres of mass, while the two axes of intermediate length, b 
and B, are also in the plane of the orbital motion. The remaining axes, a and A, 
are accordingly in directions parallel to the axis of rotation of the system as 
a whole. 

The problem studied by Darwin was to find the limiting least value of the 
separation of the components, denoted by R, consistent with their remaining 
stable configurations. In every case it is found, as might be expected, that the 
smaller mass, w, becomes unstable first as the distance apart decreases. The 
actual cases investigated by Darwin were for A = 0-4 to 1-0 at intervals of ()!, 
while the value A = corresponds to Roche's problem of an infinitesimal satellite 
in the presence of a primary affected only by its rotation, which is equal to that 
of the satellite in its orbit. Since the values of A that turn out to be of chief interest 
for our present purpose are slightly less than 0-4, the calculations for A = 0-1, 
0-2, and 0-3 have been carried out to complete Darwin's discussion and generally 
exhibit the continuous character of the series of configurations as A varies. 

Table V shows the values of the dimensions of the system in the configuration 
of limiting stability. 

TABLE V 

m = A37 M = mlX 
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It is obviously not necessary to consider values of A greater than 1, since m 
would then simply become the larger mass and only an interchange of m and M 
would be needed for the purpose of investigating the stability of M for A> 1. 
For A = the mass m is infinitesimal, but its shape, which is what determines its 
stability, is indicated by giving the ratio of its axes in limiting form. 
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By means of Table V we can now draw up a further table of relevant quantities. 
First, the quantity R (c + C) will be the shortest distance between the surfaces 
of the components and is of obvious importance in considering the geometrical 
possibility of the orbital motion. Second, may be tabulated the complementary 
quantity B + (c + C), which measures the overall extent of the binary system, 
in order to compare it with the greatest extension of the critical Jacobian ellipsoid 
having the same volume as the combined volume of the components. The fourth 
and fifth columns in Table VI give the angular momentum h and the angular 
velocity in the form eo 2 /27r(r/> in each of the limiting configurations. Below these 
are given the corresponding quantities for the critical Jacobi ellipsoid, while in 
the last line are given the interpolated dimensions of the double-system having 
the same angular momentum as the critical Jacobi figure, which corresponds to 
the value A = 0-36 approximately. 

TABLE VI 

A R-(c+C) R+(c+C) h a>*j2irGp 

0-0 1-427 3-487 0-110 0-0440 

0-1 0-800 4-130 0-227 0-0444 

0-2 0-712 4-230 0-807 0-0441 

0-3 0-673 4-281 0-363 0-0438 

0-4 0-653 4-315 0-402 0-0435 

0-5 0-636 4-334 0-428 0-0434 

0-6 0-633 4-347 0-446 0-0482 

0-7 0-688 4-356 0-457 0-0428 

0-8 0-643 4-361 0-464 0-0426 

0-9 0-650 4-366 0-468 0-0423 

1-0 0-660 4-368 0-469 0-0420 



twice longest axis = 3 ' 772 ' 389 ' 1420 

A = 0-36 0-660 4-303 0-389 0-0436 

It is seen from Table VI that for all values of A greater than about 0-36 the 
binary system will have greater angular momentum than could be stored by 
their combined mass in the form of a single body just on the verge of rotational 
instability. A larger orbit involving greater stability for the smaller mass would 
require even more angular momentum than the critical configuration, that is, a 
smaller value of A would be involved and hence a higher mass-ratio still. It is 
to be noticed that for no value of A does the limiting position of stability of the 
smaller mass allow the two components to be nearly in contact, and in fact the 
overall extent of the double system (4-303) is considerably greater than the 
overall length of the critical Jacobian figure (3-772), so that there can be no 
possibility of the latter configuration evolving gradually into the former. 

Finally may be noticed the considerable disparity between the angular velocities 
of the Jacobian figure and the critical double system of the same angular mo- 
mentum, or indeed of any more stable double system of the same total mass. 
For the Jacobian figure o> 2 /27r(5 ! p = 0-142, whereas for A = 0-36 the orbital angular 
velocity is given by aj^lZnGp = 0-0436 and is little different from this for any 
other value of A. If the separation of the components exceeded the limiting value 
R, the corresponding value of o> would clearly be still smaller. There is therefore 
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a difference of angular velocity between the critical Jacobian form and the 
double system measured by a factor of about 1-8, the former rotating much the 
faster. Now it is well known that if the angular velocity of two particles in a 
circular orbit is increased by a factor ^2 ( = 1-414), the particles would just escape 
from each other, so that there is further strong evidence here in this factor of 
1-8 that the outcome of fission must be to endow the resulting components with 
a relative velocity sufficient for complete escape from each other, and moreover 
that the necessary energy may well be carried almost entirely by the initial 
transverse component of relative velocity, as assumed above in considering the 
angular momentum balance. 

The formation of satellites 

In accordance with the foregoing dynamical arguments it seems safe to con- 
clude that the smaller portion must escape entirely from the gravitational field 
of the larger remaining mass, and the final relative velocity will be shared between 
the two parts in the ratio of their masses, so that relative to their common centre 
of mass the smaller part carries almost all the relative motion. It is suggested 
that the great planets represent the surviving main portions of such rotational 
disruptions of a number of larger primitive planets. The velocities likely to be 
involved can be estimated in order of magnitude by values for the actual planets. 
At the surface of Jupiter the escape velocity is about 60 km. sec." 1 whereas the 
orbital velocity of Jupiter round the sun is about 13 km. sec." 1 and the escape 
velocity from the sun at that distance about 18-5 km. sec." 1 . Thus a particle 
projected from Jupiter at 63 km. sec." 1 would have sufficient energy to escape not 
only from the field of the planet but also from the solar system. Similarly, for 
Neptune the escape velocity at its surface is 23 km. sec." 1 while the orbital velocity 
is 5-43 km. sec." 1 and the parabolic speed 7-65 km. sec." 1 , so that a particle pro- 
jected from Neptune at 24-3 km. sec." 1 or higher would escape entirely. Thus the 
additional energy required for escape from the solar system as well as from the 
planet is such a small fraction of the energy of parabolic or hyperbolic ejection 
that the smaller component resulting from a rotational break-up will usually 
escape entirely from both the remaining mass and the sun, though there is 
clearly a small probability that this might not happen in any particular case. 
On the other hand, the recoil velocity of the larger component will be so much 
smaller that the disruption is not likely to endow it with a velocity approaching 
a hyperbolic value, at any rate at distances from the sun comparable with those 
of the great planets. 

But considering the disruption itself in more detail, it is not to be expected 
that this would take place simply as a clear division into two separate parts, 
for the material of the body near where the main break occurred would necessarily 
be situated in the neighbourhood of a neutral dynamical point from which the 
main masses would be receding in opposite directions. It seems likely that a 
stream of material would form stretched between the masses and having velocity 
ranging between those of the main masses. Such a stream would not be a stable 
arrangement and as the main masses separated would break up into small masses 
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under its internal gravitation. The ones near the centre of the line of small 
bodies would escape from both the main bodies, while some of those nearer the 
ends might have insufficient speeds to escape the adjacent main mass, but the 
general rotation associated with the original mass would of course be maintained 
throughout and endow each of these small masses with orbital angular momentum 
about it. In this way the larger residual portion might come to possess a system 
of small satellite masses moving about it at distances comparable with the size 
of the body itself and in the same general direction as that determined by the 
rotation of the main body. 

The details of such a process inherently preclude mathematical analysis, for 
the removal of material from one gravitating body by the action of another, 
and with the addition here of a rotational field, necessarily involves motion at 
and near neutral points, and in such circumstances arbitrarily small changes in 
conditions can produce arbitrarily large differences in the resulting motions. 
This mathematical difficulty appears to be one that must always surround any 
process involving important redistribution of matter under gravitational forces, 
and it is possible therefore that many of the key cosmogonical processes must for 
ever remain shrouded in this kind of mathematical uncertainty, and perhaps be 
capable of being established only by means of considerations of an indirect but 
circumstantial character. An instance of the kind of argument meant is pro- 
vided by the differences of composition that exist between the four inner planets, 
and by their far smaller masses than those of the great planets. Both these 
features suggest that these terrestrial planets cannot have formed by the same 
process that produced the great planets, but there is the possibility that they 
originated from the central regions of the stream of material between separating 
components of the rotational break-up of a single primitive planet. Their com- 
bined mass is less than 1 per cent of that of Jupiter and therefore would represent 
an insignificant portion of the material involved in the main features of the 
disruption. Since the material would have already reached the liquid state in 
the primitive planet the difficulty concerning the formation of small bodies 
from gaseous matter at high temperature is avoided. Thus the development of 
the fission process in accordance with dynamical indications leads naturally to 
this possible mechanism for the formation of satellites, a process that itself arises 
naturally as part of the evolution of a primitive planet to a state of rotational 
stability. 
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